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Chapitre 1 

Introduction 



1.1 Rappels sur les matrices 



Soit V un espace vectoriel de dimension finie n sur un corps K (K = M 
ou C). Soit B = {ei,e2, • • • ,e„} une base de V, un vectcur v € V admet une 
representation unique dans la base B : 

n 

v = ^ViCi, 
»=i 

En notation matricielle, le vecteur v s'identifie a un vecteur colonne de K n : 

( v, \ 

\v n J 

on notera par v T et v* les vecteurs lignes : 

v T = {v\v 2 ■■■ v n ), v* = (vi v 2 ■■■ v n ). 



V = 



On definit lc produit scalaire dans le cas 
lc cas K = C) par : 



= E (resp. lc produit hcrmiticn dans 

n 

(u,v) = u T v = ^UjVj, 



(resp. 



(u,v) = u*v = y^^UjVj). 



i=i 



On dit qu'une base B — {e\, e 2 , • • • , e„} est orthonormee si : 

e*ej = Si tj 

ou Si j est le symbole de Kronecker : Sij = 1 si i = j, Sij — si i ^ j. 
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Soient V et W deux espaces vectoriels sur le meme corps K, munis de bases 
{ei, 62, ■ ■ ■ , e„} et {/i, /2, • • • , / m } respectivcment. On rappcllc qu'unc applica- 
tion lineaire C de l'espace vectoriel V dans l'espace vectoriel W est definie d'une 
manicrc unique par l'image de la base de V, c'est-a-dirc 



i=l 



Si on note A la matrice m x n ( m lignes et ra colonnes) definie par : 



A = 



( an ai2 

0.21 0,22 
\ «ml a m2 



Clin \ 
amn / 



on voit que la jeme colonne de la matrice A represente le vecteur Cej et on 
vcrific facilcment que pour tout v e V le vecteur Lv de W est represente par le 
vecteur colonne note Av suivant : 



(Av)i — ( ieme ligne de ^4) v, 1 < i < m. 

On definit la matrice transposee de A, notee A T , (resp. la matrice adjointe de 
A, notee A*) : 

(A )ij = (A)ji = a,ji 

(resp. 

(A*) ij = {Aji i = a]i). 



Dans toute la suite, sauf mention du contraire, les matrices considerees seront 
supposees carrees d'ordre n, c'est-a-dire le nombre de lignes est egal au nombre 
de colonnes qui est egal a n. 

Une matrice A carree est par definition dite inversiblc s'il existe une matrice 
notee A~ x , telle que : 

AA^ 1 = A^ 1 A = (Sij) = I = matrice idcntitc. 

On vcrific facilcment que si A et B sont deux matrices inversibles, alors : 

(BA)- 1 = A^B- 1 , {A T )- 1 = {A- 1 ) 7 : (A*)- 1 = (A- 1 )*. 

Une matrice carree A = (a^) est par definition : 

• symetrique si A est reelle et A T = A, 

• hermitienne si A* = A, 

• orthogonale si A est reelle et A T A = AA T = I, 

• unitaire si A* A = A A* = I, 

• normale si A* A = AA* , 

• diagonale si djj = pour i =/= j, on note A — diag(au), 

• triangulaire superieure (resp. inferieure) si = pour i > j 
(resp. dij — pour i < j), 

• semblable a une matrice B s'il existe une matrice P inversiblc, 
dite matrice de passage, telle que P~ 1 AP = B, 
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• diagonalisable si elle est semblable a une matrice diagonale. 
La trace d'une matrice A — (a,ij) carree d'ordre n est definie par 

n 

tr ( A ) = X! a "' 

et son determinant est defini par 

det(A) = ^2 £<ja C r(i)ia (T (2)2 ■ ■ ■ a a (n)n 

oh Q n est l'ensemble de toutes les permutations de l'ensemble {1,2, •• • ,n} et 
e a designe la signature de a. 

Les valeurs propres Xi(A), 1 < % < n, d'une matrice A carree d'ordre n sont 
les n racines de son polynome caracteristique : 

P A (x) = det(xl - A). 

Le spectre de la matrice A est defini par 

sp(A) = {X 1 (A),X 2 (A),... ,X n (A)}. 

Le rayon spectral d'une matrice A carree d'ordre n est defini par 

p(A)=max{\ Xi(A) |}. 

On rappclle que pour toute valeur propre A € sp(A), il cxistc au moins un 
vecteur v ^ tel que Av — Xv appele vecteur propre. 

THEOREME 1.1.1 Soit A une matrice carree d'ordre n, il existe une matrice 
unitaire U telle que U* AU soit triangulaire. Si de plus les coefficients de la 
matrice A sont reels et ses valeurs propres sont reelles, alors, il existe une 
matrice orthogonale O telle que T AO soit triangulaire. 

DEMONSTRATION. Soient E un espace vectoriel de dimension finie n sur 
le corps C et {ei, e 2l • • • ,e n } une base de E orthonormee au sens du produit 
hermitien : 

e i e j — 1 — i'J ' — n 

Soit A = (a,j) une matrice carree d'ordre n a coefficients complexes, on peut 
toujours l'associer a une application lineairc A: E E relativement a la base 
{ei, e2, • • • , e„} par l'image de la base : 

n 

-4( e j) = a%je%, 1 < j <n 

et si A est une valeur propre de la matrice A, alors, il existe un vecteur non nul 
v e E tel que : 

A(v) = Xv 
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En utilisant l'application lineaire A, la premiere partie du thcoreme consiste a 
prouver l'existence d'une nouvelle base orthonormee (au sens du produit hermi- 
ticn) {/i, / 2 , • • • , /„} telle que : 

VI < j < n, Afj G< h, f 2 , ■■■ ,fj> 

oil, < /i, / 2 , • • • , /j > represente le sous-espace de -E engendre par les vecteurs 
fiifi,-" ! /j sur 1° corps C. Ce qui prouve l'existence d'une matrice T tri- 
angulaire supcricure associee a l'application lincairc A rclativemcnt a la base 
{fit J2, • • • , fn}- D'autre part, la nouvelle base {/1, / 2 , • • • , f n } est orthonormee 
au sens du produit hcrmitien, ce qui prouve que la matrice de passage, qu'on 
note U, est unitaire ( U*U = I ) ct on a : T = U*AU. 

Montrons par recurrence l'existence d'une base orthonormee {/1, / 2 , ••• , fn} 
telle que : 

Vl<j<n, Af j e<f u f2,--- ,fj> 

Pour n — 1 le rcsultat est evident. Supposons qu'il est vrai jusqu'a l'odre to. Soit 
A est une application lineaire de E — > E avec E espace vectoriel de dimension 
n = m + 1 sur le corps C. Soit A G C une valeur propre de la matrice associee a 
l'application lineaire A et v\ £ E un vecteur propre associe a A, done : 

Avi = Xv\ 

Quitte a diviser par v\vi, on peut supposer que v\vi = 1. Soit v 2 ,v 3 , ■ ■ ■ ,v n , 
n — 1 vecteurs de E tels que {v\, V2, ••• ,v n } soit une base orthonormee de E. 
Done 

Avi = Xvi 

Av 2 = X)fe=l a k,2Vk 
Av n = J2k=l a k,nVk 

Soit B l'application lineaire de l'espace vectoriel F =< v 2 ,v 3 , ■ ■ ■ ,v n >, sur le 
corps C, dans lui me me definie par : 

n 

V 2 < j < n, Bvj = ^ o^k,jV k 

k=2 

d'apres l'hypothese de recurrence, il existe une base orthonormee {/ 2 , fs, ■ ■ ■ , f n } 
de F telle que : 

V2< j<n,{ t{ = ^= 2 J'° Vl , 

~ J ~ '\ Bfj G </ 2 ,/3,-" Jj > 

On pose /1 = v\. II est clair que {/1, / 2 , • • • , /„} est une base orthonormee et 
on a : 

Ah = A/i 

= (E"=2 7»,jai,i)/i +&fi 

G < /l,/2,-" ./j > 
ce qui termine la demonstration par recurrence. La mcmc demotranstion reste 
valable si on remplace le corps C par le corps R et on suppose que les valeurs 
propres de A sont reelles ; dans ce cas tous les coefficients seront dans K ce qui 
demontre la dcuxicmc partie du theoreme. 
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COROLLAIRE 1.1.1 1) Toute matrice normale est diagonalisable et admet 
une base orthonormee (pour le produit hermitien) de vecteurs propres. En par- 
ticulier, les matrices unitaires, orthogonales, hermitiennes et symetriques sont 
diagonalisables. 

2) Les valeurs propres d'une matrice hermitienne ou symetrique sont reelles. 

3) Toute matrice symetrique admet une base reelle orthonormee de vecteurs 
propres, c'est-a-dire il existe une matrice orthogonale O telle que O 7 AO soit 
diagonale. 

DEMONSTRATION. 



1/ D'apres le theoreme 1.1.1 il existe une matrice U unitaire telle que : 
U* AU = T = (tij) = une matrice triangulairc. Si on note par p%, p2, ■ ■ ■ ,p n , 
les colonnes de U, alors, p*, p\, ■ ■ ■ ,p* sont les lignes de U* . Par consequent, U 
est une matrice unitaire ( U*U — I ) est equivalent a 

Vi = l,n, Vj = l,n, (jPi)*(pj) = 5y 

On va demontrer que T est une matrice diagonale lorsquc la matrice A est 
normale. Montrons d'abord que T est normale : 

t* = (u*Auy = u*A*{u*y = u*a*u 

T*T = U* A*UU* AU = U'A'AU = U*AA*U = TT* 

ce qui prouve que T est normale. Pour demontrer que T est diagonale, on va 
comparer les coefficients de T*T et TT* d'indice 11 : 



(T*T) n = J2taki - |*n | 2 

i=i 



' Ln| 2 



(TT*)u = ^t M tu = |*ii | 2 + |*i2 1 2 + • ■ • + |*i 
i=i 

ce qui donne : *i2 = *i3 = • • • = *i« = 0. On refait la meme chose avec les 
coefficients d'indice 22 pour demontrer que les coefficients de la deuxiemc ligne 
de T en dehors de la diagonale sont nuls, puis avec les coefficients d'indice 33, 
etc ... 

2/ On a : 



U*AU = D = diag(\ 1 ,\ 2 , ■■■An) 
comme A est hermitienne, alors : 

D* = diag(\i,\2, ■■' An) 
= (U*AU)* = U*A*U 
= U*AU = D 
= diag(\i,\2, ■■ • An) 

c.q.f.d. 

3/ On utilise la deuxieme partie du theoreme 1.1.1 
Exercice d'application : 



a) Montrer que si une matrice carree A est triangulaire et normale, alors, 
elle est diagonale. 

b) Montrer les relations 

(i) de£(matrice triangulaire) — 011022 • • • s„„, 

(ii) sp(matrice triangulaire) = {an, i = l,n}, 

(iii) tr(A + B) = tr(A) + tr(B), tr(AB) = tr(BA), tr(A) = 

ELiMA), 

(iv) det(A T ) = det(A) = X 1 (A)X 2 (A) ■ ■ ■ X n (A), det(AB) = det(BA), 

(v) sp(A T ) = sp(A), sp(AB) = sp(BA), 

(vi) k e N, sp(A k ) = {(X i (A)) k 1 i = l,n}, (utiliser le theoreme 



1.1.11 



Reponse 



a) Voir la demonstration du corollaire 1.1.1 

b) (i) Supposons que A = (a^) est une matrice triangulaire superieure, c'est a 
dire : a^j = pour i > j, (meme raisonnement dans le cas triangulaire inferieure, 
en utilisant (iv)). Soit a une permutation. Si le produit <V(i)ia CT (2)2 • • • <V(n)n 
est non nul, alors, necessairement a(i) < i pour tout i = l,n . D'oii : 

cr(l)<l =*ct(1)=1 
cr(2) ^ <x(l) = 1; <t(2) < 2 => a (2) = 2 



cr(n) a(i) = i, 1 < i < n — 1; cr(n) < n =>■ cr(n) = n 

done, toutcs les permutations donncnt un produit nul, sauf peut-etre l'idcntitc. 
Ce qui donne : 

det(A) = ana 2 2 • ■ ■ a nn 

(ii) Si A est une matrice triangulaire, alors, A — xl est aussi triangulaire et 
les coefficients de la diagonale sont : an — x , i = 1, n . D'apres (i), le polynomc 
caracteristique de A est egal a : 

Pa(x) = (an - x)(a 2 2 - x) ■ ■ ■ (a nn ~ x) 

c.q.f.d. 

(iii) (iv)(v) A = (a^), B = (bij), AB = (<kj), BA = (dij), par definition du 
produit de deux matrices, on a : 

n n 

c-ij = ^ OifeOfcj, dij = 2^ bikakj 

k=l k=l 



d'ou 



tr{A + B) =ELiK» + M 

= Sfc=l Qii + Sfc=l b% 

= tr(A)+tr(B) 
tr(AB) =EILi^ 

En tt-~ \7i 1 
i=l l^k=\ a ikOki 

En r^il 1 
k=l 2^ii=l °kia%k 

= Efc=l d kk 

= tr(BA) 
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D'apres le theoreme |1.1.1| toute matrice est semblable a une matrice triangu- 
laire T. D'apres (ii), les valeurs propres d'une matrice triangulaire T sont les 
coefficients de la diagonale de T. Montrons que deux matrices semblables ont les 
memes valeurs propres : soit A € sp(P~ 1 AP), par definition, il cxiste v ^ tcl 
que P _1 j4Pi; = Xv, d'ou APv = XPv et comme P est inversible, alors Pv ^ 0, 
par suite Pv est un vecteur propre de A associe a A, ce qui montre que A € sp(A). 
De la me me maniere, on montre que les valeurs propres de A sont des valeurs 
propres de P~ Y AP. Done, la matrice A et la matrice triangulaire associee a A 
par le theoreme 1 1 . 1 . 1 1 ont le meme polynome caracteristique puisqu'elles ont les 
memes valeurs propres. D'autre part, un calcul simple nous montre que : 

P B (x) = det(B - xl) = (-l) n [x" - tr(B)x n ~ 1 + ■■■ + (-l) n det(B)] 

ce qui montre que tr(A) — tr(T) et det(A) = det(T) , soit encore : 

HA) =E? =1 A,(A) 
det(A) =n? =1 Ai(A) 

Montrons que : det(A) = det(A T ) . II suffit de remarquer que si une matrice 
B est semblable a une matrice C , alors, B T est semblable a C T . Done, A T 
est semblable a T T qui est une matrice triangulaire inferieure ( ses valeurs 
propres sont sur sa diagonale, meme raisonnement que le cas d'une matrice 
triangulaire superieure). Or, la diagonale de T est la meme que la diagonale 
de T T , done, T et T T ont les memes valeurs propres. Par consequent, A T 
et A ont les memes valeurs propres. Ce qui prouve que : det(A) — det(A T ) et 
sp(A) = sp(A T ). Montrons que : sp(AB) = sp(BA) (ce qui prouve en particulier 
que det(AB) — det(BA) ). II suffit de prouver que sp(AB) C sp(BA) et par 
symetrie on a l'autre inclusion. Soit A € sp(AB), alors il existe v ^ tcl que 
ABv = Xv : 

1. Si Bv , alors, BA(Bv) = XBv (on applique B de deux cotes) d'ou Bv 
est un vecteur propre de BA associe a la valeur propre A , ce qui montre 
que A S sp(BA). 

2. Si Bv — , alors necessairement A — 0. D'autre part, si A est inversible, il 
existe w ^ tel que Aw = v, sinon, il existe w ^ tel que Aw = 0. Dans 
les deux cas BAw = 0, ce qui montre que w est un vecteur propre de BA 
associe a la valeur propre A = 0, d'ou A — G sp(BA). 

(vi) On verifie facilement que si B est semblable a C alors B k est semblable a 
C k , k S N. D'o u, A k est semblable a T k (T est la matrice triangulaire donncc 



par le theoreme 1.1.1 1. Un calcul simple nous montre que T , k € N, est aussi 
triangulaire superieure et les coefficients de sa diagonale sont les t k i; i — 1, n. 

REMARQUE 1.1.1 D'apres le corollaire \l . 1 . 1\ si A est une matrice hermi- 
tienne, il existe une matrice unitaire P = U* = U^ 1 telle que : 

A= P* diag(Xi) P 

d'ou, pour tout v <EV : 



■j*Av = v*P* diag(X t ) Pv = (Pv)* diag(XA (Pv) = J] X l \v 



i=i 
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avec Vi la ieme composante du vecteur Pv, d'ou : 

n n 
X max \\Pv\\l = X m ax^2\Vi\ 2 > V* Av > A TO j„||Pv||2 = X min ^ 
i=l i=l 



Vi\ 2 



avec \ m ax la plus grande valeur propre de A et X m i n la plus petite valeur propre 
de A. D' autre part, la matrice P est unitaire ce qui donne ||-Pu||2 = 1Mb, d'oii : 

n n 

X m ax | M 1 2 = ^max ^ ' \ v i\ 2 v -^V > \ m in\ Ml 2 = X m i n ^ ' \ v i\ 2 
i=l i=l 

ce qui prouve en particulier que : 

Xmax = max u*Au; X min = min u* Au 

ll«l|2=l ||u|| 2 = i 

he maximum est atteint pour un vecteur propre de A associe a X max de norme 2 
egale a 1 et le minimum est atteint pour un vecteur propre de A associe a X min 
de norme 2 egale a 1. 

DEFINITION et REMARQUE 

On dit qu'une matrice hermitienne est definie positive (resp. positive) si : 

v*Av > 0,Vw e V - {0}, (resp.v*Av > 0,Vt; e V) ; 

d'apres ce qui precede, on a 

v* Av > \ m i n v*v, pour tout v e V ; 

d'ou, une matrice hermitienne est definie positive (resp. positive) si et seulement 
si X m in > (rcsp.A mm > ). 

1.2 Normes et suites de matrices 

Soit V un espace vectoriel de dimension finie. On rappcllc les trois normes 
vectorielles suivantes : 

IMIi = Ei^i 

ii ^ ib = d>,i 2 ) i/2 

i 

II v Hoc = max I Vi I 



DEFINITION 1.2.1 Une norme matricielle, notee comme la norme vecto- 
rielle par \\ . \\, est par definition une norme vectorielle : 

(i) || A ||> et || A ||=0<^=> A = 0, 

(ii) || 7 A || = | 7 HI A || pour tout scalaire 7, 
(m)\\A + B\\<\\A\\ + \\B\\, 

de plus elle verifie : 

(iv)\\ AB\\<\\A\\\\B\\. 
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Etant donne une norme vectorielle |j . ||. on lui associe une normc matricicllc, 
appclee normc matricicllc subordonncc, dc la manicre suivante : 

II A II 

A = max ■ 



x^O || x | 

(noter que les normes a droite sont vectorielles) . 

Exercice d'application : Verifier que 

(i) la norme subordonnee est bien une norme matricicllc, 

(ii) || A ||= max —. n- = max II Au | 

x#0 || X || ||u||=l 

Reponse : 



|L4||=0 ^sup^ ^ = ^V.x^0,^ = 



1. (a) 

|L4||=0 ^ 

Vx, \\Ax\\ = ^^ = 
(b) \\XA\\ = sup^ = sup^ = |A|sup^ ^=\X\\\A\\ 



(c) Act B deux matrices de meme ordre : 

W T / n ll( A +- B ) a; ll - l|Ax+Bx|| < IjAxjl+jlgxjj 

V:E ^ U ' INI ~ llxll - llxll 



<^ + ^f <W + ||B|| 

d'ou :\\A + B\\= su P;c/0 M < || A|| + ||B||. 

(d) De la definition dc ||A|| = sup^g , on dcduit que : \\Ax\\ < 
||A||||x| pour tout x. Soient A et B deux matrices de meme ordre, 
on a : 

\\ABx\\ \\A{Bx)\\ \\A\\\\Bx\\ \\A\\\\B\\\\x\\ 

V ^°' INI INI " INI " INI ~ mm 

ce qui prouve que : ||AB|| < ||A||||S||. 

2. 

\\A\\ =sup^ ^=sup^ P(^)|| 
< sup||„|| =1 \\Au\\ = su P | H | =1 ^ 
<su Pu/0 ^ = ||A|| 

Les normes matricielles subordonnees aux normes || . ||i, || . || 2 et || . ||oo 
sont donnees par lc thcorcmc suivant : 

THEOREME 1.2.1 (1) || A ||i= max || a,- ||i, (a, : jeme colonne de A), 

(2) M|| 2 = y/rt&Aj, 

(3) || A \\oo= max || a- ( o! i : ieme ligne de A), 

i 

DEMONSTRATION. 
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(1) u G R", || u || 1= 1, alors 



II AU ||l = I X ak <3 U 3 I 

fe=l,n j= l,n 

k — l,n j — l,n 

j — l.n A; — l,n 

< max( ^ | a fcJ |) ^ | Uj | 

fc— l.n j — l z n 

ce qui prouve que 

|| A ||i< max || Oj ||i 

j 

D' autre part, soit jo tel que : 

max || aj || i = || a jo ||i 

on pose 

Vj = 1, = 0, V j ^ j 

on obtient que || v ||i= 1 ct || Av ||i= max^ || ||i<|| A ||i, d'ou l'autre 
inegalite. 
(2) 

|| A \U = max || Au || 2 



max V u* A* Au 

Il«ll2 = l 



car la matrice A*^4 est hermiticnnc positive (voir remarque 1.1.1 1. 
(3) u G K n , || u 11^= 1, alors : 



1 1 co = max | a k,jUj 

j=l,n 

< max | a^ jUj 



k 

j = l,n 



II u | |oo max 



fc 

ce qui prouve que 

II A ||oo< max || a*. ||i 
D' autre part, soit fco tel que : 

max || a' k ||i = || a' ko ||i 

on pose 

signe(a ko j) , si a feoJ - ^ 0, 
, sinon 



on obtient que || v \\oc— 1 ct || Av ||oo= maxfe || a' fe ||i<|| A Hoc, d'ou l'autre 
inegalite. 
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On verifie sans peine que la norme matricicllc subordonncc a la normc 
2 d'unc matricc normalc est egale a son rayon spectral (Indication : utiliser le 
corollaire 1.1.1 et le fait que la norme 2 est invariante par transformation uni- 
taire). Une question s'impose, le rayon spectral dcfmit-il une norme matricicllc? 
La reponse est donnee par l'exemple suivant : 



.4 



1 




on a que A est une matrice non nulle et son rayon spectral est nul. 

Le theoreme suivant compare le rayon spectral et unc norme matricicllc 
quclconquc, 

THEOREME 1.2.2 (1) Soit A une matrice quelconque et || . || une norme 
matricielle quelconque. Alors 

P{A)<\\A\\ . 

(2) Etant donne une matrice A et e > 0, il existe au moins une norme matri- 
cielle subordonnee telle que 

\\A\\<p(A) + e. 

DEMONSTRATION. 

(1) Soit A une valeur propre de A telle que p(A) =\ X |. Soit v ^ un 
vecteur propre de A associe a A. Soit B la matrice carree d'ordre n definie par : 



On a : 

D' autre part 



Bx = (v*x)v, Vx e C" 
AB ||<|| A mi B ||, et || B \\^ 



ABx = A(v*x)v = (v*x)\ Vl Vx e C" 
par consequent AB = XB, d'ou || AB || = | A ||| B ||. Commc || B ||^ 0, alors 

p(A)=\X\<\\A\\ 

(2) D'apres le theoreme 1.1. 1| il existe une matrice unitaire U telle que 

T = XJ- X AXJ 
est une matrice triangulaire superieure : 







tin \ 



V o • ■ ■ o t nn J 

avec sp(A) = {tn,t 2 2, t 33} ■ ■ ■ , t nn }. Soit S > 0, on pose 

P s = diag{l,5,5 2 ,--- ,S n ~ 1 } 
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On verifie facilcmcnt que 

\ in l p i lTP s = diag{t n ,t 2 2, • ■ ■ ,W 



<5->0 

d'ou 

lim || P^U-'AUPs || 2 =|l diag{t n ,t 2 2,--- ,t nn } h= p(A) 
Soit e > 0, on choisit 5 > tel que : 

|| P^U-'AUPs h<p{A) + e 
et on definit l'application de l'espace des matrices carrees d'ordre n dans M + 
5 HI B || = || P- X U- X BUP S ||a 

on verifie facilement que . est une norme matricielle subordonnee a la norme 
vectorielle suivante : 

v HI v ||=|| PT 1 ^ || 2 

La notion de convergence d'une suite de matrices n'est autre que la notion 
classique de convergence dans les espaces vectoriels normes. 

THEOREME 1.2.3 Soit B une matrice carree. Les conditions suivantes sont 
equivalentes : 

(1 ) lim B k = 0, 

k — >OQ 

(2) lim B k v = 0, pour tout vecteur v } 

k — >OQ 

(3) p(B) < 1, 

(4) || B ||< 1 pour au moins une norme matricielle subordonnee. 

DEMONSTRATION. 

(1) ==> (2) evidente. 

(2) => (3) Soit A une valeur propre de B telle que p(B) —\ X |. Soit u^O 
un vecteur propre de B associe a A. On a : 

B k v = X k v 

d'ou 

|| B k v ||= (p(B)) k || « || 

ce qui prouve que lim^oo p(B) k = car || v || ^ 0, d 'ou, p(B) < 1 

(3) => (4) on applique (2) du thcoremc 



avec e = — i^— 1 



1.2.2 

(4) (1) evidente en utilisant la propnete de la norme matricielle 
suivante : 

II B k \\<\\ B \\ k 

THEOREME 1.2.4 Soit B une matrice carree, et \\ . \\ une norme matricielle 
quelconque. Alors 

lim || B k \\V k = p(B). 

k — >oo 
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DEMONSTRATION. 

En utilisant le theoreme on verifie facilement que 



d'ou 



sp(B k ) = {X k ■ A G sp(B)} 
p(B k ) = (p(B)) k 



d'apres le theoreme 1.2.2 



on a 



p(B k ) <\\ B k 



d'ou 



p{B) <|| B k \\ 1,k 

Soit e > 0, d'apres (2) du theoreme 1.2.2| il existe une norme matriciellc || . || e 
telle que : 

\\B\\ e <p(B)+e 

D'autre part, toutes les normes sont equivalentes car l'espace vectoriel des ma- 
trices est un espace vectoriel de dimension finie et la norme matricielle et en 
particulier une norme vectorielle, done : 

3c £ >0 ; || A ||< c s || A \\ e , VA 

Par consequent : 

p{B) <\\ B k || 1 / fe < c i/ fc || B k \\ 1 J k <c 1 J k {p{B) + e) 
Par passage a la limite sur k puis sur s, on obtient ce qu'il faut. 



Une importante norme matricielle non subordonnee a une norme vecto- 
rielle est la norme de Frobenius definie pour toute matrice A G Al m „(K) par : 

m n 

miif= (EEi^i 2 ) 172 ' 

i=l 3=1 

on verifie que 

|| A \\p= trace(A*A) 

ce qui prouve que la norme de Frobenius est invariante par transformation uni- 
taire. La norme de Frobenius est essentiellement la norme Euclidienne appliquee 
a un vecteur de mn composantes. II est facile de voir que la norme de l'identite 
I est toujours egale a 1 pour une norme subordonnee. D'autre part, il est clair 
que || / || f— \fn. D'ou, pour n > 2, la norme de Frobenius ne peut pas etre 
subordonnee a une norme vectorielle. 



1.3 Methode de Gauss pour la resolution des 
systemes lineaires 

L'idee de base derriere les methodes de resolution de Ax = b est la 
transformation de ce probleme en un probleme facile a resoudre. Considerons 
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l'exemple d'un systeme en trois dimensions 





xi + x 2 + 2x 3 = 3 

2x\ + 3x2 + x 3 = 2 

3a;! — x 2 — x 3 =6 

ce qui corresponds a 



A 



L'inconnu x x peut-etre eliminer de la deuxieme equation et de la troisieme 
equation en retranchant de la deuxieme equation deux fois la premiere equation 
et de la troisieme equation trois fois la premiere equation, on obtient alors, le 
systeme suivant 

x\ + x 2 + 2x 3 = 3 

x 2 - 3x 3 = -4 . 
- 4x 2 - 7x 3 = -3 

(II est clair que ces transformations ne changent pas la solution.) 

L'inconnu x 2 pcut-ctrc eliminer de la troisieme equation du systeme en 
ajoutant a la derniere equation quatre fois la deuxieme equation, on obtient, lc 
systeme triangulaire suivant : 

x\ + x 2 + 2x 3 = 3 

x 2 - 3x 3 = -4 . 
- I9x 3 = -19 

La solution de notre systeme de depart peut-etre determinee directement 
a partir des equations du dernier systeme. Dans la derniere equation il y a 
seulement x 3 , et on a X3 = 1. En rcmplacant x 3 par sa valeur dans la deuxieme 
equation, on obtient 

x 2 = -1. 

Enfin, en remplagant x 3 et x 2 dans la premiere equation, on obtient 

x x = 2 

soit 




On a done transforme notre systeme en un systeme qui a la mcmc solution 
et qui est facile a resoudre, e'est le systeme triangulaire. 

1.3.1 Systeme lineaire triangulaire. Methode de remontee. 

Supposons que nous voulons resoudre 

Ux = L 
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ou U est une matrice triangulaire superieure n x n invcrsiblc. On a alors, n 
equations sous la forme : 

unxi + ui 2 x 2 + ■■ ■ + u\ n x n = bi 
U22X2 + •■■ + u 2n x„ = b 2 

( La matrice U est inversible si et seulement si les elements de la diagonale sont 
non nuls.) 

La n eme equation depend uniqucmcnt de l'inconnu i„, on a 

_ , .. _ K 

^nn^n — Oni SOIL X n 

linn 

La (n — l) eme equation 

^71— 1,71— 1^71— 1 ^71—1,71^71 bn — 1 

depend uniqucment de x n et x n -i, or, x n est connu, d'ou 

^Tl—l (Pn—1 u n—l,n%n)- 

^71—1,71—1 

Pour k > 0, Xk est determine de la meme maniere que les inconnus 

^715^71—1)'*' l%k+l 

par 

Xk = (bk - Uk.k+l^k+l - U k _ k+2 X k+2 - ■ ■ ■ - U kin X n ). 

U kk 

Pour tout k = I,-- - ,n, le calcul de x k necessite une division, n — k 
additions et n — k multiplications. D'ou, lc nombrc necessaire d'operations 
clcmcntaircs pour resoudre par la methode dc remontee un systcmc triangu- 
laire est : 

n divisions 

(n - 1) + (n - 2) + ■ ■ ■ + 2 + 1 = rj ^ r 11 w ^ additions 

(n - 1) + (n - 2) + ■ ■ ■ + 2 + 1 = ^^^^ w %- multiplications 

Dans le cas d'un systeme lincaire triangulaire inferieure 

Lx = b, 

on utilise les memes techniques de la methode de remontee, au lieu de commencer 
par l'inconnu x n et on monte a x\, on commence par x\ puis on descend a x n . 
On appelle cettc procedure la methode dc descente. 

Exercice d'application 

Montrer en utilisant la methode de remontee que l'inverse d'une matrice 
triangulaire T = (tij) inversible est une matrice triangulaire de meme nature. 
De plus, les elements de la diagonale de la matrice inverse sont les inverses des 
elements de la diagonale de la matrice T. (Indication : la jieme eolonne de la 
matrice T~ x est egale a la solution du systeme triangulaire Tx = ej avec ej est 
le jieme vecteur de la base canonique.) 
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1.3.2 La methode d'elimination de Gauss 

La resolution d'un systeme lineaire triangulaire est facile. Exactcmcnt 
comme le travail a la main de l'exemple, on va transformer le systeme lineaire 
Ax = b en un systeme lineaire triangulaire ayant la meme solution. 

Soit A une matrice carree d'ordre n 





( 


an 


aii ■ 


■ a ln ^ 




f 


a'l \ 


A = 




an 


a 2 2 ■ 












\ 


a n i 


a n 2 ■ 


a nn ) 




V 


a 'n J 



ou aij e R pour 1 < i,j < n et a! i represente la i eme ligne de la matrice A, 
l<i<n. 

La premiere etape de la resolution de Ax = b consiste a elimincr xi de 
toutes les equations sauf la premiere. Dans le cas ou an ^ 0, on applique la 
technique de Texemplc, c'est a dire : 

la ligne a' 2 est remplacee par a! 2 — ^ a i 
la ligne a' z est remplacee par a' 3 — ^a[ 

la ligne a' n est remplacee par a' n — ^a[ 

de meme pour le vcctcur b : 

la composante b 2 est remplacee par b 2 — ^761 
la composante b 3 est remplacee par 63 — 

la composante b n est remplacee par b n — ^bi 



L 'clement an s'appelle le pivot. Si an = 0, on cherche un coefficient 
non nul an, i = 2, • • • , n, (un tel coefficient existe, sinon la matrice est non 
inversible) et on permute la ligne i avec la premiere ligne pour que le nouveau 
pivot qui est le coefficient a la position 1, 1 soit non nul. 

A Fetape k, la matrice A^ a la forme suivante : 



/ an 




A (k) = 



\ 



ai2 
"22 







aik 

a 2k 



(k) 
1 kk 



(fe) 
2 nk 



ain 
(2) 
2)1 



\ 



a 



(fe) 

ann / 



( 



,(2)' 



V a^' ) 



L'etape k de la resolution consiste a eliminer l'inconnu Xk de toutes les equations 
sauf les fc-prcmieres. De la meme maniere que la premiere etape, si le pivot 
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a< kk 7^ 0, on fait les transformations suivantes : 

la ligne est remplacee par — ak f k ) k 

Ik) 

la ligne aj^_ 2 est remplacee par ajK 2 — "^jwr ^ 

a kk 

la ligne est remplacee par — "^a^ 

a lk 

de meme pour le vecteur 

(fc) 

la composante est remplacee par ^wr^k^ 

a kk 
(fc) 

la composante es ^ remplacee par &j^ 2 — °' k ^ ) k 

a kk 

la composante b^ est remplacee par b^ — -f^b^ 

a lk 

Si lc coefficient = 0, on chcrche un coefficient ^ 0, i = k + 1, ■ ■ ■ ,n, (un 
tel coefficient existe sinon la matrice est non inversible) et on permute la ligne 
i avec la ligne fc pour que le nouveau pivot qui est le coefficient a la position kk 
soit non nul. 

Au bout dc n — 1 etapes, on obticnt un systcme triangulaire. 



1.4 Calcul de l'inverse d'une matrice 

Dans la pratique, on evite le calcul de l'inverse A" 1 d'une matrice inver- 
sible A. Dans le cas particulicr ou on a vraimcnt besoin de l'expression dc la 
matrice A^ 1 , on utilise l'algorithme dc Gauss Jordan. Le principc dc la methode 
de Gauss- Jordan est le meme que celui de la methode de Gauss. On initialise 
une matrice B a l'identite, cette matrice va jouer lc role du second membre de 
la methode de Gauss. 

La premiere etape de la methode de Gauss- Jordan est la meme que celle 
de la methode de Gauss, de plus, on applique les memos transformations a la 
matrice B. Supposons que le resultat de la (fc — l)-eme etape est : 



A (k) 



(fc) 
li 



V 



(fc) 

l fc-l,fc-l 



(fc) 



(fc) 



(fc) 

l nk 



(fc) 
l ln 



(fc) 

(Inn 



( *s> 



B (k) 



V K 



b (k) 6 (fc) 



(fc) Jfc) 



J n1 



On cherche un pivot non nul (par exemple par la strategie du pivot partiel) : 



(fc) 
1 ik 



I 



(fc) I I (fc) 

k<i<n 
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Si k ^ I, on cchange la ligne k avec la ligne I dans les deux matrices et 
B^ k \ On notera encore par A^ la matrice obtenue apres permutation de la 
matrice A^ k \ Puis, on fait l'elimination dans les deux matrices A^ et de 
la maniere suivante : 

(la ligne i) est remplacee par 

{(la ligne i) — , p ^ ot (la ligne de pivot)} 
Xi = lc coefficient a la position (i, k) de la matrice A^ 

avec i allant de 1 a k — 1 et de k + 1 a n. Au bout de n etapes, on obtient : 

A™ = diag(a$) flW 

Pour finir, on divise les lignes de k = 1, ■ ■ ■ , n, par a^ k , on obtient alors : a 
gauche la matrice identite et a droitc la matrice A^ 1 . 

1.5 Conditionnement d'un systeme lineaire 

Soit A une matrice carree d'ordre n inversible et b un vecteur de R n . On 
considere le systeme lineaire 

Au = 6, 

de solution exacte et unique u — A~ x b. 

Dans un premier cas, supposons que le second membre du systeme Au = b 
est perturbe en b + 5b et A restant inchangee. Soit u + Su la solution exacte du 
systeme perturbe 

A(u + 5u) = b+Sb. 
Comme Au = b, la relation precedente implique que ASu = 5b, d'ou 

5u = A^Sb. 

Utilisant une norme subordonnee, nous obtcnons une majoration de || 5u \\ : 

II Su \\<\\ A- 1 mi 5b\\. 

La relation precedente montre que la norme de la perturbation de la solution 
exacte du systeme Au = b due a une perturbation du second membre 5b est 
au plus egalc a || A^ 1 || multipliee par || 5b ||. D'autre part, la relation Au = b 
implique que 

IIMI<MIIII«II- 

Combinons les deux incgalites precedcntes, nous obtenons l'inegalite importante 
suivante : 

m<uAUA-i um. 



Par consequent, lc rapport d'amplification des erreurs relatives est au plus egal 
a || A mi A- 1 \\. 

EXEMPLE 
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Soit A la matrice 

/ 10 7 8 7 \ 

7 5 6 5 

8 6 10 9 
\ 7 5 9 10 / 

elle est symetrique, son determinant vaut 1 et sa matrice inverse 



/ 


25 


-41 10 


-«\ 




-41 


68 -17 


10 




10 


-17 5 


-3 


V 


-6 


10 -3 


2 J 



On considere le systeme lineaire 

23 
33 
V 31 J 



Au = 



de solution 



/ 1 \ 

1 

1 

V 1 J 



et on considere le systeme perturbe, ou le second membre est legerement modific, 
la matrice A restant inchangee : 

/ 9.2 \ 

, . -12.6 
, de solution 



A(u + 5u) = 



/ 32.1 \ 
22.9 
33.1 
\ 30.9 J 



4.5 

ce qui donne un rapport d 'amplification des erreurs relatives de l'ordrc de 2000 
(avec || . ||i). Thcoriquement, on obtient avec la meme norme un maximum pour 
le rapport d'amplification cgale a 

|| A HiH A- 1 ||i= 136 x 33 = 4488. 



Maintcnant, nous perturbons la matrice A et nous laissons b fixe. Soit 
u + Au la solution exacte du systeme perturbe : 

(A + AA)(u + Au) = b. 

On suppose que AA est assez petit pour que la matrice perturbee reste inver- 
sible. De la meme maniere que le premier cas on montre l'inegalite suivante : 

IIA«|| <M n M -i.. IIAAH 



At 



A 



La quantite || A \\\\ A 1 || apparait a nouveau pour controlcr l'amplifica- 
tion des erreurs relatives. 

EXEMPLE 

Considerons le meme systeme que l'exemple precedent oil, cette fois on 
perturbe la matrice A 

( 10 7 



8,1 7,2 \ 

7,08 5,04 6 5 

8 5,98 9,89 9 

\ 6,99 4,99 9 9,98 j 



(u + Au) = b, de solution 



f-Sl\ 
137 

-34 

V 22 J 
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Nous definissons le nombrc conditionnement de la matrice inversible A 

par : 

cond(A) = \\ A mi A- 1 || . 

Pour n'importe quelle normc matriciclle subordonnee, la norme de la 
matrice identite I est egale a un. D'autre part, / = AA^ 1 ct || AA^ 1 \\<\\ 
A mi A~ x ||, ce qui montre que cond(A) > 1. On dit qu'un systeme est bicn 
conditionne si le conditionnement de sa matrice est de l'ordre de 1 et il est mal 
conditionne si le conditionnement de sa matrice est tres grand par rapport a 1 . 



1.6 Exercices : chapitre 1 



Exercice 1 

1. Soit A € .Mn(R). Montrer que si rang(A) — 1 alors il existe a, « £ I" tels que 
A = uv T . 

2. On considere la matrice elementaire 

E = /„ — a uv T 

(a) Montrer que E est inversible si et seulement si a u T v / 1. 

(b) On suppose que a u T v / 1, montrer que 

E' 1 =I n -f3uv T ou P = 



a u T v — 1 

(c) Determiner les valeurs propres de E. 

(d) En deduire que si M £ M n (R) est inversible et u £ 1", alors 

det(M + uu T ) = (1 + w t M _:l m) det M. 

Exercice 2 

Dans A4„(M,), on considere la matrice tridiagonale suivante : 



A(a,b) = 



(a b 
b '■■ 

'•• 




V 











\ 




b 

a ) 



a,beR 



, n sont des 



Jusqu'a la troisieme question, on supposera que a = et b = 1. 

1. Verifier que : u T k = (sin(^), sin(f^), • • • ,sin(^)),fc = 1,- 
vecteurs propres de ^4(0, 1). Deduire le sp(A(0, 1)). 

2. En deduire le spectre de A(a, b) et des vecteurs propres associes. 

Exercice 3 

Etant donne une norme vectorielle |j.||, la norme matricielle subordonnee as- 
sociee est definie par : 

u^Q \M\ 

(Les normes a droite sont vectorielles). 
Verifier que 
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1. \\A\\ = sup^ii^! \\Au\\ = sup||„n =1 Pu|| 

2. P|| = inf{fc > 0; ||Aw|| < fc||u||,Vit} 

Exercice 4 

Soit A — (oij) G .M n (K). Montrer les proprietes suivantes 

1. La norme ||.|J2 est invariante par transformation unitaire. 

2. Si A est une matrice normale alors || ^4|| 2 = p(A). 

Exercice 5 

On munit K n dc la norme euclidienne notee ||.||2- On munit M n (R) de la norme 
induite (notee aussi ||.||2 ) et le conditionnement associe est note cond2(A). 
Soit A une matrice inversible de 

1. Montrer que si A est normale alors 



cond'2(A) — 



max; I A, (A) I 
mi^ \K{A)\ 



ou les nombres \i(A) designent les valeurs propres de A. 

2. Montrer que cond2(A) = 1 si et seulement si A — aQ, a G K* et Q est une 
matrice orthogonale. 



Exercice 6 

Soit A G Mr 



la matrice definie par 
/ 1 2 



A = 



\ 



2 
1 / 



1. Soit x = (xi ■ ■ ■ x n ) T G R" la solution du systeme lineaire Ax 
(bi ■ ■ ■ b n ) T est un vecteur donne de R". 

(a) Montrer que x k = ^"=0* ''{-^fbk+i, k = l,...,n. 

(b) En deduire A^ 1 et un vecteur v =fc tel que 1 1 ^4 1 1 1 = 

(c) Calculer condoo(A) et condi( J 4). 

2. Soit (ei, • • ■ , e„) la base canonique de E n . 

(a) Calculer 1 1 ^4 1 1 1 2 - En deduire un minorant de 1 1 ^4 1 1 1 2 - 

(b) Calculer || 4e 2 1| 2 • En deduire que cond 2 (A) > 2 n . 



Exercice 7 

Resoudre, par la methode de Gauss, le systeme lineaire AX = b avec : 

2 

V 2 / 



b ou b 



( 2 1 





4 \ 


-4 -2 


3 


-7 


4 1 


-2 


8 


V -3 


-12 


-1 J 



Exercice 8 (Algorithme de la methode de Gauss) 

1. Ecrire un algorithme d'echange de deux vecteurs de R". 

2. Ecrire un algorithme de resolution de Ax = b dans le cas ou la matrice A est 
triangulaire. 
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3. Ecrire l'algorithme de la methode de Gauss pour la resolution de Ax = 6 (sans 
strategie de pivot) ; puis avec pivot partiel. 

Exercice 9 (Methode de Gauss-Jordan) 

Appliquer la methode de Gauss-Jordan pour calculer l'inverse de la matrice 



A = 



{ 10 


i 


4 





\ 


1 


10 


5 


-1 


4 


5 


10 


7 




V o 


-1 


7 


9 


I 



1.7 Corrige des exercices : chapitre 1 

Reponse 1 

1. Si rang(A) = 1, alors, il existe u G E n , u / 0, tel que : 

Vx G R", 3a x G R/Ax = a x u 

en particulier : 

Vi € {1, • • • ,n},3vi e R/Aei = ViU 

oil (ei)i=i, n est la base canonique de M n . D'ou, si on ecrit x = ^2™ =0 Ziei, on 
obtient : 

Ax = ^2i—0 ZiA6i 

= (Er=o^) u 

= (v T x)u 
= (uu T )a; 

ce qui prouve que A = uv T . 

2. On remarque que pour z G R™ : 

Ez = => 2 — a(« T 2)u = =>- 2 G< u > 

d'oii : KerE C< u >— \e sous-espace de W 1 engendre par u. On supposera que 
u 7^ ct v 7^ 0, sinon E = I n . 

(a) _E inveresible <=> fcerS = {0} -<=> Sit / (car fceriS C< if >) 
£u = (1 - a(v T u))u / 1 - a(v T u) / 0. 

(b) 

Mx G R", EE~ 1 x = E{x - (3(v T x)u) 
= Ex- f3{v T x)Eu 

— x — a(v T x)u — (](v T x)(u — ct(v T u)u) 

— x — a(v T x)u — (v T x)/3(l — a(v T u))u 

— x — a(v T x)u + a(v T x)u 

— x 

(c) Soit V =< v > ± =le sous-espace R n orthogonal au sous-espace engendre 
par v. On a : 

Va; eV,Ex = x 

ce qui prouve que 1 est une valeur propre de multiplicity au moins n — 1 = 
dimV. Si a(v T u) / 0, alors, la nieme valeur propre est egale a A„ = 
1 — a(v T u) / 1, car Eu = (1 — a(v T u))u. Si a(v T u) — 0, alors : 

Ex = \„x => x — a(v T x)u = \ n x => (1 — \„)x = a(v T x)u 
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done \ n = 1 ou x est colineaire a u. Comme Eu = u, alors A n = 1. Dans 
les deux cas X n = 1 — a(v T u). D'ou : 

sp(E) = {1, 1 - a(v T u)} 

(d) On ecrit M + uu T = M(I n - (-l)(M~ 1 u)u T ) = MS, avec a = -1 et 
« = M~ 1 u. On utilise les proprietes du determinant : 

det(M£) = det(M) det(S), 
det(-B) = Ai x A 2 x ■ ■ ■ x \ n 

ou Aj , i = 1, n, sont les valeurs propres de _B. On obtient ce qu'il faut. 

Reponse 2 

1. Pour tout k = 1, • • • , n, il faut verifier que 

, . , fcvr . . . 2kir . . . nfc7r ..y 

Ufc = (sm( — ),sm( — ),..., sm( — )) 

est un vecteur propre de A(0, 1). Done, pour tout j = 1, • • • , n, il faut calculer : 



avec a; j = sin(^). 

En utilisant la relation trigonomctriquc 



a — b a + b 
sin(a) + sin(o) =2cos( — - — )sm( — - — ) 



on obtient pour tout j = 1, • • • , n 



Xj+i+Xj-i =sin( u w+ ; )+sin( u n+ ; ) 

= sin(^+ii^) + sin(it^) 

V n+1 ' v n + 1 / 

= 2cos(-5 T )sin(^-) 

\ n + 1 ' ^ n+1 ' 

= 2cos(-^):r 7 - 

d'ou, Mfe est un vecteur propre associe a la valeur propre Afe = 2cos(^j-). 
2. On a A(a, b) = aI„ + bA(0, 1). Si b = 0, A(a, 0) = a/„ cc qui donne sp(A(a, 0)) = 
{a}. Supposons que b 7^ 0, alors : 

★ AG sp(A(a, 6)) => Jv / 0/A(a, b)v = \v 

=> A(0, i)» = V»^x£ s p(^(o, i)) 

★ AG sp(A(0, 1)) 3u / 0/4(0, 1> = \v 

A(o, 6)u = a/„?j + 64(0, l)v = (a + b\)v a + bA G s£>(A(a, &)) 
D'ou, sp(A(o,6)) = {a + 26cos(^ I ),fc = 1, • • ■ ,n}. 

Reponse 3 Pour tout 1 , on a :||Ar|| < || A\\ \\x\\ , d'ou : 
f. 

|| A|| = S up^ ^ = S u P ^ ||A(^)|| 

< SU PIMI=1 \\ AU \\ < SU PIMI<1 \\ AU \\ 

< sup^n^! IIAIHMI = ||A|| 

2. 

inf{fc > 0/Vx, \\Ax\\ < fe||a;||} < ||A|| 
Soit k G]0,+oo[ tel que pour tout x, on a : ||Ac|| < k\\x\\ alors : 

ce qui donne que : ||j4|| < fc. Done : 

\\A\\ < inf{fc > 0/Vx, ||Ac|| < fc||x||} 
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Reponse 4 

1. Soient U et V deux matrices unitaires (U*U = UU* = V*V = VV* = I n ), 

\\U*AV\\ 2 = ^/ p((U*AV)*(tFAV ) 
= sj p{V* A*UU* AV) 
= l/ P (V*(A *A)Vj 
= ,Jp(A*A) = \\A\\ 2 

car V* (A* A)V et {A* A) sont semblables et par consequent elles ont les memes 
valeurs propres. 

2. Si A est une matrice normale, alors A est diagonalisable : D — diag(Xi) — U* AU 
avec U une matrice unitaire. D'ou D* D = diag(\\i\ ) = U* A* AU ', ce qui prouve 
que : 

P|| 2 =y /p(ATA) = y/p( U'A*AU) 

= y/mzxi | Ai| 2 = sj (max, Ai|) 2 
= max, | Ai| = p(A) 

Reponse 5 

1. Si A est une matrice normale, alors A est diagonalisable : 

D = diag(Xi) = U'AU et \\A\\ 2 = max |A,| 

i 

D'ou 

D- 1 = diag(^) = U- 1 A- 1 {U*)- 1 = U* A~ l U 
Ai 

( U- 1 = U* et (U*)- 1 = U car U*U = UU* = I), ce qui prouve : 
||j4. _1 ||2 = p(A~ L ) = max- 



|Aj| min,|A,| 

car A^ 1 est aussi normale (pour le prouver, on utilise : 
(A* A)^ 1 = A~ 1 (A*)~ 1 = A' 1 ^' 1 )*). 
2. On a : || A\\ 2 = >fp{A*A) = ^/p(Ai4*) (car sp(AB) = sp(BA)). D'ou : 

cond 2 (A) = PlhP" 1 ^ 

= ^(AM)p((A»)-M-i) 

orp(B) = max, |A,(B)| etp(-B _1 ) = max, lA^B" 1 )! = max, = min . \\ i{B)r . 

d'oii : 



1 = cond 2 (A) = 



maxj |Aj(^*A)| 
mini |Ai(AM)| 



ce qui donne max, | Aj {A* A) | =mini | Ai (A* A) | , or, la matrice A* A est hermitienne 
definie positive ( x* A* Ax = (Ax)* {Ax) — \\Ax\\ 2 > 0, pour x / 0), done ses 
valeurs propres sont strictement positives, ce qui donne : 

\i(A*A) = M(A*A) = • • ■ = \ n (A*A) =r>0 

Par consequent : A* A = U*(rI n )U = rl n , d'ou : {A=A)*(^A) = I n . Ce qui 
prouve que -^.A — Q est une matrice unitaire. Comme A est reelle, alors, Q est 
reelle, d'ou, Q est orthogonale. 

Reponse 6 
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1. (a) La fc-eme equation du systeme Ax = b, k = 1, • • • , n — 1, s'ecrit : 

x k + 2x k+1 = b k 

recurrence lorsque k decroit de n a 1. Pour k — n : 
x„ = Z)i=o( -2 )* & "+« = bn- Supposons que : 

n — fc — 1 

»=o 

alors : 

x k = b k — 2x k +i = b k + Y^7=o ~ (— 2) I+1 &fe +4+ i 
= bk + E^(-2) j bk +j 

(b) Si on note par A~ x = (/%), alors, la Z-ieme colonne (/3i; /3 2 ; • • • /3 n ;) T de 
A -1 est egale a A~ 1 ei, done e'est la solution du systeme : Ax = e;. D'ou, 
en remplagant le vecteur b par le vecteur e; dans la question precedente, 
on obtient : 

z fc =/? fci = ^(-2)6 fc+i = j fc>z 

Ce qui prouve que : 

Halloo =max fc ^ =1 |/3 H | = max fe JXi 2'~ fc 
= (Er= 1 2 ; )max fc 2-^|^ 
= 2™ - 1 

On prend v = (-1, 1, -1, • ■ ■ , (-1)™), on a H^oc = 1 et : 

p-^lloo =max fc |E; =1 (-l) 3 /3 fcJ | 

= max fc |E; =1 (-l) 3 (-2) 3 - fc | 
= max fc 2- fc |E; =1 (-l) 3 (-2) J | 
= |E^ =1 2 J =2"-l=p- 1 || 0O 

(c) amd^A) = PIUp- 1 ^ = 3(2™ - 1) car PIU = 1 + 2 = 3. 

condi(A) =P||ip- 1 ||i 

= P1|oc||(^- 1 )*||oo 

= 3max fc= i,„Er=il(-2) fc " i | 
= 3(E?=i 2_i ) max fc= i, n 2 k 

= 3( 2 ~ 1 1 :^:r 1 )2" 

= 3(2™ - 1) 

2. (a) p-^ih = ||ELi^ne fc || 2 =V^C = VE fc 4»-*=2"VET4^ > 

2 „-i 

(b) Pe 2 ||2 = ||2ei + e 2 || 2 = ^5 > 2. D'ou : 

cond 2 (A) = PI^P" 1 ^ > Pe 2 || 2 p- 1 e n || 2 > 2.2™" 1 = 2™. 

Reponse 7 Utiliser l'exemple interactif du chapitre 1. 
Reponse 8 

1. La procedure d'echange de deux vecteurs : u = (wiw 2 «3 • • • u n ) T et v = 

(viV 2 «3 ' ' ' Vn) T 
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Debut de la procedure echang(u,v) 
Pour i allant de 1 an faire 

X — Ui 
Ui = Vi 
Vi = x 
Fin faire i 
Fin procedure 

2. La procedure de remontee pour resoudre un systeme lineaire triangulaire : 
A = (a»j) et b = (&») 

Debut de la procedure remonte(A,b,x) 

Xn — b n 

Pour k allant de n-1 a 1 faire 
%k = b k 

Pour i allant de k+1 a n faire 

x k = x k — a k i x Xi 
Fin faire i 
Fin faire k 
Fin procedure 

3. a) Le programme de la methode de Gauss sans strategie de pivot ( dans ce 

on suppose qu'a chaque etape k de la methode d'elimination de Gauss 
le coefficient a la position (k,k) est non nul ) : 
Debut du programme Gaussl 
Pour i allant de 1 an faire 
lire bi 

Pour j allant de 1 an faire 

lire Oij 
Fin faire j 
Fin faire i 

Pour k allant de 1 a n-1 faire 

appel procedure : elimin(A,b,k) 
Fin faire k 

Appel procedure : remonte(A,b,x) 
Pour i allant de 1 a n faire 

affiche Xi 
Fin faire i 
Debut de la procedure remonte(A,b,x) 

Xn — b n 

Pour k allant de n-1 a 1 faire 
Xk = b k 

Pour i allant de k+1 a n faire 

Xk — Xk Oiki X Xi 

Fin faire i 
Fin fair e k 
Fin procedure 

Debut de la procedure elimin(A,b,k) 

pivot = a kk 

Pour i allant de k+1 a n faire 

Pour j allant de k+1 a n faire 

a,ij = Oij — aik x -^^i 
Fin faire j 
Fin faire i 
Fin procedure 
Fin programme Gaussl 
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b) Le programme de la methode de Gauss avec pivot est le meme que le pro- 
gramme Gaussl, sauf que juste avant de faire appel a la procedure d'elimination de 
Gauss : elimin(A,b,k), on fait appel a la procedure : pivot(A,b,k) qui recherche le pivot 
de l'etape k et permute la ligne de pivot avec la ligne k et il faut ajouter la procedure 
pivot(A,b,k) au programme. 

Debut de la procedure pivot(A,b,k) 
pivot = \a k k\ 
I = k 

Pour i allant de k+1 a n faire 

Si |o»fc| > pivot alors 
pivot — \dik\ 
l = i 

Fin si 
Fin faire i 
Si pivot = alors 

affiche : A non inversible 

quitter le programme 
Fin si 

Si I ^ k alors 

Pour j allant de k+1 a n 

X — dkj 

aij = x 
Fin faire j 
x = b k 
bk = bi 
h = x 
Fin si 
Fin procedure 

Reponse 9 Utiliser Pexemple interactif du chapitre 1. 
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Chapitre 2 

Methodes iteratives pour la 
resolution d'un systeme 
lineaire 



Lcs methodes d'elimination ou de factorisation sont utilisees surtout lorsquc 
Fordre de la matrice est petit ( matrice 100 x 100 par exemple ) ou lorsque la 
matrice est pleine (i.e. peii de coefficients nuls). 

Dans la pratique, beaucoup de problcmcs ncccssitcnt la resolution d'un 
systeme Ax = b d'ordre assez important, avec, A une matrice creuse (i.e. beau- 
coup de coefficients nuls) . Des systemes de ce type sont donncs par exemple par 
1' application de la methode des differences finies ou la methode des elements 
finis. 

Pour ce genre de problemes, on utilise lcs methodes iteratives. Etant donnc 
un vecteur initial arbitraire x°, on construit une suite de vecteurs 

T° T 1 ■ ■ ■ T k ■ ■■ 

qui converge vers la solution x du systeme lineaire Ax = b. 

2.1 Construction d'une methode iterative 

On considere le systeme lineaire 

Ax = b (1.1) 

avec A une matrice carree d'ordre n inversible et b un vecteur de R n . Pour toute 
matrice M carree d'ordre n inversible, le systeme (1.1) est equivalent a 

Mx — (M — A)x = b (1.2) 

ou encore, en posant N — M - A, B = M~ 1 N et c = M~ 1 b 

x = Bx + c. (1.3) 
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Ce qui nous permet de dcfinir la methode iterative : 

x° e R", vecteur initial 
x k+1 = Bx k + c 



(1.4) 



Soit x la solution de (1.1), si on note e k = x k — x le fcicme vecteur erreur, 
on obticnt 



x k - x = {Bx k ~ 1 + c) - (Bx + c) = B^" 1 
Be^ 1 = J B fe e° 



x) 



(1.5) 



On dit que la methode iterative (1.4) converge si la suite de vecteurs (e fe ) 
converge vers zero independamment du vecteur initial x°, ce qui est equivalent, 
d'aprcs le theoreme 2.3 du chapitre precedent, a Tune des deux propositions 
equivalentes : 

(1) p(B) < 1 

(2) || B || < 1 pour au moins une norme matricielle. 

REMARQUE 2.1.1 Pour chaque choix de la matrice inversible M, on ob- 
tient une methode iterative. Le meilleur choix de M doit satisfaire les conditions 
suivantes : 

(a) la matrice M est facile a inverser, 

(b) p(B) = p{M^ 1 N) est le plus petit possible. 

Dans la condition (a), on n'a pas besoin de l'inverse de M, il suffit que le calcul 
de x k+1 en fonction de x k , en utilisant (1.2) : 

Mx k+1 = Nx k + 6, 

soit facile. La condition (b) est une consequence du comportement asymptotique 
de F erreur e k ; en effet, 



lim fe — [max e o^ || e k \\ / \\ e° \\] 1/k = linife^oo [max e o^ || B k e° \\ / || e° ||] J 

= Iim fc _ 00 || B k \\ 1 ' k = P {B) 

la methode est done d'autant plus rapide que p(B) est plus petit (voir theoreme 
1.6 du chapitre precedent). 

On considere la decomposition suivante : 

/ an ai2 ■■■ Oi„ \ 



A = 



«2r. 



= D~E-F, 



D = diag(au,a 2 2, ■■■ , a nn ) 



E = — 



( o 

021 
\ 0„i • • • O n ,n-l / 
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On suppose dans toute la suite que D est inversible et on pose 

L = D- X E, U = D^F, 
(1) la methode de Jacobi 

M = D, N = E + F, J = M~ X N = L + U 

lc calcul de x k+1 a partir de x k se fait dircctcmcnt par lc systcme 



onxf +1 = 
a 22 x k+1 



dn.n.X 



fe+1 



-a X2 x% -ai 3 x% 



-a n \x\ —a n2 x 2 



fl nn-l^n-l 



-a Xn x\ +&i 
-a 2n x k +b 2 

+b n 



On remarque que la methode de Jacobi necessite n memoires pour le vecteur 
x k et n memoires pour le vecteur x k+1 . La matrice J s'appcllc la matrice de 
Jacobi. 

(2) la methode de Gauss-Seidel 

M = D — E, N = F, H=M- 1 N=(I-L)- 1 U 
le calcul de x k+l a partir de x k se fait dircctcmcnt par lc systcme 



a 22 x 2 



-a 21 x 1 



-a 12 x% 



-a 13 x% 
-a 23 x% 



a T k+1 — —n 1 r fe+1 -n or fc+1 



-a nT fc+1 



-a ln x n +&i 
-a 2n x k +b 2 



+b r , 



la methode de Gauss-Seidel necessite seulement n memoires, la composante x k+1 
prend la place de x\ qui ne sera pas utilisee pour le calcul de x k +\ , x k +\ , • • • , x k+1 . 

(3) la methode de relaxation 



M 



D -E, N - 



1 - w 



D + F, 



H(lu) = M- 1 N=(I-ujL)- 1 [(1-lu)I + ujU}, 
le calcul de x k+1 a partir de x k se fait directement par le systeme 
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anx k+1 = clux\ -lo{ a lx x\ +a X2 x\ H +a ln x k - b x } 

a 22 x\ +1 = a 22 x\ -u>{ a 21 x k+1 +a 22 x\ H +a 2n x k - b 2 } 

, QnnX,^ &nnX n Ld{ d n \X^ -\- • ~~\~Q"n.n—lX n _i "^dnnX n ^n} 

la methode de relaxation necessitc sculcmcnt n mc moires. Lc paramctrc lo s'ap- 
pelle le parametre de relaxation et pour lo = 1 on retrouve la methode de 
Gauss-Seidel. 

Exemple : 





La matrice de Jacobi J et la matrice de Gauss-Seidel H sont donnees par : 

1 / -1 -1\ / 2 \ _1 / -1 -1 
J=- -1 -1 , if= 1 20 00 -1 

2 \ 2 -2 0/ \ -1 1 1 / \ 

Le systeme qui donne la suite de la methode de Jacobi est donne par : 

2x k+1 = -x\ - x\ + 1 

T fc +l _ _ k _ k 

ZlJb^ 1 3 

™^~t~l ^fk ^k 

x 3 — x l - h 2 

Le systeme qui donne la suite de la methode de Gauss-Seidel est donne par : 

2x k+1 = -x\ - x% + 1 



O^-k+L _ _ rr k+L _ k 

T k + l _ T k+l _ k+1 
x 3 — x l x 2 



Pour les deux mcthodes, on demarre de x° = (0 0) T . A chaque iteration k de 
la methode de Jacobi, le vecteur x — (xi x 2 x 3 ) T joue le role de x k et le vecteur 
V = (j/i U2 Vz) T joue le role de x k+1 . Par exemple, si on veut calculer x 3 : 
k = : x = (0 0) T et y est donne par : 

yi = (-X2 - x 3 + l)/2 =1/2 
2/2 =(-Xi-x 3 )/2 =0 
J/3 = xi - x 2 =0 

1 : x = (1/2 0) T et y est donne par : 

= (-x 2 -x 3 + l)/2 =1/2 
= (-an - x 3 )/2 =-1/4 

= .Tl - x 2 = 1/2 

(1/2 - 1/4 1/2) T et y est donne par : 

(-s 2 - s 3 + l)/2 = (1/4 - 1/2 + l)/2 = 3/8 
(-xi-x 3 )/2 = (-1/2 -1/2) /2 = -1/2 
xi-xa =1/2+1/4 = 3/4 
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Ce qui donne : 




A chaque iteration k de la methode de Gauss-Seidel, le vecteur x = (x\ xi x 3 ) T 
joue le role de x k et au meme temps le role de x k+1 . Par exemple, si on veut 
calculer x 2 : 

k = : x = (0 0) T et x est donne par : 

xi = {-x 2 -x 3 + l)/2 = 1/2 
x 2 = (-x 1 -x 3 )/2 = (-l/2-0)/2 = -1/4 
x 3 = xi - x 2 = 1/2 + 1/4 = 3/4 

= 1 : x = (1/2 — 1/4 3/4) T et x est donne par : 

X! = {-x 2 -x 3 + l)/2 =(l/4-3/4+l)/2 = l/4 
x 2 = (-an - x 3 )/2 = (-1/2 - 3/4)/2 = -5/8 
x 3 =xi-a;2 =1/4 + 5/8 = 7/8 

Ce qui donne : 




On remarque que les iterations de la methode de Jacobi necessitent 6 va- 
riables : le vecteur x et le vecteur y, done 6 memoires et les iterations de la 
methode de Gauss-Seidel necessitent seulement 3 variables : le vecteur x, done 
3 memoires. 



2.2 Convergence 

Dans le cas d'un systemc Ax = b, avec A une matrice hermitienne definie po- 
sitive, le theoreme suivant donne une condition suffisante pour qu'une methode 
iterative associee a une decomposition A = M — N converge. 

THEOREME 2.2.1 Soit A une matrice hermitienne definie positive, decomposee 
sous la forme 

A = M — N, M : matrice inversible. 

Pour que la methode iterative associee a cette decomposition converge, e'est a 
dire 

p{M~ 1 N) < 1 

il suffit que la matrice hermitienne M* + TV soit definie positive. 

DEMONSTRATION : Vcrifions d'abord que M* + N est hermitienne : 

M* + N = A* + N* + N = A + N + N* = M + N*. 
Soit A e sp{M- x N) et soit v ^ tel que 

M~ 1 Nv = \v 
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On pose 



soit encore 



on a alors 



u = M~ x Av = {I- M _1 A/> = v - \v 
Xv = v — u 



A | 2 v*Av = (v — u)*A(v — u) 

= v* Av — v* Au — u* Av + u*Au 

= v* Av — (Av)* u — u* Av + u* Au 

= v*Av-(Mu)*u-u*Mu + u*Au 

= v*Av- u* (M* + M - A)u 

= v*Av -u*(M* + N)u 



d'ou 



(1- I A \ 2 )v*Av = u*(M* + N)u 
comme v ^ 0, u ^ ct A, M* + N sont definies positives, on deduit que 

| A|<1. 

COROLLAIRE 2.2.1 Soit A une matrice hermitienne definie positive, alors 
la methode de relaxation converge pour < uj < 2. En particulier, la methode 
de Gauss-Seidel est convergente lorsque A est hermitienne definie positive. 

DEMONSTRATION : La matrice A est hermitienne definie positive, d'ou, 
les coefficients de la matrice diagonale D sont strictcmcnt positifs ct E* = F. 
D'autre part 

M=—-E, N=^—^D + F 

LO UJ 

d'ou 

2 - u 
M* +N = D 

LO 

done, pour < lo < 2, la matrice hermitienne M* + N est definie positive. 
THEOREME ( Kahan (1958) ) Le rayon spectral de la matrice de relaxation 

H(uj) = (I- ljL)~ 1 {loU + (1 - 

verifie pour tout w^O : 

p(H(u)) >| w-1 I . 

Ce qui montre, en particulier, que la methode de relaxation ne converge pas pour 
wg]0,2[. 

Demonstration : On a 

n w . l i^.ii ^y =l .-r 

or 

inwM)i<p(#Mr 

i 

d'ou 

p(H(to)) >| w-1 | . 



35 



DEFINITION 2.2.1 Soit A une matrice carree et J = L + U la matrice de 
Jacobi associee a A. On pose J{a) = aL + a U. On dit que la matrice A est 
correctement ordonnee si sp(J(a)) est independant de a. Autrement dit 

sp(J(a)) = sp(J), Va^O. 

THEOREME ( Young (1950) ) Soit A une matrice correctement ordonnee 
et u ^ 0, alors 

(a) p G sp(J) =^> — fi G sp(J) 

(b) p € sp(J) et 

(\ + uj-1) 2 = AwV 2 , (2-1) 

alors A G sp(H(u>)). 

(c) ^ A e sp(H(u>))) et jj, verifie (2.1), alors p G sp(J). 

Demonstration : 

(a) J(-l) = -L-U = -J, d'ou 

ap(J) - ap(J(-l)) = sp(-J). 

(b) -(c) Soit A G C*, on a : 

det(XI-H(u)) = det{(I - uL)(XI - H(u))} 

= det{(X + uj-l)I -XluL-luU} 

soient a, [3 G C tcls que 

a 2 = Xuj 2 , et (3 = a/u 

alors, 



dettXI - H(u)) = a n det{^^ — - (/?L + (i^U)}. 

a 

D'autre part, p est une solution de (2.1) si et seulement si 

(A + w-1) 



p = ±- 



Done 



Xesp{H{uj)) _ f dei{^/-(±/3i + (±/3)- 1 C/)} = 

/i = ± (A+a>-l) | M = ± (A+^l) 

J M G sp(J(±/3)) - sp(J) 

^ 1 „2 _ (A+^-l) 2 

COROLLAIRE 2.2.2 £oii A une matrice correctement ordonnee, alors 

p(H) = (p(J)) 2 . 

En particulier, si A est correctement ordonnee la methode de Gauss-Seidel 
converge si et seulement si la methode de Jacobi converge, et si les deux methodes 
convergent et p(J) ^ 0, alors : 

< p(H) < p(J) < 1. 
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DEMONSTRATION : Pour u=l, l'equation (2.1) dcvicnt : 

A 2 = Xp 2 

et d'apres les proprietes b) et c) du theoreme de Young, on a : 

(a) p e sp(J) =>■ A = e sp{H) et A = p 2 e sp(H) 

(b) A e sp(H) =^> ±p e sp(J) avcc p 2 = A 

par consequent, p(J) — si et seulement si p(H) — et si p(J) ^ alors 
p(H) = (p(J)) 2 . 

2.3 Determination theorique du parametre de 
relaxation optimal 

Soit A = D — E — F une matrice correctement ordonnee, avec D inversible. 
Nous supposons que la methode de Jacobi est convergente. D'apres le corollaire 
precedent la methode de relaxation converge pour ui = 1, et par continuitc, il 
y a convergence sur un intervalle contcnant 1. Lc probleme que nous voulons 
resoudre est de determiner lc parametre de relaxation optimal ujb tel que : 

p{H{uj b )) = min p{H{u)). 

UJ 

THEOREME 2.3.1 Soit A une matrice correctement ordonnee, si toutes les 
valeurs propres de J sont reelles et < p{J) < 1, alors, la methode de relaxation 
converge si seulement si < u) < 2 ; de plus, on a : 

(!)<*= ; =i + ( ; (J) V 

(2) p(H(iJ b )) = LJ b - 1 

(3) p(H(u))>p{H(u b )), Vlo^u;,. 

DEMONSTRATION : On peut limiter l'etude a to €]0,2[ (voir theoreme de 
Kahan). D'apres le theoreme de Young, les valeurs propres non nulles de J sont 
des paires ±/ij, avec : 

< m < p{J) 

et a chaque p, £ sp(J) correspond une paire de valeurs propres de H(u>) 

{X m (uj,fi), A m (w,m)} 

solutions de (2.1) (on suppose que | X m (u, p) \<\ Am(w, p) |), et a chaque valeur 
proprc non nulle A G H(w) correspond une valeur propre < p, e sp(J) telle 
que 

A = X m (u), p) ou A = A M (w, p). 

D'ou 

p(H(u)))= max \X M (u,fi)\ 

D'apres l'equation (2.1) 

(u> - l) 2 = X m (ui, p)X M {oj, p) <\ A M (w, p) | 2 
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d'ou 

I Am(w,a«) |>| oj - 1 | 

ct on a : 

| X m (u>,n) |=| A M (w,/i) |=| w - 1 I 

dans le cas d'une racinc double ou dans lc cas de deux racines complexes 
conjuguees. D'apres (2.1), si A € sp(H(u>)) est reelle, alors, necessairement A > 
ct on a, pour u> ^ : 

X + W ~ 1 =±»V\. (3.1) 

On definit 

ai A + w-1 

et 

m M (A) = < it < p(J) < 1. 

Le graphe de g w est la droite qui passe par (1, 1) et de pente ^. Geometriquement, 
les solutions reclles de (2.1) sont donnccs par l'intersection des deux courbes as- 
sociees a g^ et to m : 




fig. 1.5 

Pour oj > u>, les racines de (2.1) sont complexes conjuguees (car la droite 
d'equation g u ,(X) nc coupe pas les courbes ±m Ai (A)). II est clair, d'apres le 
graphe ci-dessus, que pour to ^ et tout V/z £ sp(J) : 



<| \ M {u,p(J)) |< 1, si\ M {u,n)£ 
=1 w — 1 |< 1, sinon 



or 



d'ou 



soit encore 



| A M (w,p(J)) |>| w - 1 | 
V/iG sp(J), | A M (^,A < ) l<l A M (w,p(J)) |< 1 

p(i/H)=|A M (^,p(J))|<l, 
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ce qui prouve que H(ui) converge pour < uj < 2. Encore d'apres le graphe 
ci-dessus, on a : 



| \ M (u,p{J)) \>\ \ M {u,p{J)) I, Vw<u> 
| w — 1 I, > uj 



D'autre part uj est la valeur de uj <G]0, 2[ pour laqucllc l'equation 
(A + uj - l) 2 = Aw 2 (p(J)) 2 A 2 + {2(uj - 1) - w 2 (p(J)) 2 }A + (w - l) 2 = 
a une racine double, c'est a dire : 

{2(uj - 1) - w 2 (p(J)) 2 } 2 - 4(w - l) 2 = 

d'ou 

w 2 (p(J)) 2 -4w + 4 = 
Cette equation admet une racine > 2 et une autre egale a : 



UJ = 



2- y/4-4(p(J))- 

(P(V)) 2 



= 1 + 



l + v /l-(p(./)) 2 



et on a 



D'ou 



^i + v /i-^) 2 ; 

p{H(u>)) =| A M (w,p(J)) |=| w - 1 |= w - 1. 



p(H(ui)); minu; — 1 



min p(H(u>)) = min 

= u) — 1 = p(H(Cj)) 

Par consequent 

CJb = UJ 

vcrifie ce qu'il faut pour lc thcoremc. 

2.4 Exercices : chapitre 2 

Exercice 1 

1. Calculcr dans chacun des cas suivants le rayon spectral de la matrice de la 
methode de Jacobi et le rayon spectral de la matrice de la methode de Gauss- 
Seidel pour la resolution du systeme Ax = b 

/ 1 2 -2 \ / 2 -1 1 \ 

A = j 1 1 1 et A = 2 2 2 
\ 2 2 1 / \ -1 -1 2 / 

Que peut-on deduire ? 
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2. Montrer que si la matrice A est a diagonale dominante stricte : 

VI < i < n, \au\ > ^ |oij| 

alors la methode iterative de Jacobi pour la resolution de Ax = b est convergente 
(on pourra montrer que ||J||oo < 1)- 

3. Etudier la convergence de la methode de relaxation (pour la resolution du 
systeme Ax = b) lorsque 

..(!!!) 

V o o 1 J 

Exercice 2 

Soit A = (aij) une matrice d'ordre n telle que an =fc pour tout i — 1, . . . ,n. On 
rappelle que la matrice de Jacobi J = L + U et la matrice de Gauss-Seidel H = 
(J - Ly 1 !] avec A = D-E- F , L = D~ X E et U = D~ l F. On dit que la matrice A 
est correctement ordonnee si : 

Va / 0,sp(J(a)) = sp(J) 

avec J(a) = ah + q _1 (7. 

1. Calculcr la matrice J = (cij) et la matrice PaJP^ 1 = (dij) en fonction des coeffi- 
cients de A avec P a = diag(l, a, a 2 , . . . , a 11-1 ). 

2. Deduire qu'une matrice tridiagonale est correctement ordonnee. 

3. En utilisant un theoreme du cours, montrer que si A est tridiagonale alors : 

sp(H) = {0} U {iS/li G sp(J)} 



Exercice 3 

Le but de l'exercice est d'etudier le probleme 

f -u"(x) = /(i),i6]0,l[, 
\ w(0) = w(l) = 

Soit h = j^q-j- avec N G N* . Le probleme discretise correspondant est : etant donne le 
vecteur F = (fi)i<i<N trouver U = (ui)i<;<jv tel que 

{-^{ui+i - 2m + m-i) = fi, l<i<N, 
(1) 
wo = itjv+i = 

1. Montrer que le systeme (1) est equivalent a la resolution de 

AU — F 

ou on explicitera la matrice A = (aij)i<i t j<N ■ 

2. Calculer les valcurs propres des matrices de Jacobi et Gauss-Seidel pour la 
resolution de AU — F en fonction de celles de A (On pourra utiliser un exercice 
du chapitre precedent). 

3. Comparer ces deux methodes. 

4. Donner le parametre optimal de relaxation pour la matrice £ w . 
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Exercice 4 
Partie I : 

Dans cette partie, on se propose de prouver quelques proprietes utiles pour la partie 
II. Soient K et M deux matrices carrees d'ordre n. On suppose que K et M sont 
symetriques et (I + K) inversible, I etant la matrice identite d'ordre n. 

1. Montrer que p(K) = \\K\\ 2 . 

2. On rappelle que p(L) < \\L\\ pour toute matrice L et toute norme matricielle 

|.||. En utilisant 1), montrer que : 

p{KM) < p{K)p{M). 

3. Montrer que 

(a) A G sp(K) => X -1 et S sp(K(I + K)' 1 ), 

(b) /3 G sp(K{I + K)- 1 ) => 13 + 1 et ^ G sp(if). 

4. Montrer que : p(K(I + K)^ 1 ) < 1 <4> \l + K est definie positive. 
Partie II : 

Soit A une matrice carree d'ordre n symetrique. Cette partie est consacree a l'etude 
d'une methode iterative de resolution du systeme lineaire Ax — b. On introduit la 
decomposition 

A = D + H + V 

ou D — cl, c > et ou H et V sont deux matrices symetriques telles que les deux 
matrices D + H et D + V soient inversibles. On considere la methode iterative suivante : 



(2.1) 



j {D + H)x (k+ ^ = -Vx {k) +b 
\ (Z? + V> (fc+1) = -Hx {k+ ^+b 

1 . Exprimer x^ en fonction de x' fc '. En deduire 

lim x (k) = x ^ p((D + V)~ 1 H{D + < 1 

k — »oo 

2. On pose B = D^H, C = D~W. 

(a) Montrer que p((D + V)~ 1 H(D + H)~ l V) = p(B(I + B)' 1 C{I + Cy l ) 

(b) Montrer que les matrices B(I + B)^ 1 et C(I + C)" 1 sont des matrices 
symetriques. (On pourra utiliser, apres justification, que B(I + B)^ 1 = 
I-il + B)- 1 ). 

(c) En deduire que p((D + T/)- 1 ff(D + //)- 1 V) < p(B(I + B)- 1 )p(C{I +C)- 1 ) 

3. Montrer que la methode iterative (2.1 1 converge des que \D -+- H et \D + V 
sont definies positives. 



2.5 Corrige des exercices : chapitre 2 




Jacobi :M = D = I,N = E + F = 
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-2 2 \ 
-1 -1 
-2-2 / 
■ 3 ^ S p(J) = {0}^p(J) = 
/ 1 \ 
Gauss-Seidel : M = D - E = 1 1 \,N = 

\ 2 2 1 / 
/ -2 2 \ 
H = M~ 1 N = 2 -3 

V 1 / 
S p( J ff) = {0,2,l}^p(ff) = 2. 

Conclusion : la methode de Jacobi converge et la methode de Gauss-Seidel diverge. 

/ 2 -1 1 \ 
b)A= 2 2 2: 

V -1 -1 2 J 
Jacobi : M = D = 21 , 

I 1 -1 \ / 1 

iV = £; + F= -2 -2 et J = M _1 iV = ± I -2 

V i i o / v 1 1 

det(A7 — J) = A 3 + 5A sp(J) = {0, ±^5} p( J) = y/5 

/ 2 \ 
Gauss-Seidel : M = D — E = J 2 2 , 

V -l -l 2 y 

/ 1 -1 \ / 0.5 -0.5 \ 

N=l -2 et// = M- 1 iV= -0.5 -0.5 

\ / \ -0.5 / 

sp(H) = {0, ±0.5} => p(H) = 0.5. 

Conclusion : la methode de Jacobi diverge et la methode de Gauss-Seidel converge. 
2. La matrice de Jacobi est egale a : 



-1 
-2 




/ 



J = M~ 1 N 



-"21 
"22 



— "12 
an 





— "13 
an 

-"23 
a22 



"11 

-"2n 



\ 



"n-ln-1 




d'ou, 



3. A 












-1 


















max; E"=i,^i|-^fl 
:m ax* i^ti E"=i,^i|oy'l < 1 
\ 

-1 = D - E - F , d'ou 
/ 



i ^ 
1 1 

UJ 


J + F) 

n 




0^-1 


i , 




Q Q 1-U 


1 - w " 







— w(l — w) 


1-w 


I) 








1-OJ I 



ce qui donne sp(H(co)) = {1 — lo}, soit p(H(uj)) = |1 — w|. Par consequent, la methode 
de relaxation converge pour to G]0, 2[. 
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Reponse 2 
1. 



D~ l E 



D~ X F = 



J = D- 1 E + D- 1 F 








_ "21 





«22 


"31 


_ £32 


»33 


£>33 


_ gnl 


_ a n2 






ai2 




an 









d'ou : 





a 13 
ail 
"23 
"22 



\ 



o J 



. "2ra 
"22 



.SI 

si 



= 3 



Si on note par dij les coefficients de la matrice P a JP a 1 , on obtient : 



— 6j -Pet J Pql &j 

■ — P^tx J Q~ ^3 

= a J ei P a 2^k=i c kj£k 

— r\, ^ 3 rv^ 1 r* ■ ■ — rJ" 3 n . . 



d'ou : 

_ f si i = j 

13 ~\ -c^^t si i + j 

2. Si A est tridiagonale, les coefficients non nuls de A situes sur la diagonale <Hj avec 
i — j = 0, la sous-diagonale atj avec i — j = 1 et la sur-diagonale avec i — j — — 1. 
D' autre part, les coefficients non nuls de la matrice — _E sont les coefficients de la sous- 
diagonale et les coefficients non nuls de — F sont les coefficients de la sur-diagonale. 
Ce qui donne que : 

PcJPa 1 =aL + a _1 U = J{a) 

Done, J(a) et J sont semblables et par consequent, elles ont le meme spectre. 

3. D'apres un theoreme du cours, les valeurs propres {A G sp(H)} et les valeurs propres 
{/i G sp(J)} sont liees par la relation suivante : 

A 2 = A/! 2 

car u> = 1. D'ou, chaque valeur propre /i G sp(J) est associee a deux valeurs propres 
de H : A = et A = fi 2 . D'autre part, det(fl') = det(J - Ly 1 det([7) = 0, done A = 
est une valeur propre de H, ce qui donne : 

sp(H) = {0}U{/i 2 , M G sp(J)} 
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Reponse 3 1. 



1 



— 2ui +U2 
Ml — 2u-2 +U3 

u k -i —2u k +u k +i = fk 



. un—i — 2un 
ce qui donne le systeme AU = F avec : 



/ 2 -1 
-1 2 -1 
0-12-1 



A = 



h 2 



\ ••• 



2. A = D - E - F avec D = £l , 



/ 
10 
10 



\ 



V 



10/ 





-1 



= h 
= h 




= fk 




= fN 


( h \ 

h 
h 


et F = 





et F = — 



V 



/ 1 
1 
1 



\ 



/ 



d'ou : 



J = D^E + D^F = ^ 



/ 1 
1 1 
10 1 



V o 



\ 



10/ 



J est done tridiagonale. D'apres un exercice du chapitre precedent : 

sp(J) = {cos(-^j-);fc = 1,. . . ,n} 
cc qui donne p(J) = cos(^q-j-). D'apres l'exercice precedent : 

sp(H) = {0} U {cos 2 (^) ; k = 1, . . . ,n} 
Ce qui donne : p(-ff) = cos 2 ( ) . 

3. Les deux methodes sont convergentes et la convergence de la methode de Gauss- 
Seidel est plus rapide que la convergence de la methode de Jacobi car p(H) = p(J) 2 < 
p(J) < 1. 

4. D'apres un theoreme du cours : 



Reponse 4 



44 



Partie I : 

1. \\K ||2 = v p(K*K) , comme K est symetrique, d'apres un theoreme du cours : 

3Uunitaire/ K — U* diag(\i)U 
avec sp(K) = {A^}. D'ou, K* K = U*diag(\Xi\ 2 )U, ce qui prouve que : 
p{K*K) = max |A t | 2 = (max |A,|) 2 = p(K) 2 

i i 

2. p(KM) < \\KM\\ 2 < ||ir|| 2 ||M||2 = p{K)p{M), car K et M sont des matrices 
symetriques. 

3. a) A G sp(K) alors, il existe u / vecteur propre de K associe a A tel que : 

Ku = Ait 

Comme I + K est inversible, alors A / -1 (sinon (J + .K^u = et u 7^ 0). 
D' autre part : 

& = \u =^{1 + K)u = (A + l)u 

Ku = (X+ \)K(I + K)- 1 u 
=> Au = (A + 1)JST(7 + if)"^ 

=^i<:(7 + ir)- 1 «= T A xU 

b) (3 G sp(K(I + K)^ 1 ) alors, il existe t)^0 vecteur propre de -ff(7 + K)^ 1 associe a 
P tel que : 

K(I + Ky 1 v = /3v 
Comme u^O, alors u = (J + K)~ 1 v / et on a v — (I + K)u, d'ou : 

K(I + K)- 1 v = f3v =^ Ku = /3{I + K)u 
=^ (1 - P)Ku = [3u 



ce qui montre que (3 =fc 1 et on a 



ce qui donne : 

l^P e SP{K) 

4. I J + K est symetrique et elle est definie positive si et seulement si : 

sp(^I + K) C]0,+oo[ 

Supposons que p{K(I + K)^ 1 ) < 1. Soit 7 G sp(\l + K), alors , il existe w / tel 
que : 

(-J + K)w = -yw 

ce qui donne Klo = (7 — done, A = 7 — | G sp(K). D'apres a), A = 7 — \ =fc — 1 
et on a : 

7-| _7" 1 



Or piKil + K)- 1 ) < 1, d'ou : 



7+5 7+5 7+5 
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ce qui donne : 

-1 < — i-p - 1< 1 

7+1 

d'ou : | < 7 + |, ou encore : 7 > 0. Supposons que \l + K est definie positive. Soit 
f3 G sp(K{I + K)- 1 ). D'apres b) : 

I 3 + let T^p G sp{K) 

d'ou I + G sp{\l + K). Par consequent : 

2 + 1-/3 2(1-/3) 

ou encore : 

ce qui donne {1 + /3 > et 1 - j3 > 0} ou {1 + f3 < et 1 - (5 < 0}. La premiere 
possibilite donne — 1 < (3 < 1 et la deuxieme est impossible, d'ou, ce qu'il faut. 
Partie II. 

1. De la premiere equation, on a : 



x 



(fc+i) 



2> = —{D + H)-Wx (k) + (D + H)-^ 



La deuxieme equation devient alors : 

x (k+i) = ( D + v yi H ( D + H)~ 1 Vx i - k) — (D + Vy^iD + H)'^ + b 

ce qui donne une methode iterative de matrice (D + V)~ 1 H(D + H)~ 1 V , done, la 
convergence de la suite (x k ) est equivalente a p((D + Vy 1 H(D + H)~ 1 V) < 1. 

2. a) On a : 

(D + Vy 1 H(D + H)~ 1 V =(I + D- 1 Vy 1 D~ 1 H(I + D- 1 Hy 1 D~ 1 V 

= (I + cy 1 B(i + b)- 1 c 

= (7 + C^Bil + B)~ 1 C(I + cy\i + C) 

d'ou, la matrice (D + Vy 1 H(D + HyW est semblable a B(I + By 1 C(I + C) _1 , 

done, elles ont le meme spectre et par consequent le meme rayon spectral. 

b) 

b(i + By 1 =(B + I- /)(/ + By 1 = (B + /)(/ + By 1 - 1(1 + By 1 
= I-(I + B)- 1 

et il est clair que I— (I+By 1 est une matrice symetrique. Meme chose pour C(/+C) _1 . 
c) On a : 

p((D + Vy 1 H(D + H)~W) = p(B(I + By L C{I + C)" 1 ) 

Comme B(I + B)^ 1 et C(I + Cy 1 sont deux matrices symetriques, d'apres la partie 
I : 

P (B(i + By 1 c{i + cy 1 ) < p(B(i + s)- 1 )p(c(i + cy 1 ) 

3. D'apres la partie I. question 4. p(B(I + B)~ r ) < 1 si et seulement si |J + B = 
D _1 (|D + H) est definie positive ce qui est encore equivalent a \D + H est une 
matrice definie positive car D = cl avec c > 0. On montre de meme que : \D + V 
est definie positive si et seulement si p{C(I + C) _1 ) < 1. Par consequent, |D + H et 
\D + V sont definies positives nous donne que p(B(I + By 1 )p(C(I + C)^ 1 ) < 1 , ce 
qui donne, d'apres 2. a) : 

p{{d + v)- 1 h{d + h)- 1 v) < 1 

D'ou, d'apres 1., la convergence de la methode iterative (1). 
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Chapitre 3 

Methodes iteratives pour le 
calcul des valeurs propres et 
des vecteurs propres 



3.1 Introduction 

II est connu, d'apres lc theoreme d'Abel, qu'il n'existe pas de methode qui, 
au bout d'un nombre fini d'operations, donne les racines d'un polynome de degre 
n > 5. D'autre part, il est facile de verifier que tout polynome de degre n > 2 
ecrit sous la forme : 

p(x) = x n + a 1 x n ~ 1 H ha„ 

est egal, a un coefficient multiplicatif pres, au polynome caracteristique de la 
matrice 

/ -ai -02 -o„ \ 

1 



V 



1 



/ 



Par consequent, les methodes de calcul des valeurs propres d'une matrice A 
sont, en general, des methodes iteratives d'approximation, qui convergent (dans 
un sens a preciser) vers les valeurs propres. 

Unc facon pour approcher les valeurs propres d'une matrice consiste a construire 
une suite de matrices scmblablcs a A : 



A k 



P h ~ X AP k 



qui convergent vers une matrice dont les valeurs propres sont faciles a calculer, 
diagonalc ou triangulaire par excmplc. 



48 



3.2 Methode de Jacobi 



La methode de Jacobi s'applique uniquement au calcul des valeurs propres 
et des vecteurs propres des matrices symetriques. 

On a vu, chapitre 1 corollaire 1.1, qu'une matrice symetrique est diagonali- 
sable dans K et possede une base orthonormee de vecteurs propres, c'est a dire : 
il existe une matrice orthogonale O et n nombres reels Ai > A 2 > • • • > A„ tels 
que 

O t AO = diag(\i, ■ ■■ , A„). 

De plus, la ieme colonne de O est egale a un vecteur propre associe a la valeur 
propre A^. 

DEFINITION 3.2.1 Soient p et q deux entiers verifiant 1 < p < q < n, et 9 
un nombre reel. On definit la matrice de rotation dans le plan (p, q) et d 'angle 
9 par : 

( 1 \ 
1 



n = 



cos9 



sin9 



—sin9 



cos9 



V 1 7 

Dans toute la suite et a chaque fois qu'on utilise une rotation d'angle 9 dans le 
plan (p, q), l'entier p est suppose strictement inferieur a l'enticr q. 



3.2.1 Resultats preliminaires 

On remarquc facilcment les points suivants : 

- La matrice de rotation f2 est une matrice ortogonalc : 

n T n = i 

- Dans la transformation A = (a^) — ► B = (bij) = QACl , avec fl la 
rotation d'angle 9 dans le plan (p, q), seules les pieme et gieme colonnes et 
ligncs de la matrice A sont modifiees : 

= aij si i^p,q et j ^ p, q, 
b P i = a P iCos9 — a q isin9 si i ^ p, g, 
b ql = a pi sin9 + a qi cos9 si i ^ p, q, 
b pp = a pp cos 2 9 + a qq sin 2 9 — a pq sin29 ^ ' ' 

b qq — a pp sin 2 9 + a qq cos 2 9 + a pq sin29 
b pq = b qp = a pq cos29 + app ~ aqq sin29 

- La normc de Frobenius est invariante par transformation orthogonale, 
d'ou, 

n n 

E h l = E 4 (3-2) 

i,j = l i,i = l 
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PROPOSITION 3.2.1 Soient A = (ay) une matrice symetrique et B = 
(bij) = £IA£1 T , oil Q, est la rotation d'angle 9 dans le plan (p,q), alors : 

(l) bpp -\- bqq Opp ~t~ Gjqq 

(ii) bpp — bqq = (a pp — a qq ) cos 29 — 2a pq sin 29 

(iii) b 2 pp + 2b 2 pq + b\ q = a 2 pp + 2a 2 pq + a 2 qq 

(iv) E & « = E4- 2 «-^) 

DEMONSTRATION : 



Los formules (i) et (ii) se deduiscnt facilcmcnt dcs formulas (3.1l. On verific 
facilement par le calcul que : 

cos 9 sin 9 
— sin 9 cos 9 

d'ou, si on applique la norme de Frobenius a cette egalite et en utilisant encore 
une fois que la norme de Frobenius est invariante par transformation ortogonale, 



[ bpp 


bpq 


H 


/ cos (9 


— sin ( 




bqq , 




, sin 9 


cos 9 



on obtient (iii). D'autre part, la formule (3.2 1 nous donne : 



£$ = (£$) + (£ 



= £« 



2 

ij 



= (£4) + (£ 

et la premiere formule de (13 . 1 h nous donne : 



2\ , 2 , 2 
i" a pp T 



^=£ 

»^P,9 «^P:9 
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il suffit done d'utiliser (iii) pour obtenir (iv). 

3.2.2 Principe de la methode de Jacobi 

La methode de Jacobi consiste a construire une suite de matrices scmblablcs 
a la matrice A, de la maniere suivante : 

Ai = A 

Ak+i = il k A k Q^ 

ou fl k est une matrice de rotation dans un plan (p, q) et d'angle 9 bien choisis. 
II est clair que (A k ) est une suite de matrices semblables a A et on a : 

A k+1 = {Q k n k -i ■ ■ • ni)i4(n*n k _i ■ ■ • «i) T 

On veut que A k tende vers une matrice diagonale (qui est necessairement sem- 
blable a A) quand k — > oo. Par consequent, si on note par A k+ \ — (bij(6,p,q)) 
et A k — (ciij ) , on a interet a choisir 9 et (p, q) qui minimisent la somme des 
carres des elements hors diagonaux de la matrice A k+ i, e'est a dire : 



£ bij{6,p, q) = min ^ b 2 } (9,p, q) < > </; 



£4 
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D'apres la formulc (iv) de la proposition 3.2.1 il faut ct il suffit dc choisir p, q 
et 9 tels que : 

done, p, q sont choisis tels que : 

| a m |= max | a pq \ 

v+q 



et 9 tel que : 



bpq(9,p,q) = 



PROPOSITION 3.2.2 Pour tout (p,q) tel que a pq ^ 0, il existe une unique 
valeur non nulle de 9, — f < 9 < ~, telle que : 

b pq (9,p,q) = 



DEMONSTRATION : D'apres les formules (|3_lj), on a : 

<>,,,, I) < • cotg(29) 



or, la fonction cotg(29) est bijective de ] — J, 0[U]0, ^] dans K. D'ou, l'existence 
et l'unicite de 9 e] - f ,0[U]0, |] tel que b pq (9,p,q) = 0. 

3.2.3 Convergence de la methode de Jacobi 

Soit Ai — A = (ay) = (oy) et = (a^) la matrice de l'iteration k de la 
methode de Jacobi. Soit pk, qk tels que 

Soit 6» fc G] - f ,0[U]0, f] l'unique solution de 

cotg(29) = (3.3) 

Za Pkq k 

et f2fc la matrice de rotation dans le plan (pj,, q k ) d'angle 9^. On definit la matrice 
Ak+\ de l'iteration k + 1 par : 



THEOREME 3.2.1 La swiie (A^) de matrices obtenues par la methode de 
Jacobi est convergente, et 

lim A k = D 

k >oo 

oil D est une matrice diagonale qui a le meme spectre que A. 
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DEMONSTRATION : On note par 

A k = D k + E kl avec D k = diagia^} 



£ k -|| E k \\p= ^2 



k |2 



II est clair que : 



|< 9fc \ 2 <s k <n(n-l)\a k pkqk | 5 



(3.4) 



D'autre part, d'apres la formule (iv) de la proposition 3.2.1 

|2 



e k+1 = £ fe - 2 | a pkqk 



car a p + q \ - b Pkqk 



ce qui donne : 
avec 

Par consequent. 
Ce qui prouve que 



0, et de (3.4) on tire que : 



a k | 2 < \ 

PkU I - n ( n _ 1) 



< £fc+i < re fe < r 2 e fc _i < • • < r fc £i 



0<r=(l-^ -)<1 



lim £fe = 

k — >oo 



Urn E k = 

— >oo 



Soit i 6 {l,2,..,n} fixe. D'apres lcs formulcs (3.1l, (3.3 > , pour tout k 6 N, si 
i ^ p k ct i ^ q k on a : 

a'L 1 " — a" 



si i = pi on a : 



,fc+i 



= -tg(0 k )a* h9h 
et enfin si i = on a : 



= < Pfc cos2 (^) + a q kqk sin2 ( k) - a k pkqk sin(26 k ) - a k pk 
= ( a L qh ~ a p kPk ) sin2 ( 9 k) ~ a* sin{29 k ) 
= 2a p k q k cot 9( 2d k)sin 2 (e k ) - a k pkqk sin{26 k ) 
[2cotg(29 k )sin 2 (6 k ) - sin(29 k )]a k 



k+i _ k 
q k q k q k q k 



ak p k p k sin2 i°k) + a g k q k cos2 ( e k) + a k kqk sin(28 k ) - a 

( a Q k q k ~ a p kPk ) COs2 ( 9 k) + a p k q k Sin ( 20 ^ 

2a k pkqk cotg(29 k )cos 2 (9 k ) + a k kqk sin{29 k ) 
[2cotg(29 k )cos 2 (9 k ) + sin(29 k )]a k kqk 



k 

p k p k 



or | 9 k \< j ce qui donne \tg(9 k )\ < 1, d'ou 



\at + } - at r , | < \a k „ I et \a k+ } — at „ I 
i p k p k p k p k > — i p k q k > i q k q k q k q k > 



< K k q k \ 
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D' autre part 

VfcGN, \a k pkqk |<V^<(^) fe_1 V^ 

par consequent 

VfceN, \a k i t-a k ii \<{^r) k - 1 ^ 
d'ou, la serie numerique J2k I a ii +1 — a ii I est convergente, ce qui prouve que la 
serie X)fc(°ii +1 ~ a u) esij convergente ou encore, la suite des sommcs partielles 



Sfc = E(°« 1 -°«) = a « +1 -«« 
1=1 

est convergente. Par consequent, pour tout i fixe, la suite (a k A converge, ce qui 
prouve que la suite de matrices (D k ) est convergente vers une matrice diagonale 
D. D'autre part, on a : 

det(A k - XI) = det(A - XI) 
et par continuity, en utilisant le fait que E k — > et D k — > D, on a : 
lim det(A fc - A/) = cfet(D - A/) 

/c^oo 

Par consequent, A ct D ont le meme polynome caracteristique et la suite de 
Jacobi (Ak) converge vers D. 

THEOREME 3.2.2 On suppose que toutes les valeurs propres de A sont dis- 
tinctes. Alors la suite de matrices orthogonales 

Uk = ^k^k-l ' ' ' "1 

converge vers une matrice orthogonale U dont les colonnes de la matrice U T 
sont des vecteurs propres de la matrice A. 

DEMONSTRATION : Soit 

5 = \ min | A j — Xj \ > 

3 i^j 

avec 

D = diag{Xi, A 2 , • • • , A„} = lim D k 

k^oo 

et Dk = diag(a k i ). Soit ko > 1 tel que pour tout 1 < i < n, on a : 

V/c>fc , \a k u -Xi\<5 

d'ou 

Vfc>fc , Vi^j, \a k i -a k j \>S 

en particulier 

Vk>k 0l \a k kqk -a k kPk \>5 

ou encore 



Vfc>fc , 2|o* || cot5r(2^ fe ) |> 5 
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cc qui prouvc que 

Vk>k , | tg(20 k ) |< 2^|^i 

D' autre part, 

I a \<\ tg(a) |, Va e] - tt/2, tt/2[ 

d'ou 





D'autre part, en utilisant le fait que la norme de Frobenius est invariante par 
transformation ortogonale et que U k+ i — fl k +iU k ■ 

|| t/fe+i - U k HI = || (fi fc+1 - 7)C/ fe HI 

= ll(fi*+i-/) III 

= 2(cos# fe+1 -l) 2 + 2(sin(^ +1 )) 2 

= 8sin 2 (0 fc+1 /2) 

< 2 1 e k+1 1 2 

Par consequent 

Vfc>fco, \\U k+1 -U k \\ F <V2^p<^(Vr) k 

d'ou, la serie ^ || — t^fc ||f est convergente, done, la serie J2(U k +i — U k ) 

est convergente. Ce qui prouve que la suite des sommes partielles 



S k = J2 Ui+i -Ui = u k+1 - u x 



1=1 

est convergente. Done, la suite de matrices orthogonales (Uk) est convergente 
vers une matrice orthogonale U et par passage a la limite on a : 

D = UAU T 

d'ou, la ieme colonne de U T est un vecteur propre de A associe a Aj. 

3.2.4 Mise en oeuvre de la methode de Jacobi 

Dans la pratique on ne cherche pas a calculer le nombre 9 k pour determiner 
les coefficients des lignes et des colonnes p k , q k de A k+ i ; on pcut les obtenir 
par les relations algebriques suivantes : soit 

_ <3k<?fc "PfcPfc 



2 <« 
et 

t = tg9 k 

on a les relations trigonometriques suivantes : 

c = cos6 k = - 
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et pour que aj" 1 . 1 soit nul, il faut que 



"PkQk 

avec la condition t e] — 1, 0[U]0, 1]. D'ou, t est egal a l'unique solution du trinome 

t 2 + 2i?i-l=0, t e] - 1,0[U]0,1]. 

On obtient alors les formules suivantes : 

a pi =ca ni-* a ni 31 i^Pk,qk, 
= sa Li + CO* i si i^p k ,q k , 

PfcPfc PfcPfc PkQk 

r .k+l n k i 

„/c+l _ k+1 _ n 

■ a p k qk — a q k Pk ~ u 

REMARQUE 3.2.1 Un element annule a I 'iteration k pent devenir non nul 
a une iteration I > k. 

On distingue trois strategies pour le choix du couple (p, q) : 

1. Methode de Jacobi classique : on choisit I'un des couples pour lesquels 

\ap q \= max I a-ij I • 

2. Methode de Jacobi cyclique : on annule successivement tous les elements 
hors-diagonaux par un balayage cyclique, toujours le meme : par exemple, 
on choisit les couples (p,q) dans I'ordre suivant 

(1,2), (1,3),--- ,(l,n);(2,3),--- , (2, n); • • • ; (n - l,n) 

3. Methode de Jacobi avec seuil : on procede comme dans le cas de la methode 
de Jacobi cyclique mais on omet d'annuler les elements hors-diagonaux 
dont le module est inferieur a un certain seuil. 

EXEMPLE 

Soit 

2 -1 1 
1=1-1 3 -4 

1 -4 3 

On pose Ai = (ajj) = A. La matrice A 2 — (of-) de la premiere iteration de 
la methode de Jacobi classique est obtcnuc comme suit : on cherche d'abord 
la position (p, q) du plus grand coefficient en valcur absolue en dehors de la 
diagonale de A\ : c'est la position (2,3), on calculc le nombre : 

_ a 33 ~ a 22 _ 3 — 3 _ 
202 3 ~& 
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on resout l'equation : t 2 + 2Rt— 1 = t 2 — 1 = 0, on prcnd l'unique solution qui est 
dans ] — 1, 0[U]0, 1] ( dans notre cas c'est t = 1 ), on pose c = l/Vl + t 2 = 1/V%, 
s = t/VTTt 2 = l/s/2 et enfin on a : 

o?i = «n = 2 

a 21 — a 12 = c X a 21 ~ s X a 31 — -1/V^ - l/\/2 = -2/y/2 
a 31 = a 13 = s X a 21 + c X a 31 = -l/\/2 + l/\/2 = 
a 32 = a 23 = 

a 22 = a 22 - * X a 23 = 3 + 4 = 7 
a 33 = a 33 + * X a 23 = 3 - 4 = -1 

ce qui donne : 



2 


-2 
V2 





-2 


7 











-1 



3.3 Transformation de Householder et la facto- 
risation QR 

3.3.1 Transformation de Householder 

A tout vecteur ti£l", on associe la matrice 

I - 2uu T T uu T . 1 1 1 - - o . , „ 

/, si u = 

la matrice H est appelee la matrice de Householder et le vecteur u est appele le 
vecteur de Householder. (On note parfois H(u) pour montrer la dependance du 
vecteur u) . II est facile de voir que les matrices de Householder sont symetriques 
(H = H T ), orthogonales (H T H = HH T — I) et dependent seulcmcnt de la 
direction du vecteur de Householder. 

Les matrices de Householder ont deux proprietes importantes : 
- pour tout vecteur a ^ et tout vecteur ^ 6 ^ a, avee || a \\2—\\ b H2, il 
existe un vecteur de Householder u tel que : 

Ha=(I-^f)a = b 

soit encore 

{-f^)u = b-a (3.5) 

d'ou u est un multiple de b — a ^ 0. D'autre part, si u = b — a, en utilisant 
le fait que b T b — a T a, on obticnt : 

(o — a) T a 

Ha =a-2 1 / db-a) 

(b-a) T (b-a) t (3 6) 



= o + {-l-2 {h 0) ° A (b-a) = b 
(b — a) 1 (b - a) 
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- tout vecteur transforme par H a une forme speciale. Si on applique H a 
un vecteur c, on obticnt : 



He 



it uu \ 



T 

,U C. 



(-)" 



(3.7) 



d'ou, lc vecteur He est la difference cntre lc vecteur c et un multiple du 
vecteur de Householder u. On deduit aussi de (3.7) que l'application de H 



a c ne change pas les composantes qui correspondent a des composantes 
nulles de u et que c est invariant par H si u T c — 0. Enfin, le calcul de He 
necessite seulemcnt lc vecteur u et le scalaire (3. 



Par exemple, soit u = (—1, —2) , on a 



= 5 



H = I 



uu 



1 
1 



2 4 
5 5 



3.3.2 Factorisation QR 

Soit A — (a.ij) une matrice carree reelle quelconque. La factorisation QR 
consiste a trouver une matrice Q orthogonale et une matrice R triangulaire 
superieure telles que : A = QR. En utilisant les proprictes des matrices de 
Householder, on va construire une suite de (n — 1) matrices de Householder 
telles que : 

H n -iH n _ 2 ■ ■ ■ H 2 HiA = R 

oil R est une matrice triangulaire superieure. 

La premiere etape consiste a construire une matrice de Householder Hi qui 
transforme a\ (la premiere colonne de A) en un multiple de e\, ce qui revient a 
annuler les composantes 021,031, • • ■ , a n \. La norme Euclidienne est invariante 
par transformation orthogonale, par consequent, on cherche U\ tel que 



Hid! - (/ 



T 

U\U{ 



)ai = ± || a\ || 2 ei 



( ri1 \ 




V / 



avec rn 



a\ H2. Nous savons que u\ doit etre un multiple de ± || a\ \\ 2 e%— a\ 
(voir (3.5)). Comme H\ depend uniquement de la direction de u\, on prendra 
ui comme suit : 

/ an - rn \ 
021 



"1 



Par definition de u\, le signe de rn peut-etre positif ou negatif. Pour eviter les 
erreurs numeriques, le signe de rn est en general choisi de la maniere suivante : 



sign(rn) = -sign(cin)- 



57 



pour que (3 = \ || u \\ 2 soit le plus grand possible, car on divise par (3. Apres 
application de la matrice dc Householder Hi, la premiere colonne de A 2 = Hi A 
est un multiple de ei : 



Hi A 



( rn 


V o 



(2) 

Hi 

l 22 



l n2 



■ a {2) \ 

"in 



1 (2) 
l 2n 



a {2) I 



Tous les elements de A sont modifies. Pour % allant de 2 a n, on a : la colonne 
di est remplacee par a,- L 



u 1 dj 



Ui. 



A la deuxieme etape, il faut que la matrice dc Householder H 2 ne change pas 
la premiere colonne ct la premiere ligne de A 2 . Pour ccla, il suffit dc choisir le 
vecteur dc Householder u 2 orthogonal a ei, e'est a dire la premiere composante 
de u 2 nulle. Avec un tel choix, l'application de H 2 a un vecteur quclconque ne 
change pas sa premiere composante ct l'application de H 2 a un multiple de e\ 
( commc le cas de la premiere colonne dc A 2 ) lc laissc inchangc, d'ou : 



/ (2) 



a = 



12 

l 22 
2 (2) 

h2 



\ 



b = 



I ^ \ 

r 22 




/ 



, u 2 = b — a = 



U 2 2 v v / 





(2) 

1 2 2 - r 22 

a (2) 

"32 



V ^ J 



/ (2) (2) 

avec r 22 = ±\J (a 22 >) 2 H h (a y n2 > ) 2 pour que ||a|| 2 = ||6|| 2 . 

Au bout de (n — 1) etapes de reduction de Householder, nous obtenons : 

H n -i ■ ■ ■ HiA = R 

oil R est une matrice triangulairc superieure. 
On pose 

Q T = H n _i ■ ■ ■ H 2 Hi 

soit 

Q = H\H 2 ■ ■ ■ H n _i 
ce qui donne la factorisation QR de A : 

Q T A = Ron A= QR, 

avec Q est une matrice orthogonalc ct R est une matrice triangulairc superieure. 



REMARQUE 3.3.1 Dans le cas d'une matrice A carree a coefficients com- 
plexes, on pent refaire les memes calculs, en remplagant u T par u* dans la 
definition de la matrice de Householder, pour obtenir une factorisation A = QR 
avec Q une matrice unitaire et R une matrice triangulaire et dans ce cas les 
matrices de Householder sont hermitiennes et unitaires. 
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EXEMPLE 

Soit 



A = 




le calcul suivant donne la factorisation QR de la matrice A. La norme 2 dc la 
premiere colonne de A est egale a \f \ + 2 2 + 2 2 = 3, d'ou, le vecteur Ui de la 
matrice Hi est egal a : 

/ 1 + 3 

ui = 2 

V 2 

ce qui donne : /?i = i ||| = 12 et 



T / -1/3 -2/3 -2/3 \ / -3 -1/3 -8/3 

Hi = h-^=\ -2/3 2/3 -1/3 , HiA=\ -8/3 -7/3 
» -2/3 -1/3 2/3 / V -5/3 -7/3 



Pour calculer le vecteur w 2 de la matrice H 2 , on doit calculer : r 22 = ± ^/64/9 + 25/9 = 
±V§9/3, d'ou : 





U 2 



cc qui donne : f3 2 = |||w2||| = 18.27 et 

T / -0.33 -0.66 -0.66 \ / -3 -0.33 -2.66 
H 2 Hi = {I z -'^^)Hi = 0.91 -0.38 -0.07 , R = H 2 H X A = 3.14 3.21 

^ 2 \ -0.21 -0.63 0.74 / \ -0.74 

et la matrice 

-0.33 0.91 -0.21 
Q = HiH 2 = | -0.66 -0.38 -0.63 
-0.66 -0.07 0.74 



3.4 Methode QR 

Soit Ai = A unc matrice quelconque ; on ecrit sa factorisation QR, soit 
Ai = QiRi, puis on forme A 2 — RiQi ; on ecrit sa factorisation QR, soit 
A 2 = Q 2 R 2 , puis on forme la matrice A3 = R 2 Q 2 - A l'etape k, on ecrit la 
factorisation QR de A^, soit A^ = QkRk puis on forme A^ + i = RkQk et ainsi 
de suite. On obtient ainsi une suite de matrices Ak qui sont toutes semblables 
a la matrice A, puisque : 

A 2 = R1Q1 - Q\AQi 

A k+ i = R k Q k = Q* k A k Q k 

= ■■■ = (QiQ2---Q k )*A(QiQ 2 ---Q k ). 
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Dans la pratique, avant d'appliqucr la methode QR, on commence par mettre 
la matrice A sous la forme d'une matrice de Hessenberg superieure : 



X 


X 


X 


X 
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X 


X 
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X 














X 


X 


X 


X 
















X 


X 


X 


















X 


X 


J 



V 

semblable a A, en utilisant les matrices de Householder (voir exercice 3). 

L'interet de la mise sous la forme de Hessenberg de la matrice A est que la 
suite des matrices Ak de la methode QR d'une matrice de Hessenberg restent 
sous la forme de Hessenberg (voir exercice 4), ce qui rcduit considerablement le 
temps de calcul. Vu la forme des matrices de Hessenberg, il suffit d'annuler les 
coefficients de la sous-diagonale. Si a l'etape k un element de la sous-diagonale 
de la matrice Ak est numeriquement nul, on partage la matrice en deux sous- 
matrices 
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X 


X 


X 


X 


X 


A\ 
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X 








X 
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X 
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X 
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X 








X 
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X 
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X 


X 


X 








X 




A\ 




X 










X 






X 












X 




X 














X 


X 



et il est clair que 



spectre(Ak) — spectre(Al) U spectre{A\) 



ce qui permet de rcduirc l'ordre de la matrice. On garde le bloc A\ en memoire 
et recommence la methode QR avec le bloc A\, ainsi de suite, jusqu'a ce que 
le bloc infcricur soit d'ordrc deux (resp. un), dans ce cas on a determine deux 
valeurs propres (resp. une valcur propre) . Puis on redemarre avec le bloc qui est 
juste avant. 



3.5 Exercices : chapitre 3 

Exercice 1 Soit 

/ 2 -1 3 1 \ 

-13 2 

3 2 1-1 
\ 1 0-12/ 

Determiner la matrice semblable a A donnee par la premiere iteration de la methode 
de Jacobi. 



GO 
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-1 


2 


2 






-2 


3 


4 


1 






1 


2 





-1 
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1 





-1 


2 





Exercice 2 Soit 



1. Determiner la matrice Ai semblable a A et qui a la forme d'une matrice de 
Hessenberg. 

2. Determiner la matrice A2 semblable a A\ donnee par la premiere iteration de la 
methode QR. 

Exercice 3 (Forme de Hessenberg) 

Soient A — (aij) = Ai = (a}j) une matrice carree d'ordre n et (ej)j = i, n la base 
canonique de E n . 

1. Montrer qu'il existe un vecteur iti / de M n orthogonal a e\ tel que : 
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\ 
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X 


X 
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X 











X 


X 


/ 



avec Hi — H(ui) = la matrice de Householder associee au vecteur ui. 
2. On pose A2 = H1A1H1 — H\A\Hi. Montrer que A2 a la forme : 

,2 2 „2 



/ z 2 2 \ 

I an a 12 a ln \ 

22 2 

a 21 a 22 a 2 



A 2 = 



a 3 2 



2 

a 3r7 



V a^ 2 "L / 

3. On suppose qu'a la suite de l'etape k — 1, on a une matrice Ak de la forme 

aL \ 



I 811 012 

a 21 a 22 

a%2 




A k = 



a 2n 



a '.in 



V 







*kk-l 






a kk 
a k+lk 



a nk 



a k / 



Montrer qu'il existe un vecteur Uk =fc de E n orthogonal a e», i = 1, fc , tel que 



Gl 



H k A k a la forme suivante 



&21 





bi2 
622 
632 



6l„ ^ 

&2n 
&3n 











bk+ik bk+ik+i 
6fc+2fc+i 



V 







avec Hk = H(v,k) = la matrice de Householder associee au vecteur Uk- 



4. On pose A k+ i = H k A k H k = H k A k Hk 



( 



A k+ i = 



ni 





7 fc+l 

x 12 

-fc+1 

- t 22 
^32 





Montrer que A^+i a la forme suivante : 



,fc+i 

,fe+l 



,fe+l 



V 







.fe+1 
l k+lk 






,fc+l 
l k+lk+l 



,fc+l 

' l fe+2fc+l 



,k+l 



,fe+l 



5. Deduire que toute matrice carree ^4 est semblable a une matrice qui a la forme 
d'une matrice de Hessenberg : 
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... X X 
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X 


HS = 
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X 



V 







c'est a dire pour tout j — 1, n— 2, le vecteur HSej €< ei, €2, • • ■ , e.3+1 > l'espace 
engendre par ei, 62, ■ ■ ■ , e,j+\- 

6. Montrer que si la matrice A est symetrique alors la matrice HS semblable a 
A obtenue par la methode ci-dessus est symetrique. Deduire que toute matrice 
symetrique est semblable a une matrice tridiagonale. 

Exercice 4 

Soit A une matrice carree d'ordre n qui a la forme d'une matrice de Hessenberg 
(voir exercice precedent). Soient Hi, H2, • • • , H n -\ les matrices de Householder de la 
factorisation QR de la matrice A, c'est a dire : Q — H1H2 ■ ■ ■ -ff„-i=matrice unitaire 
et R = Q* J 4=matrice triangulaire. 

1. Montrer que pour tout k = 1, 2, • • ■ , n — 1, on a : 

{H k e k €< e/fe, e k +i > 
H k e k +i G< e k ,e k +i > 
Mi =fi k,i ^ k + 1, H k ei = 



G2 



2. Montrer que pour tout k = 1, 2, • • • , n — 1 : 

(a) Qe k = H!H 2 ■ ■ ■ H k e k 

(b) Qe k €< ei, e 2 , • • • , e fe+ i > 

3. Deduire que la matrice B = RQ a la forme d'une matrice de Hessenberg. 

4. Conclusion : Montrer que si une matrice A a la forme d'une matrice de Hessen- 
berg, alors, les matrices de la suite (A k ) de la methode QR ont la forme d'une 
matrice de Hessenberg. 



3.6 Corrige des exercices : chapitre 3 

Reponse 1 On trouve que (p, q) = (1, 3), d'ou : R = — 1/6, ce qui donne le trinome : 

t 2 + 2Rt - 1 = t 2 - t/3 - 1 = 

l'unique solution qui appartient a ] — 1, 0[U]0, 1] est egale a t = (1 — V37)/6 w —0.847, 
ce qui donne : 

s = t/\/l + t' 2 w -0.646, c= l/\/l + i 2 w 0.76 

par consequent, si on note par B = la matrice donnee par la premiere iteration 
de la methode de Jacobi classique, on obtient : 



&12 


= &21 


— c x a\2 — s x 


032 ~ 


0.76 x (-1) - (-0.646) x 2 


&14 


= &41 


= cx ai4 — s x 


a 34 ~ 


0.76 x (1) - (-0.646) x (-1) 


&32 


= &23 


= s x 012 + c x 


032 ~ 


(-0.646) x (-1) + (0.76) x 2 


&34 


= 643 


— s x 014 + c x 


a 34 « 


(-0.646) x (1) + (0.76) x (-1) 


622 


= 0,22 


= 3 






644 


= 044 


= 2 






631 


= &13 


= 






6ll 




= on — t x ai3 


w 1 - 


(-0.847) x 3 


633 




= £133 + t X ais 


w 2 - 


(-0.847) x 3 



Reponse 2 

1. On trouve : 



2. On trouve : 



2 


2.45 


1.22 


1.22 \ 


2.45 


-0.33 


-3.96 


-0.7 





2.09 


-2.65 


0.005 








1.39 


2.5 / 


-3.16 


-1.29 


-2.3 


0.2 \ 





2.97 


-0.38 


1.08 








-4.16 


2.31 











2.27 / 


-0.63 


0.55 


0.53 


-0.11 \ 


-0.77 


-0.45 


-0.43 


0.09 





0.7 


-0.69 


0.15 








0.21 


0.97 / 



A 2 =RiQi = 

V 



3 


-2.77 


0.51 


0.09 


2.3 


-1.6 


-0.79 


1.28 





-2.93 


3.38 


1.62 








0.49 


2.22 
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Reponse 3 (Forme de Hessenberg) 

/ all \ / xi \ 



1. On veut que Hi 



all 



%2 





( voir cours ) il sufBt que 



( xi \ 


V J 



V all J \ / 

/ an \ f Xl \ 

xi 



( a " \ 

ahi 

all 



\ aii j 



k et ui = 



ahi 
all 



V oil / 







V J 



Pour que ui soit orthogonal a ei, il faut que xi = a}i et pour avoir 
la meme norme 2, il faut que 

n 

J2(ali) 2 = (aii) 2 + {x 2 f 



d'ou 



X2 — ± 



/ 

<4i ± V^U(«Ji) 



ce qui donne ui 
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V a} nl 



J 

( a " \ 



2. A 2 ex = HiAiHiei = HxAxex = Hi 



a\i 
all 



f xi \ 

X2 





V aii j 



\ / 



2 

a 2 i 




V / 



i?iei = ei, (wi _L ei). 
3. Meme raisonnement que la premiere etape, on cherche 

/ a k kl \ 



u k _L et,i = 1,2, • • • ,k / i?* 



a kk 
k 

a k + lk 



\ <fc J 



( bki \ 

bkk 

bk+lk 



V J 



on trouve 



/ 



u k = 



a k+ik ± yjj2i=k+i( a ik) 2 



a k+2k 



G4 



On a H k d = ei, i = 1, 2, • • • , k , d'ou : 

Vai G< ei, e 2 , • • • , e k >, H k x = x 

Par consequent : 



et 



Vi = 1, 2, • • ■ , fc - 1, H k A k ei = A fc ei 



/ <4 \ 



H k A k e k — H k 



l kk 



l k+lk 



\ <* J 



b kk 




V / 

4. A k+ iei = H k A k H k ei = H k A k ei = A^ei, i = 1, 2, • • • , fc — 1 et 



A k+1 e k = H k A k H k e k = H k A k e k = 



b kk 
bk+ik 




a 



fc+i 

kk 
fc+1 





V / V o y 

5. A la suite de l'etape n — 2, on obtient la forme HS de Hessenberg. 

6. Si A est symetrique, alors A k est symetrique pour tout fc = 1, 2, • • • , n — 2. Done, 
ffS = A n _ 2 est symetrique et il est clair qu'une matrice symetrique qui a la 
forme de Hessenberg est necessairement tridiagonale. 



Reponse 4 



1. Par recurrence sur fc. Pour fc = 1, on a Ui = 



car A\ = A 



\ J 

a la forme d'une matrice de Hessenberg, d'ou, H\e\ G< ei,e2 > car mi G< 
ei, e 2 > et on a bien : Hid = ei pour tout i > 3. Supposons que le resultat est 
vrai pour 1 < fc < I — 1 et montrons qu'il est vrai pour k = I. On a : 

Aiei = H1-1H1-2 ■ ■ ■ H\AH\ ■ ■ ■ Hi- 2 Hi^iei 

G H1-1H1-2 ■ ■ ■ HiAHi ■ ■ ■ H1-2 < ei-i,ei > 



G Hi-^Hi-2 ■ ■ ■ HiA < ei,e 2 , ••• ,ei> 
G Hi-iHi-2 ■ ■ ■ Hi < ei,e 2 , • • • , e l+1 > 
G< ei, C2, • • • , e; + i > 

Done, mj 6< e;,e; + i > et par consequent -ff; verifie : 

Hid = e; - ( "' )3 e ' ) Mi G< e;, e !+ i > 
Hiei+i = e ;+ i - ( "' e p +1> ui G< ej, e i+ i > 
^e, = a - ( -^p 1 u l = e, Vi^l,Vi^l + l 
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2. a) D'apres (1) Hie k = e k pour tout I > fc+1, d'ou : Qe k — H1H2 ■ ■ ■ H n -ie k =HiH2 ■ ■ ■ H k e k 
b) 

Qe k = H 1 H 2 ■ ■ ■ H k e k 

€ H l H 2 ■ ■ ■ -Hfc-i < e k , e k+ i > 
€< ei, e2, • • • , efe+i > 

3. La matrice i? est triangulaire superieure, d'ou : 

Vi = 1, • • • ,n, ii < ei,e 2 , • • • , e; >C< ei,c 2 , • • • , ei > 
ce qui donne : 

Vfc = 1, • • • , n - 2, RQ < ei, • • • , e k >C R < ei, • • • , e fc+ i >C< ei, • • • , ejt+i > 

ce qui prouve que B = RQ est une matrice de Hessenberg. 

4. A chaque etape k de la methode QR, on applique la factorisation QR k A k = 
Q k R k et on pose A k+ i = R k Q k . D'apres les questions precedentes, si A k a la 
forme d'une matrice de Hessenberg alors R k Q k = A k+ i a aussi la forme d'une 
matrice de Hessenberg. Comme A\ = A a la forme d'une matrice de Hessenberg, 
un raisonnement par recurrence nous donne ce qu'il faut. 
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Chapitre 4 



Methodes d'interpolation et 
d'approximation 



4.1 Formule d'interpolation de Lagrange 

Soicnt (xi,fi),i — 0, • • • , n, avec x^ ^ Xj pour i ^ j, (n + 1) -couples reels 
ou complexes. Lc problemc d'interpolation polynomiale consiste a trouver un 
polynome P(x) de degre inferieur ou egal a n tel que : 

P(xi) = fi, * = 0,-.. ,n (4.1) 

THEOREME 4.1.1 Formule d'interpolation de Lagrange 

Soient (xi, fi),i = 0, • • • , n, avec Xi ^ xj pour i ^ j, (n+ l)-couples reels ou 
complexes alors, il existe un unique polynome P{x) de degre inferieur ou egal a 
n verifiant 

DEMONSTRATION: Unicite : Supposons l'existence de deux polynomes 
Pi(x) et P2(x) de degre inferieur ou egal a n et verifiant (4.1). Par consequent 
P\(x) — Pz(x) est un polynome de degre inferieur ou egal a n et possede (n + 1)- 
racines ce qui donne Pi — P2 = 0. 
Existence : On pose 



Li(x) 



avec 



Le polynome 



(x - x )(x - xi) ■ ■ ■ (x - Xj-i)(x - x i+ i) ■■■(x- x n ) 

(xi -x )--- (xi ~ Xi-i)(xi - x l+ i) ■ ■ ■ (xi - x n ) 
v(x) 

(X - Xi)v'(Xi) 

v(x) = {x- x )(x - Xi) ■ ■ ■ (x - x n ) 



P(x) = J2fMx) 



est un polynome de degre inferieur ou egal a n et verifie les conditions (4.1 1. 
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EXEMPLE: On chcrche un polynome de degre inferieur ou egal a 2 tel que : 

P(0) = 1, P(l) = -1, P(4) = 5 
D'apres le theoreme precedent, il y a un scul polynome donne par : 
P(x) = 1 x L (x) + (-1) x Li (a:) + 5 x L 2 (aO 



avcc 



d'ou 



L (a;) : 
Li(ar) = 
L 2 {x) ■■ 

P(x) 



(x 


-l)(x 


-4) 


(0 


-1)(0 


-4) 


(x 


-0)(x 


-4) 


(1 


-OKI 


-4) 


(X 


-0)(x 


-1) 


(4 


-0)(4 


-1) 


x 2 


-3x + 


1 



En general les couples (xi, fi), i = 0, • • • , n sont les resultats d'une experience. 
Par consequent les nombres fi, i = 0, • • • ,n sont les valeurs d'une fonction / 
aux points Xj, i = 0, • • • , n et le polynome donne par la formule d'interpolation 
de Lagrange coincide avec cette fonction aux points Xj, i = 0, • • • , n. Done la 
formule d'interpolation de Lagrange est un moyen pour approcher unc fonc- 
tion / lorsquc on connait les valeurs prises par / aux points Xj, i = 0, • • • , n. 
La question qui se pose maintenant est de savoir estimer l'erreur commise en 
approchant /(x) par P(x). 

THEOREME 4.1.2 Si f est (n + V)-fois derivable, alors, pour tout x el, il 
existe c e I ( I est le plus petit intervalle contenant Xo, x\, ■ ■ ■ ,x n ,x) tel que : 

f(x)-P(x)=v(x)L^ T y- (4.2) 
avec v(x) = (x — x )(x — Xi) • • • (x — x„). 

DEMONSTRATION: Si x est egal a l'un des x i; il n'y a rien a demontrer. 
Supposons que x ^ Xj, i = 0, • • • , n. Soit K e R tel que : 

F(x) = f(x) - P(x) - Jf«(x) 

vcrifie F(x) = 0. Done, la fonction F(x) possede (n+ 2)— zeros dans l'intervalle 
/. En appliquant le theoreme de Rollc a F(x) puis a P'(x) puis a P"(x) • • • , on 
obticnt que F(" +1 )(x) possede au moins un zero dans l'intervalle /, d'oii : 

3ce //P ( " +1) (c) = 

or = et w( ,l+1 )(x) = (n + 1)!, cc qui donne : 

/ (n+1) (c)- K(n+l)\ = 

d'ou 

(n + 1)! 

En remplacant la valeur de X dans l'expression de P(x) et en utilisant le fait 
que P(x) = 0, on obtient le resultat donne par le theoreme. 
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4.2 Algorithme de Neville 

L'algorithme de Neville est une methode tres pratique pour le calcul sur 
ordinateur de P(x), i £ R fixe, oil P(x) est le polynome d 'interpolation de 
Lagrange qui verifie (4.1 1. 

Pour tout < i < j < n, on note par Pij(x) le polynome d'interpolation de 
Lagrange de degre inferieur ou egal a j — i tel que : 

Pi,i{x k ) = f k , t<k<j (4.3) 

En particulier le polynome P,,, n'est autre que la constante /, et le polynome 
PoM = P(x) . 



PROPOSITION 4.2.1 Pour tout < i < j < n, on a 

Pi A*) 



(x - Xi)Pi+ij(x) - (X - Xj)Pij-\(x) 



DEMONSTRATION: Soit < i < j < n, d'apres l'unicite du polynome 
d'interpolation de Lagrange, il suffit de prouver que le polynome : 

= (x - Xi)P i+hj {x) -{x- X^PjJ^ljx) 



est un polynome de degre inferieur ou egal a j — i et verifie (4.3 1. Le degre de 
Pi+i,j(x) est inferieur ou egal a j — i — 1 et le degre de Pij-i(x) est inferieur ou 
egal hj — i — l, ce qui donne que le degre de R(x) est inferieur ou egal a j — i. 
D' autre part, par definition on a : 



et 



d'ou 



Pi+i,j(x k ) = Pij-i(x k ) = f k Vi + 1 < k < j - 1 

Pi+lji^j) fj ) Pi,j — l\^i) ft 

R{x k ) =.f k Vi + 1 < k < j - 1 



R( Xi ) = x 'j) p ^-i^) = j. 

X •! " X 2 

R{ Xj ) =kzWiyM = /i 

x j x 1 1 

Soit x € K fixe. Pour calculer P(x) = P 0i „(x), l'algorithme de Neville est 
constitue de n etapes : 

Etape 1 : On calcule tous les Pi^+\{x), i = 0, • • • , n — 1 par la formule 

D / — \ _ — Xi)fi+1 — (x — X i+ i)fi 

— 

Xi-\-\ X{ 

Etape 2 : On calcule tous les P, j + 2(cc), i = 0, • • • , n — 2 par la formule 

(S - Xj)P i+ i )i+ 2(S) - (S - Xi +2 )Pi,i+l(S) 

P M+2 (x) = 

Xi+2 - Xi 
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Etape n—1: On calculc tous lcs P iii+ „_i(x), i = 0, 1 (sculcment deux 
valours a calculcr) par la formulc 

_ (x - x^P t+M+ „_i (g ) - (x - x i+ra _i)P i)i+n _ 2 (x) 

*i,i+n— 1 (^J 

Xi-\-n— 1 

Etape n : Une seule valeur a calculcr Pi, i+ „(x), « = par la formule 
p _ (x - xq)Pi^(x) - (x - x n )P , n _i(s) 

On peut resumer ces etapes dans un tableau, par exemple dans le cas n = 3 










1 


2 


3 


x 


Jo 


> 








Xl 


A 


> 


Pi A*) 


> P ,2(S) 


> Po,s(S) 


x 2 




> 


P2,3 (X) 


> A, 3 (x) 




x 3 


h 











EXEMPLE: Soit lc polynome P(x) = x 3 - 2x 2 + x + 1. On veut calculer 



! — 


P(3) = 


13 : 









1 


2 





1 








> 


1 




1 


1 




> 3 




> 


-3 




2 


3 




> 21 




> 


-27 




3 


13 







> -3 = P(-1) 



4.3 Methode des differences divisees 

L'algorithmc dc Neville permet de calculer P(x) pour une valeur donncc 
x — x. Si on veut calculcr P(x) pour plusieurs valeurs de x, on doit appliqucr 
l'algorithme de Neville plusieurs fois. Dans ce cas, il vaut micux calculcr toute 
la fonction polynome P(x). 

Soient (n + l)-couples reels ou complexes (x,, fi),i = 0, • • • , n, avec x, ^ Xj 
pour j et soit P(x) le polynome de degre inferieur ou egal a n d'interpolation 
de Lagrange tels que : 

P{xi) = fi, i = 0,- •• ,n 



REMARQUE 4.3.1 // est facile de voir que la famille : l,(x — x ),(x — 
Xq)(x — Xi),- ■ • ,(x — x )(x — xi) • • • (x — x„_i) est une base de I'espace vectoriel 
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des polynomes de degre inferieur ou egal a n, meme sans I'hypothese : xi ^ Xj 
pour i =/= j . En effet, le nombre de vecteurs est egal a la dimension de Vespace 
vectoriel, il suffit done de prouver que cette famille de vecteurs est libre. Soient 
ao, (Xi, • • ■ ,a n , (n + 1) -reels tels que : 

a + ai(x — x ) + • • • + a n (x — x )(x — x\) ■ ■ ■ (x — x n -i) = 

C'est une egalite entre polynomes, done, elle est vraie pour tout x € K. En 
particulier, si on remplace x = xq on obtient que ao — puis on divise par 
[x — xq) la nouvelle egalite et on remplace x = x% on obtient que ct\ — 0, etc..., 
jusqu'a a n — 0, ce qui prouve que la famille en question est libre. 

Par consequent tout polynome de degre inferieur ou egal a n s'ecrit : 

P(x) = a + a\{x — x ) + a 2 (x — x )(x — x\) + ■ ■ ■ 

+a n (x - x a )(x - xx) ■ ■ ■ (x - x„_i) 

La methode des differences divisees permet de calculer les coefficients a^, k — 
0, • • • , n du polynome P(x) avec le minimum d'operations et pour calculer P(x) 
pour une valcur x = x donnee, on utilise le schema de Horner : 

P(x) = (((• • • (a n (x - x n -i) + a n -i)(x ~ x n - 2 ) 

+a n - 2 ) ■ ■■) + ai)(x - x ) + a Q 

On rappelle que, pour < i < j < n, Pi j{x) est le polynome de degre 
inferieur ou egal a j — i d'interpolation de Lagrange qui verific : 

Pi,j(xk) = fk, k = «,••• ,j 



THEOREME 4.3.1 Pour tout < i < j < n, le polynome Pij(x) est egal a 
Pi,j{x) = Ci t i + c it i + i(x - Xi) + Ci^ i+2 (x - Xi)(x - x i+ i) H 

+ Cij(x — Xi)(x — Xi + i) ■ ■ ■ (x — Xj_i) 

ou les coefficients Cij sont definis par 

fi: i , • • • ; Tl 

<i < j <n 



En particulier 



P{x) = P ,„(x) = c .o + Co,i{x - xq) + c 0t2 (x - x )(x - Xi)+ 
h ca :n (x - x )(x - xx)- - • (x - x n -x) 



DEMONSTRATION: D'apres la Proposition 



4.2.1 



Pi,i(x) = f u i = (),••• ,n 



et pour tout < i < j < n 



(x - Xi)P i+ xj(x) -(x- x j )P iJ _ 1 (x) 



(4.4) 
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d'ou pour tout < i < j < n 

Pi,j{x) = —^Pi+x^x) - X _ X " Pij-i{x) + Pi,j-l(x) 
= Pij-^x) + Ri,j(x) 

avec 

RiA x ) = — (fl+uW ~ Pij-iW) ( 4 -6) 

"j J **J % 

On verifie facilement que le degre de Rij(x) est inferieur ou egal a j — i et que 
Rij(xi) = R itj (x i+1 ) = ■■■ = Rtj^j-i) = 

d'ou 

Rij(x) = Ci,j{x - Xi)(x - x i+ i) ■ ■ ■ (x - Xj-i) (4.7) 

ce qui donnc 

Pij(x) = Pij-x[x) + Ci,j(x - Xi)(x - x i+1 ) ■ ■ ■ (x - xj-i) (4.8) 

On applique le meme raisonnement a Pij-i(x) puis a Pi t j-2(x), etc..., on ob- 
tient : 

Pij(x) = Ci ti + Ci ti+1 (x - Xi) + Ci,i +2 {x - Xi)(x - x i+ i) H 

+Cij{x - Xi)(x - X i+ i) ■■■(x — Xj-l) 

II est clair que Ci^ = ft (il suffit de remplacer x par Xi). D'apres les formulcs 



(4.5) et (4.7 1, le coefficient du monome x 3 1 du polynome Pij(x) est egal a Cij 



car le polynome Pi j_i(x) est de degre inferieur ou egal a j — i — 1. D'apres la 



formule (4.4), le coefficient du monome x 3 1 du polynome Pij(x) est egal a 

X j x^ 

d'ou la relation, dite des differences divisees, entre les coefficients Cij. 

Dans la pratique, on utilise unc methode analogue a la methode de Neville 
pour calculer les differences divisees : 

1 2 ■■■ n 

x o fo = c 0,0 

> 

Xl fl > 0^2 

> Ci, 2 \ 

: : : : : : co, n 

> C n ~2,n-1 
Xn— 1 fn-1 > 

•£n fn 

Done, on lit les coefficients de P(x) sur la diagonale superieure. 

EXEMPLE: On cherche le polynome d'interpolation de Lagrange par la methode 
des differences divisees pour approcher la fonction 
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sur l'intervalle [—2,2] en utilisant les couples : (— 2, /(— 2)), (— 1.5, /(— 1.5)), 
(-1, /(-!)), (-0.5, /(-0.5)), (0, /(0)), (0.5, /(0.5)), (1, /(!)), (1.5, /(1.5)), (2, /(2)). 





n 
U 


l 


g 


j 


A 
H 


b 


b 


1 Q 
f □ 


-2.0 


0.200 


n 91 c 














-1.5 


0.308 


U.JDb 


0.169 


n mi 
U.Uj I 










-1 n 

I . u 


n 5nn 

U. JUU 


0.600 


U. L 1 J 


-0.277 


U. I 31 


0.037 






-0.5 


0.800 


0.100 


-0.200 


-0.100 


-0.062 


0.185 


0.019 


-0.019 


0.0 


1.000 


0.-100 


-0.800 


0.100 


0.100 


-0.185 


-0.123 


0.025 

0.019 


0.5 


0.800 


-0.600 


-0.200 


0.277 


-0.062 


-0.037 


0.019 




1.0 


0.500 


-0.385 


0.215 


-0.031 


-0.151 








1.5 


0.308 


-0.215 


0.169 












2.0 


0.200 

















P(x) de Lagrange 
f(x)=1/(1+x A 2) 




Lc theoreme suivant sera utilise dans le prochain paragraphe pour generaliser 
la methode des differences divisees. 

THEOREME 4.3.2 Soit f : R — ► K une fonction suffisamment derivable 
telle que : fk = f(xk), < k < n. Pour tout < i < j < n, on a : 
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oil I est le plus petit intervalle de K contenant Xi, Xi+i, ■ ■ ■ , xj et les Ci_j sont 



donnees par le theoreme 4-. 3.1 



DEMONSTRATION: On note par v id (x) = (x - x l ) ■ ■ ■ (x - Xj-i). D'apres 
le theoreme 14.1.21 : 

Vx g R,3c e I/f{x) = P id -i(x) + v itj (x) f ^. 

avec I est le plus petit intervalle de R contenant En particulier 

pour x = Xj, on obtient : 

/W-0( c ) 

3c G I/f{Xj) = I', , i ir ,*! + V it j{Xj)—^ 7\~|~~ 



(J-*)! 



or, par definition du polynome Py, on a : 

Pij{Xj) = fj = f(Xj) 



et d'apres la formule (4.8) 



d'ou 



J(i-i){c) 



c ij v ij( x j) — v ij( X j)~ 



or Vij{xj) ^ 0, ce qui prouve le theoreme. 



4.4 Interpolation de Hermite 



Soicnt /• , a«, i = 0, 1, • • • , m et fc = 0, 1, • • • , ri, — 1, des nombres reels tels 



que 



ao < ai < • • • < a r , 



Le probleme d'intcrpolation de Hermite consiste a trouver un polynome P{x) 
de degre infer ieur ou egal a n = — 1 + YlT=o ni verifie : 



P ( - k \a i ) = /* V0 < i < m,V0 < fe < m - 1 



(4.9) 



THEOREME 4.4.1 // exist e un unique polynome de degre inferieur ou egal 
a, 



a n = — 1 + ^""^ rij qui verifie (4-9) 



DEMONSTRATION: On note par R n [X] l'espace vectoriel des polynomes 
de degre inferieur ou egal a n. Soit C l'application lineaire de K„[X] dans K n+1 
definie par : 

C{P) = {P{a ),P'{a ),--- .P^Koo),--- , 

P{a m ),--- > P("™- 1 >(a m )) 
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Montrons que C est injective. Soit P(x) un polynome de degre inferieur ou egal 
a n tel que C(P) = 0, alors le polynome P(x) verifie : 

V0 < k < m - 1, P {k) {oti) = V0 < i < to 

Par consequent, pour tout i — 0, 1, • • • , m, (Xi est une racine de P(x) de mul- 
tiplicitc rii. Si on comptc lcs racincs avcc lcurs multiplicitcs, on obticnt n + 1 
racines , or P(x) est de degre inferieur ou egal a n, done, necessairement P = 0. 
D' autre part la dimension de l'espace vectoriel K n [X] est egale a n + 1, par 
consequent, C est bijective, d'ou l'existence et l'unicite de P £ M. n [X] qui verifie 



(4.9) 



Si on suppose qu'il existe une fonction (n+ l)-derivable / telle que pour tout 
i = 0, • • • , to et tout k = 0, • • • , le nombre /* = f^(ai), alors, le theoreme 
suivant permet d'estimcr l'erreur commise en approchant /(x) par le polynome 
de Hemite P(x) : 

THEOREME 4.4.2 Pour tout x £ [a, b] il existe un unique £ £ I tel que : 

v(x)f^ +1 HO 



f(x) = P(x) + 



(n+1)! 



avec v(x) = (x — ao) n °(x — ai)™ 1 • ■ • [x — a m ) nm et I est le plus petit intervalle 
contenant x, ao, a%, ••• , a m . 

La demonstration du theoreme precedent est exactement la meme que la 
demonstration du theoreme |4.1.2| du meme chapitre. 

Dans ce qui suit, on se propose de generaliser la mcthode des differences 
divisees pour determiner le polynome d 'interpolation de Hermite. Pour ccla, on 
ecrit les couples (an, i = 0, • • • , m, k = 0, • • • , — 1 sous la forme : 

(x J ),(x 1 ,f 1 ),--- ,(x n ,f n ) 



avec 



Xi = a , fi = fo i = 0, l,--- ,n -l 
Xi+n a = ax, f i+no = fl i = 0, 1, • • • , n x - 1 

Xi+no+m = 012, fi+n a +m = fl « = 0, 1, • • • , n 2 — 1 



Ji+n— n m +l fm ^ , 1 , • • - , n r 



REMARQUE 4.4.1 Pour tout < i < n, il existe un entier < r(i) < m 
tel que : Xi — a r ^y 

En utilisant la remarque 4.3. 1[ on ecrit le polynome d'interpolation de Her- 
mite P(x), associe aux couples (oti, ff) sous la forme : 

P(x) — a + ai(x — x ) + a 2 (x — xq)(x — x\) + ■•• 

+a n (x - x Q )(x - xi) ■ ■ ■ (x ~ x n -\) 

Le but est de calculer les coefficients ao, ai, ■ ■ • ,a n . On ne peut pas appliquer la 
mcthode des differences divisees car on n'a pas que Xi ^ Xj pour i ^ j. Soient 
£i, i = 0, 1, 2, • • • ,n, des nombres reels tels que : 

£o < a < • • • < fn 
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Le polynome d'interpolation de Lagrange de degre < n associe aux couples 
P(£i)), i = 0, 1, • • • , n coincide avec le polynome P(x) grace a l'unicite et la 
mcthode des differences divisees nous donne alors : 

P(x) = rf ,o + d ,i(x - £ ) + d , 2 (x - £ )(x - £1) H 

+do. n {x - £ )(x - Cl) • • • (X - 

avec 

di,i = P(£i),i = 0,1, • • • ,n 

et 

>/,,, ' / "'; ; '!.-' '. 0<i<j<n (4.10) 
si 4 « 

THEOREME 4.4.3 Pottr tout < z < j < n, la limite de dij lorsque £, tend 
vers Xi iend vers Xi+i,...et £j tend vers Xj existe et si on note cette limite 

par Cij on obtient 

P(x) = c , + co,i(x - x ) + c , 2 (x - x )(x - X\) + 

h C 0iTl (x - X )(x -Xl)"'(x- X n -l) 



et 



frii) * = 0,l,-",n 

fj-i 
•> r(i) 



{j-iY- , < i < j < n 
— , sinon 

X j " ~ 

avec < r(i) < m tel que ct r u\ = x\. 

DEMONSTRATION: Montrons par recurrence sur k = j — i avec k — 
0, 1, • • • , n, que les limites des dij existent. Pour k — 0, d i:i = P(£*), i — 
0, 1, • • • , n, done, par continuite de la fonction polynome P(x) on obtient que dij 
admet unc limite lorsque £j tend vers Xi et cette limite est egale a P(xi) = f®,^ 
avec < r(i) < m tel que a r ^ — x%. Supposons que les limites des di >s existent 
pour tous les I, s tels que < s — I < k et montrons que les limites des dij pour 
les i,j tels que j — i = k + 1 existent. En effet, d'apres le theoreme |4.3.2| et la 
formule (4.101, on a : 

d 



'•j 



di+ij — dij-i 
pO-i)( c ) 



< i < j < n 



avec c G I = le plus petit intervalle de R contenant £j, ■ • • , £j. D'apres 

l'hypothese de recurrence (resp. d^j-i) admet une limite. Par consequent, 

si Xi Xj la limite de dij existe et on a : 



D' autre part, SI X .j — X n tous les • • • , £j tendent vers la mcmc limite n 

et par suite c tend vers xt car l'intervalle / a la limite est reduit a {x{\. Dans ce 
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cas on utilise la continuity de l \x), on obticnt quo la limite de dij existe 
et on a : 



P^iXi) _ fr(i) 



(j-iY- (j-i)l 



avec < r(i) < m tel que a r (i) = Xi. 

EXEMPLE: On chcrchc lc polynomc d'intcrpolation dc Hermite par la methode 
des differences divisees generalisee pour approcher la fonction 



i 



l + x 2 



sur l'intervalle [-2,2] en utilisant les couples : (-2, f(-2)), (-2, f (-2)), (-1, f(-l)), 
(-1, /'(-!)), (0,/(0)), (0,/'(0)), (1, /(!)), (!,/'(!)), (2,/(2)), (2,/'(2)). 
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0.000 




0.000 




0.000 
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P(x) de Hermite 
f(x)=1/(1+x A 2) 

1 

0.8 
0.6 
0.4 
0.2 


-2 0. 2. 



4.5 Polynomes de Tchebychev 

Soicnt x € [—1,1] et a — Arccos(x) € [0,7r]. Pour tout rt € N, on a : 

cos(na) + isin(na) = (cos(a) + isin(o;)) n 

= (cos(Arccos(a)) + i sin(Arccos(a))) n 
= (x+ Wl - x 2 ) n 

La partie reelle de cette egalite nous donne : 

cos(nArcos(x)) = x n + C 2 n x n - 2 (x 2 - 1) + C^x n ~\x 2 - l) 2 + • • • 

DEFINITION 4.5.1 Pour tout n G N, on appelle polynome de Tchebychev 
de degre n le polynome T n (x) tel que sa restriction a Vintervalle [—1,1] est egale 
d la fonction cos(nArccos(x)) . 

REMARQUE 4.5.1 En utilisant la propriete : cos[(n + 2)a] + cos(na) = 
2cos[(n + l)a]cos(a) , on obtient la relation de recurrence : 

T n+2 (x) = 2xT n+1 (x) - T n {x) (4.11) 

et de cette relation, on verifie facilement par un raisonnement par recurrence 
que le coefficient du mondme du plus haut degre du polynome T n est egal a 2 n ~ 1 . 
On a egalement : 

T n (cos(0)) = cos(n0). 
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Le tableau suivant nous donne les polynomes de Tchebychev T n pour n = 
a 7 : 



n 


T n (x) 





1 


1 


X 


2 


2x z - 1 


3 


4a; 3 - 3x 


4 


8a; 4 - 8a; 2 + 1 


5 


16x 5 - 20x a + 5a; 


6 


32x B - 48a; 4 + 18a; 2 - 1 


7 


64a; Y - 112a; 5 + 56x d - 7x 



THEOREME 4.5.1 Pour tout n e N, le polynome de Tchebychev T n de degre 
n possede n racines : 

2k + 1 

Xfe = cos( — 7r), k = 0, 1, • • • , n — 1. 

2n 

£a fonction T n sur [—1,1] possede un extremum local en (n + 1) -points : 

k 

x' k = cos(— 7r), k = 0, 1, • • • , n. 
de plus, pour tout k = 0, 1, • • • , n, T„(a4) = (— l) fe . 

DEMONSTRATION: II est clair que 1 > x a > x x > ■ ■ ■ > x n > -1 et en 

utilisant la propriete T„(cos(#)) = cos(n#), on obticnt : 

2k + 1 7T 
T„(x fe ) = T„(cos(^^7r)) - cos[(2fc + l)-] = 

D' autre part, : 

Ti 

T' n (x) = (cos(nArccos(x)))' — sm(nArccos(x)) 

VI — x 2 

d'ou 



V* 


-K) 2 

n 


V* 


-K) 2 

n 




-K) 2 



sin(nArccos(a4)) 
sin(n^4rccos(cos(^7r))) 
sin(fc7r) = 



et on a : T n (x' k ) = T„(cos(|tt)) = cos(fcTr) = (-l) fe . 

PROPOSITION 4.5.1 Pour tout n e N : 

max |T„(x)| = 1 
xe[-i,i] 

DEMONSTRATION: 

On a : |T" n (x)| = \cos(nArccos(x))\ < 1 pour tout x e [—1,1] et T n (xk) — 
(-l) fc - 
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4.6 Meilleure approximation au sens de Tche- 
bychev 

Soit / une fonction continue sur 1 'inter valle [a, b]. On se propose de determiner 
les points Xi, i = 0, 1, • • • , n, pour lcsqucls lc polynome d 'interpolation de La- 
grange P n (x) associe aux (n+ l)-couples {x u yi)i=o, n , avec yi — f(xi), realise la 
meilleure approximation de la fonction /. D'apres le theoreme de l'estimation 
d'erreur de l'interpolation : 

/(»+*) ( C ) 

Mx £ [a,b],3c £\a,b[/f(x) - P n (x) = (x - x )(x - xi) ■ ■ ■ (x - x n )- — — — 

(n+lj! 



DEFINITION 4.6.1 On appelle meilleure approximation au sens de Tche- 
bychev de la fonction f par un polynome de degre < n , le polynome de l'inter- 
polation de Lagrange associe aux points d 'interpolation x,- L , i = 0, 1, • • • ,n pour 
lesquels la quantite : 

E(x ,xi, ■ ■ ■ ,x n ) = max \(x - x )(x - x{) ■ ■ ■ (x - x n )\ 

x£ [a, 6] 

est minimale. 

LEMME 4.6.1 Soit P un polynome de degre n tel que le coefficient de x n est 
egal a 1. Alors : 

max |P(a;)| > - — r 
xe[-i,i] 1 v n ~ 2™- 1 

DEMONSTRATION: On pose : 

= P{*) - 

Les points extremums (x' k )k=o,n de T n verificnt : 

1 = x' > x[ > ■ ■ ■ > x' n = -1 

Si on suppose que : 

3k £ {0, 1, • • • , n - l}/Q(x' k )Q(x' k+1 ) > (4.12) 
on obtient l'cxistcncc de k £ {0, 1, • • • , n — 1} tel que : 

Q(4) - P(x' k ) ^ > et Q(4+i) = P{x' k+ i) > (4.13) 

ou 

Q{x' k ) = P(x' k ) < et Q(x' k+1 ) = P(x' k+1 ) < (4.14) 

Si k est un entier pair, alors : 
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et si k est un entier impair, alors : 

Dans tous les cas, il existe i € {k,k + 1} tel que : 

1^)1 > i 



Si on suppose que la proposition (4.121 est fausse, alors : 

Vfc e {0, 1, • • • , n - l}/Q(x' k )Q(x' k+1 ) < 

Dans ce cas, la fonction continue Q(x) change de signe n-fois dans l'intervalle 
[-1,1], done, elle possede n-racines dans l'intervalle [—1,1]. D'autre part, Q(x) 
est un polynome de degre inferieur ou egal a n — 1, car le coefficient de x n du 
polynome de Tchebychev T n est egal a 2™" 1 . D'oii : 

Vx e R, Q(x) = P(x) - ^ = 

ce qui donne : 

max |P(a;)| = max | "^" , 1 = — — r- 
Dans les deux cas, on a ce qu'il faut pour demontrer le lemme. 

THEOREME 4.6.1 La meilleure approximation de la fonction f au sens de 
Tchebychev est realisee par le choix suivant : 

a + b b - a 
Xi = — 1 — s i: i = 0, 1, • • • ,n 

ou Si, i = 0, 1, • • • , n sont les racines du polynome de Tchebychev T n+ \ de degre 
n + 1, a savoir : 

,2i + l , 
Zn + 2 

DEMONSTRATION: Soient x , x u --- , x n , (n + l)-points de l'intervalle 
[a, b] et P n (x) le polynome d'interpolation de Lagrange associe aux (n + 1)- 
couples (xi,yi)i = o. n , avec yi = f(xi). II est clair que l'application : 

[-1,1] — [a, b] 
* ' t — b -^t+ *±s 

est une bijective et que : 

max |/(a;) - P„(x)| = max \f(<p(t)) - P„(p(t))| 

xG[a,b] 1,1] 

Si on note par = i = 0, 1, • • ■ ,n, Q n (i) = P n (cp(t)) et ,g(i) = f(<p{t)), 

on verific facilement que le polynome d'interpolation de Lagrange associe aux 
(n + l)-couples (^,j/i)i=o,n, avec yi = g(U) = f(xi), est egal a Q„(t). Par 
consequent, la recherche des points (xi)i = o,n qui realisent la meilleure approxi- 
mation de / sur l'intervalle [a, b] est equivalente a la recherche des points (U)i = o. n 
qui realisent la meilleure approximation de g sur l'intervalle [—1, 1]. 
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Etant donne des points (ti)i=o,n de l'intervalle [—1, 1], lc polynome 
P(t) = {t-t ){t-t x )---{t-t n ) 
est de degre n+1 et le coefficient de t n+1 est egal a 1, d'apres le lemme precedent : 

1 T 1 (t) 

max \(t - t )(t - h) ■ ■ ■ (t - t n )\ > — = max I "+ U ' \ 

D' autre part, on a : 

T n+1 (t) = 2 n (t - s )(t - si) • • • (t - s„) 

avec (si = cos(^q^7r))i = o.n les (n + l)-racines du ploynome de Tchebychev 
T„+i. D'ou : 

max |(i - t )(t - ti) ■ ■ ■ (t - t n )\ > max \(t - s ){t - s x ) • • • (t - s n )\ 
te[-x,i] «6[-l,l] 

Ce qui prouve que les points (sj)j = o,n realisent la meilleure approximation de g 
et par consequent les points : 

a + b b — a 
Xi = ip(si) = — 1 — Si, i = 0, 1, • • • ,7i 

realisent la meilleure approximation de /. 

4.7 Approximation au sens des moindres carres 

4.7.1 Cas general 

Etant donne u>i G]0,+oo[, i = 1, 2, • • • , m, on definit le produit scalaire sur 
M m comme suit : 

m 

yeR m ,ze M m , < y, z >= ^ ^y,z 4 = y T Dz 

t=i 

avec I? = diag(u>i, u>2, ■ ■ ■ ,ui m ) et on note par ||.||d la norme associee a ce 
produit scalaire : 



V 6 K m , Hylic = v<y,y > = VtFdv 

Soient A une matrice reelle m x fc et y S M m . On considere le probleme de 
minimisation : 

/= ™ in 4 Jly- z llc (4-i5) 



REMARQUE 4.7.1 II est connu que le probleme de minimisation (4- 15) ad- 
met un point minimal unique (un point qui realise le minimum de I) zq — Axq € 
Im(A) qui est egal a la projection orthogonale de y sur le sous-espace vectoriel 
Im{A). Le vecteur xo est unique si et seulemnet si la matrice A est injective (ce 
qui est equivalent a : les colonnes de A sont lineairement independantes). 
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THEOREME 4.7.1 Etant donne A une matrice m x k reelle et un vecteur 
y £ R m , alors : 

1. zq = Axo, xo £ M. k , est un point minimal du probleme de minimisation 
(4-15) si et seulement si xq est une solution du systeme lineaire : 

A T DAx — A T Dy (4.16) 



en particulier, le systeme lineaire (4-16) admet au moins une solution et 
cette solution est unique si et seulement si A est injective. 

2. la matrice symetrique A 1 DA est toujours positive et elle est definie posi- 
tive (en particulier inversible) si et seulement si la matrice A est injective. 

DEMONSTRATION: 



1. Le probleme de minimisation (4.15 1 est equivalent a la minimisation de la 
fonction : 

x £ R k , F(x) = \\y - Axf D = (y - Ax) T D(y - Ax) 

qui est une fonction derivable et on a : 

x £ R k , F'(x) = -2A T D(y - Ax) 

D 'autre part, il est facile de verifier que la fonction F est convexe : 

AG [0,1], (x u x 2 ) £R k xR k , F{X Xl + (1 - X)x 2 ) = \\y - A(\ Xl + (1 - A)z 2 )||d 

- \\\{y - A Xl ) + (\ - \)(y - Ax 2 )\\l 
< (A||y- J 4a; 1 || I) + (l-A)||2;-^ 2 || D ) i 
<A|| I /-A B i||2, + (l-A)||i/-AF a ||i, 
= XF( Xl ) + (1 - \)F(x 2 ) 

(la derniere inegalite est donncc par la convexite de la fonction reelle a 
variable reelle : x — > x 2 ). La fonction F est minoree sur K fe , done son 
minimum global est atteint si et seulement si l'equation : 

F'{x) = 

possede une solution dans R k , ce qui est equivalent au systeme lineaire : : 

A T Dy - A T DAx = 

et d'apres la remarque precedente, le minimum global de F est atteint en 
un point Xo £ R k et ce point est unique si et seulement si A est injective. 

2. II est clair que la matrice A T DA est une matrice symetrique et on a : 

Va; G R k , x T {A T DA)x = {Ax) T D{Ax) = \\Ax\\ 2 D > 
et Ax 7^ pour tout x ^ si et seulement si A est injective. 
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4.7.2 Approximation au sens des moindres carres d'une 
fonction par un polynome 

A la suite d'une experience, on obtient des resultats qui sont en general sous 
la forme : 



i 


1 


2 


3 




100 









0.3 


0.1 




-0.5 


Vi 


-2.3 


1.6 


5.2 




3.1 



Dans cette experience, on refait une meme manipulation rt-fois ( n =100 dans le 
cas du tableau ci-dessus) et a chaque fois on introduit une valeur x*, i = 1, • • • , n, 
Xi =/= Xj pour i ^ j, et on mcsurc lc resultat yi. Le probleme d' approximation 
consiste a trouver une fonction g qui approche le micux possible lcs points 
experimentaux, c'est a dire, une fonction g telle que sa courbe passe le plus 
proche possible des points experimentaux (x^ 2/i)i=j. ... n - Soit V un sous-espace 
vectoriel de 1'espace des fonctions continues ( Vest en general de dimension 
finie) et uji g]0,+oo[, i = 1,2, ••■ ,m, l'approximation au sens des moindres 
carres consiste a determiner g e V qui realise le minimum de : 
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min^c^s^) 

1=1 



Vif 



(4.17) 



Si on prend V = p e N, alors g(t) — a + ait + • • ■ + a p t p est un 

polynome de degre < p et le probleme de minimisation de l'approximation au 



sens des moindres carres (4.17) devient 



min 7 Wjfdo + aiXi + a 2 x\ H V a v x\ - yA 1 



(4.18) 



avec x T — (a a\ ■ ■ ■ a p ) e E p+1 . 

Si on note par D — diag(u>i, • • • , u> n ) et par 

x\ \ 



Ap — 



( 1 xi x\ x\ 

2 x 2 



1 X2 Xn x\ 



\ 1 



/ °0 \ 

ai 

«2 



et y = 



V a P ) 

lc probleme de minimisation (4.181 s'ecrit sous la forme gencralc 



min || j/ - ApX|||, 



2/2 



(4.19) 



On applique, alors, les resultats du paragraphe precedent pour calculer les co- 
efficients du polynome g. 



Exercice - Montrer que si p = n—l, alors la matrice A n _\ du probleme de mi- 



nimisation de l'approximation au sens des moindres carres (4.191 est inversible. 
Dcduire que pour tout p < n — 1 , les colonncs de la matrice A p sont lincairement 
independantes. (Ind. On utilisera le theoreme d 'interpolation de Lagrange.) 
Reponse : Pour p = n — l, le polynome d'interpolation de Lagrange de degre 
< n — 1 associe aux couples (xj,j/j), i = 1,2, ••• ,n, realise le minimum de 
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(4.181 car il donne une somme egale a zero. Done, tout polynome de degre 



< n — 1 qui realise le minimum de (4.181 donne une somme egale a zero ct par 



consequent, il vcrifie les hypotheses du theorcme de Lagrange. D'apres l'unicitc 



du polynome d'interpolation de Lagrange, le probleme de minimisation (4.181 



possede un point minimal unique. D'apres le theoreme 4.7.1 la matrice A„_! 
est injective, done, ses colonnes sont lineairement independantes, ce qui prouve 
qu'elle est inversible car e'est une matrice carree. On remarque que les colonnes 
de la matrice A p associee au probleme de minimisation pour p < n — 1 sont les 
p premieres colonnes de la matrice A n _\ associee au probleme de minimisation 
pour p = n — 1 qui est une matrice inversible. 



REMARQUE 4.7.2 L'exercice precedent montre que dans le cas p < n — 1, 
le polynome de degre p qui realise V approximation au sens des moindres carres 
est unique et que pour p > n—1, le polynome d'interpolation de Lagrange associe 
aux couples (a;,-, yi)i=\. n est I'un des polynomes qui realise le minimum de (4.. 18). 



Exemple 1. On considere l'experience suivante : 

[(Au module i de coefficient on) — >(Salah obtient , j/i=note du module 



I 


1 


2 


3 


4 


Xi 


modulel 


module2 


moduleS 


moduleA 


Vi 


Vi 


Vi 


2/3 


2/4 



Le probleme est de trouver une note globale de Salah. On cherche une note x 
la plus proche possible de ses quatre notes : y±, j/2, 2/3 e t 2/4- L 'approximation au 
sens des moindres carres permet de calculer cette note x, il suffit de minimiscr 
les carres : 



miny^a,(a; - y l ) 2 
xGM l — J 

i=l 

(le coefficient represente le poids du module i). Done p = 0, n — 4 et 



A = 



( 1 \ 
1 

1 

V 1 J 



x = x 6 



1)2 
2/3 

V ^ / 



et D — diag(cti, 012, «3, 0(4) 



ce qui donne A T DA = otx +02 + 0:3 + «4 et A T Dy = ot\yx + otiDi + caVz + «42/4- 
D'apres ce qui precede, le point x £ K qui realise le minimum est la solution du 
systeme lineaire : 

A T DAx = A T Dy 



d'ou 



aiyi + a 2 y2 + ct 3 y 3 + a^y^ 



Ot\ + «2 + CY3 + 0:4 

ce qui donne la moyenne. 

Exemple 2. Lors d'une experimentation, on a obtenu les resultats suivants : 



i 


1 


2 


3 


4 


5 







0.1 


0.2 


0.3 


0.4 


Vi 


-0.5 


0.2 


0.6 


1.2 


1.3 



<S(> 



ct on sait que la courbe de la fonction associee a cette experience est une droite 
(y = ax + b). Lc problcmc dc minimisation de l'approximation au sens des 
moindres carres est commc suit : 



min > [iii — axi — bY 

i—1 



(4.20) 



(les u)i = 1). Dans cet exemple p = 1 et n = 5 et 

/ 1 \ 

1 0.1 

1 0.2 

1 0.3 

V 1 0- 4 J 

D'apres le paragraphe precedent, le vecteur x € 



A = 



, x = 



y = 



(4.201 est une solution du systeme lineaire : 

A T Ax = A T y 



( -0.5 \ 
0.2 

0.6 et D = I 
1.2 

V 1-3 / 

2 qui realise le minimum de 

(4.21) 



Un calcul simple nous donnc 

,/ . I 5 1 



A 1 A = 



1 0.3 



et A T y = 



2.8 
1.02 



d'ou, la solution du systeme lineaire (4.211 : a = 4.6 et b = —0.36. Done, la 



droite la plus proche des resultats experimcntaux dc cct exemple au sens des 
moindres carres est : 

y = A.6x - 0.36 

1.6 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



Exemple 2 ° 
y=4.6x-0.36. 







0.05 



0.1 



0.15 



0.2 



0.25 



0.3 



0.35 



0.4 
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Exemple 3. 

A la suite d'unc experience, on a obtenu les resultats suivants : 



% 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Xi 





0.5 


1 


1.5 


2 


2.5 


3 


3.5 


4 


4.5 


Vi 


3.1 


2.5 


2.4 


2.3 


2 


1.7 


1.6 


1.5 


1.2 


1.3 



et on sait que la courbe de la fonction associee a cette experience est de la 
forme : y = (3e~ ax avec a et (3 deux reels positifs. On remarque que Log(y) = 
—ax + Log(/3) = ax + b qui est une fonction affinc. Au lieu de chercher une 
fonction qui approche y, on cherche une fonction qui approche Log(y), ce qui 
revient a appliquer la meme methode de l'exemple precedent avec les resultats 
suivants : 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Xi 





0.5 


1 


1.5 


2 


2.5 


3 


3.5 


4 


4.5 


Log{yi) 


1.13 


0.92 


0.88 


0.83 


0.69 


0.53 


0.47 


0.41 


0.18 


0.26 



ce qui permet de calculer les nombres a et b et par suite a = — a et (3 = e . 



4.8 Exercices : chapitre 4 

Exercice 1 

1. Calculer le polynome d'interpolation de Lagrange p de la fonction cosinus aux 
points suivants : 

7T 7T 
XX = 0, X 2 = -, X 3 = - 

Tracer soigneusement les courbes representatives de cosinus et de P dans le 
meme repere . Quelle conclusion en tirer ? 

2. Calculer le polynome d'interpolation passant par les points 

(0,0), (1,3), (3,1), (5,2), (8,2) 
par la methode des differences divisees. 
Exercice 2 

Nous savons que l'erreur pour Interpolation de Lagrange de / aux points xo,xi 

est 

f"(£(x)) 

f(x) - p(x) = (x - x )(x - x-i) , x < x < Xl, 

si / est de classe C 2 sur [xo, x\\. Determiner la fonction £(.) explicitement dans le cas 
oil f(x) — 1/x, xo = 1, xi = 2 et trouver le maximum (respectivement le minimum) 
de£(.)sur[l,2]. 

Exercice 3 

On interpole la fonction sinus aux points equidistants Xj = jh, j > 0. 

1. Majorer l'erreur d'interpolation dans l'intervalle [a;»,a;i+i] lorsqu'on fait passer 
un polynome de degre 3 par Xi-i, Xi, x i+i , x i+2 . 

2. Quel h peut-on prendre pour que cette erreur soit < 10 -8 , pour tout i > 1 ? 
Exercice 4 
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1. Determiner explicitement un polynome p de degre 3 tel que : 

p(l)=4, p(2) = 5, p'(l) = 3, p'(2) = 2 

2. Trouver un prolongement de classe C 2 par un polynome de la fonction 



/(*) = 



(x 2 — l)sin(irx) si x > 1 



Exercice 5 

Soient / de classe C" 1 sur l'intervalle [0, 1] et xo = < xi < xi < ■ ■ ■ < x n = 1 
une subdivision de l'intervalle [0, 1], n € N*. On note par B^n+ipf] l'espace vectoriel 
des polynomes de degre < 2n + 1 et par Pf le polynome d'interpolation de Hermite 
tel que : 

P f (x i ) = f(x i ),P' f (x i ) = f'(x i ) Vt = 0,l,"- ,n. 
1. Montrer que l'application C : R271+1 [X] — > E 2n+2 definie par 



C(P) = (P(x ), P(xi), ■■■ , P(x n ),P'(x ),P'( Xl ),- • • , P'{x n )) 

est un isomorphisme. 
2. Deduire l'existence de (2n+2) polynomes de degre < 2n+l, Po(x), Pi(x), 
Qo(x), Qi(x), ■ ■ ■ , Q n (x), independants de / tels que : 

n n 

Pf{x) = £ f(Xi)Pi(x) + J2 f'(xi)Qi(x) 



3. On suppose que xo = et x\ = 1. Montrer que pour tout a; G [0, 1], il existe 
c G [0, 1] tel que 

= 2/i =\2 

-/ (4) ( C ) 



/(i) = P/(x) + 



x 2 (l-a :) 2 
24 



4.9 Corrige des exercices : chapitre 4 



Reponse 1 

1) xi — 0, fi — cos(a;i) = 1, xi = f , /2 = cos(a;2) = 0.5, xz 



f 3 = cos(a;3) = 








1 


2 





1 






7T 


1 


2ir 


-3 


3 


2 


-3 




7T 

2 





7T 





P(x) 
2) 



27T *^ 


^2 ^(^ 


3 > 




7T 2 ^ 











1 


2 


3 




4 








3 










1 


3 


-1 


-4 
3 


-23 






3 


1 


+ 3 


15 


-131 




i 

2 


8 


-19 
280 


1680 




5 


2 





-1 
10 








8 


2 
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P(x) = 3x - f x(x - 1) - ffx(x - l)(x - 3) - ^x(x - l)(x - 3)(x - 5) 

Reponse 2 p(x) = § - f , /"(x) = Jr => £(x) = v^, 

1) V:r £ [xi,x i+ i],3c €]xi_i, x i+2 [/ 

sin (4) (c) 

sin(x) — P>(x) = (x — Xi_i)(x — Xi)(x — x i+ i)(x — x i+2 ) — — 

Pi(x)=\e polynome de degre < 3 qui passe par (xi_i, sin(xj-i)), (xi, sin(xj)), (xj+i, sin(xi+i)), 
(xj+2, sin(xj+2)), avec Xj = i/i. Comme pour tout x £ [xi,x;+i] on a \x — Xj_i| < 2/i, 
|x — Xi\ < h , \x — Xj+i| < h, \x — Xj+2| < 2/i, on obtient : 

| S in(x)-P(x)| < ^ 

2) II suffit que h < 

Reponse 3 

1) 








1 


2 3 


1 


4 










3 




1 


4 




-2 






1 


3 


2 


5 




1 






2 




2 


5 







P(x) = 4 + 3(x - 1) - 2(x - l) 2 + 3(x - l) 2 (x - 2) 

2) /(0-) = 0, /'(0-) = 0, /"(0-) = 2, /(!+) = 0, /'(!+) = 0, /"(!+) = -4tt 








1 


2 


3 


4 5 


























1 


-1 



















1 

-2tt - 1 


1 











-2tt 


-2tt 


1 








-2tt 






1 














P{x) = x 2 - x 3 + x 3 {x - 1) - (2tt + l)x s (x - l) 2 

{x 2 e x six < 
x 2 -x A + x\x - 1) - (2tt + l)x 3 (x - l) 2 si0 < x < 1 
(x 2 — 1) sin(7rx)six > 1 

Reponse 4 

1) dim(K 2n+2 ) = dim(R 2n+ i[X]) = 2n + 2 et il est facile de verifier que £ est 
lineaire. II suffit done de prouver que C est injective : 

C(P) = P(xi) = P'(xi) = OVi = 0, 1, • • • , n 
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d'ou, pour tout i = 0,1, ■ ■ ■ ,n, Xi est une racine double de P, ce qui prouve que P 
possede 2n + 2 racines et comme d°P < 2n + 1, necessairement PsO. 

2) On note par (ej)o<j<2„+ila base canonique de R 2n+2 . II suffit de prendre : 

Vi = 0,1, ■■;n,P i (x) = C~ 1 (e i ) 
Vz = 0, l,--- ,n,Qi(x) = £ _1 (e i+n+ i) 
II est clair que les Pj et les sont independants de / et comme : 

£{Pf) =(Pf(X0),--- ,Pf(x n ),P' f (x ),--- ,P'f{Xn)) 

= (/(*(>),••■ ,/(4/V«),- ,/'(*»)) 

= ^-'(/(lo), • • ■ , /(*„), f'(xo),- ■ ■ , /'(*„)) 

= £ 1 (Er=o /(aJi)ei + £™=o /"O0 e ;+™+i) 
en utilisant la linearite de on obtient ce qu'il faut. 

3) On applique le theoreme 4.4.2. 
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Chapitre 5 

Methodes numeriques 
d'integration d'une fonction 
et d'integration d'un 
systeme differentiel 



5.1 Formule d'integration de Newton-cotes 

Soit / unc fonction continue d'un intervalle [a, b] dans R. Pour calculer, d'une 
manicrc approchec, la valeur de 



la methode de Newton-cotes consiste a utiliser un polynome d'interpolation P(x) 
pour approcher la fonction f(x) sur l'intervalle [a, b] et a prendre 



comme valeur approchec de 1 'integrate de f{x) entre a et b. 

On pose h = , n e N* et on considcrc la subdivision uniforme : 



Soit P(x) un polynome d'interpolation de degre infericur ou egal a n tel que : 





Xk = a + kh, fc = 0,l,---,n 



f{x k ) = P{xk) = fk, fc = 0,l,--- 



n 



on obticnt alors : 
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Si P{x) est le polynome d'intcrpolation de Lagrange, on obticnt : 

hj a {x-x k )v'{x k ) 

avec v(x) = (x — x )(x — x\) ■ ■ ■ (x — x n ). En utilisant le changement de variable 
x = a + ht, on obticnt 

- 1 t V ^ dx - 17 — dt 

ak hj a (x-x k )v'{x k ) X J Ji^k-i 

par consequent les nombres a k , k = 0, l,--- , n sont des nombres rationncls 
independants de la fonction / et de l'intervalle [a, b] , ils dependent seulement de 
n. Si / est la fonction constante 1 et a = 0, 6=1, on obtient que le polynome 
d'intcrpolation de Lagrange P est aussi egal la constante 1, d'ou : 



«1 n n 

/ dx = h fkak = h^a k 

i — n l — n 



k=Q k=0 

ce qui donnc 

n 
k=0 

car dans ce cas particulier h— 1/n. 
Pour n = 1 , on obticnt 

ao = a i = 1/2 



d'ou 



6 f(x)dx « £(/(a) + f(b)) = b -^(f(a) + /(&)) 



Pour n = 2, on obtient ao = 1/3, a\ = 4/3, a 2 = 1/3. D'ou 

* « ^(/( a )+4/(a + />) + /(&)) 

= _Z_(/( a ) + 4/(a + /0+/(6)) 



5.2 Methode du trapeze 

Soit JV £ N*. On considere la subdivision uniformc : Xi = a + ih, i = 
0, 1, • • • ,N, avec h = (b — a)/N. La methode du trapeze consiste a appliqucr 
la methode de Newton-cotes avec n = 1 sur chaque intervalle i = 

0, 1, • • • , N - 1. On obticnt 

,6 N-l r x i+1 

/ f{x)dx = ]T / « T W 
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avec 

JV-l 

2 



= H^r + f(a+h) + -- + f(a + (n - + ^) 

THEOREME 5.2.1 La methode du trapeze est d'ordre deux, c'est a dire : 

,6 

/ f(x)dx-T(h) = 0(h 2 ) 

J a 

Plus precisement, si f est suffisamment derivable et h = (b — a)/N, on obtient : 
3c e]a, b[ I J f(x)dx - T(h) = - b -j^h 2 f" (c) 

Demonstration : 

LEMME 5.2.1 On suppose que la fonction f : [a,b] — » E est de classe C 3 . 
Alors 

3c €]a,b[ I f f{x)dx = b —^(.f(b) + f(a)) ^-^/" (c) 

J a, — ■ 

Demonstration ( du lemme ) 

On definit la fonction g : [a, b] — > E par : 

= £ f{t)dt- X -^{f{x) + f(a)) 
Le calcul des derivees de g> nous donnc : 

g(a) = g'(a) = g"(a) = et g^(x) = -\(.r(x) + (x a)f^(x)) 
D'apres la formule de Taylor avec reste integrate, on a : 



9(b) = jf {b ^gW(t)dt 

= ~\[( [b-tfr{t)dt+ [ (b-tf(t-a)fW(t)dt] 

^ J a J a 

en integrant par parties le deuxieme terme, on obtient : 
g{b) = - l -j\b-t){t-a)r{t)dt 

D'apres la formule de la moyenne et le fait que (6 — t)(t — a) nc change pas de 
signc sur l'intervalle [a, b], on obtient : 

3c e]a, b[/g(b) = J\b t){t a)dt = ^^(b a) 3 
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Revenons a la demonstration du theoreme. 

fb N-i Xi+1 

/ f(x)dx — T(h) = J2( -f^ dx 

Ja i=Q Jx, 

_x J± ip i{f{xi+i) + f{xi))) 

N-1 
i=0 

avec q €]xi, [c]a, 6[, i = 0, 1, • • • , N — 1. 

D'apres la continuity de /" sur l'intervalle [a, 6] et le fait que : 

m in fix) < ±TM + ■- + /" (^-i) < max r{x) 

xe[a,b] W ~ TV - xe[a.b] K ' 



on obticnt 



d'ou 



3c e ]a , h[/ r^) + r^) + - + r(cN-i) = f , (c) 



Bc€]a,b[/ I f{x)dx-T{h) = - h —^h\r{c) 

5.3 Methode de Simpson 

Soit N G N*. On considerc la subdivision uniformc : xi = a+ih, i = 0, 1, • • • , 2N, 
avec /i = (b — a)/2N. La methode de Simpson consiste a appliqucr la methode de 
Newton-cotes avec n — 2 sur chaque intervalle [x2i,X2i+2], i = 0, 1, • • • , N — 1. 
On obtient 

/•f> JV— 1 pX2i+2 

/ /(*)<&; = £ / /(z)da; w 

Ja j=0 >/x2i 

avec 

= ^ o (/( a + 2 ^) + 4 /(« + ( 2i + + /(« + ( 2 * + 2 W) 

= | (/(a) + 4/(a + h) + 2f(a + 2h) + if (a + 3h) + 2f(a + Ah) 
+ ■■■ + 2f(a + 2(N - 2)h) + 4/(a + (2N - l)h) + f(b)) 

THEOREME 5.3.1 La methode de Simpson est d'ordre quatre, c'est a dire : 

r b 

/ f(x)dx - S(h) = 0(h 4 ) 

J a 

Plus precisement, si f est suffisamment derivable et h = (b — a)/2N 7 on obtient : 
3c e]a,b[ I J f(x)dx S(h) = -^feV (4) (c) 
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Demonstration : 

LEMME 5.3.1 On suppose que la fonction f : [a, b] — > M est de classe C 5 . 
Alors, il existe c G]a, b[ tel que : 

f(x)dx = *Z£(/( a ) +4/(^±*) + /(&)) - ^/ (4) ( C ) 
Demonstration ( du lemme ) 

On pose A = a+ B = ^ ct D = On definit la fonction h : [a, D] -» K 
par 

= - ^(/(^ - x) + 4/(£>) + /(S + x)) 

JA-x 6 

le calcul des derivees de h nous donne : 

h(a) = h'(a) = h"{a) = h^\a) = h^\a) = 

et 

h^{x) = - 2 -i,{x)- X ^ a ^{x) 

avec 

iP{x) = f (4) (A-x) + f (4 \B + x) 
D'apres la formule de Taylor avec reste integrale, on a : 

h(D) = J°Vz^ h V> m 

~kj a {D ~ t)4{t ~ a ^'^ dt 

en integrant par parties le deuxieme terme, on obtient : 

h{D) = ~ [ D [3(D - t) 4 - 4(D - a){D - tf]^{t)dt 

J a, 

Comme t £ [a, D], on a : 

3(L> - t) 4 - 4(D - a)(D - t) 3 < —(D - t) 4 < 
D'apres la formule de la moyennc : 

3d e]o, b[/h{D) = [ D [3(D - t) 4 - A{D - a)(D - tf]dt 

'2 J a 

Apres avoir calcule l'integrale, on obtient : 

3c 1&] a A IKD) = - ^(^) 5 
D'autre part, en utilisant la continuitc de et le fait que 

min f^\x) < hfW(A- Cl )+f^(B + Cl )) < max /< 4 >(s) 

x£[a,6] Z x£[a,b] 
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on obtient 

3c€KM/^ = / (4) (c) 

d'ou 

3ce]a,b[ I h(D) = -1(^)5/^) 
Revenons a la demonstration du theoreme. 



f(x)dx - S(h) 



£(/ /(^ 

„-_n ^a;2i 



i=0 rfa! « 

- ^e"^' (/(»») + 4/(^1) 

+/(*a+a))) 



iV-l 
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avec Ci €]a, b[, i = 0, 1, • • • ,N — 1. D'apres la continuite de /( 4 ) sur l'intervalle 
[a, 6] et le fait que : 

• f(4)/ N < / (4) (co) + - + / (4) (^-l) . .(4), x 

mm j y '(x) < — < max j x '{x) 

x£[a,b] N x£[a,b] 



on obtient 



d'ou 



3c e]a, 6[ / / (4) (co) + /^(c 1 ) + - + ^( C .- 1 ) = /(4)(c) 



< =]«•'<[ / ' f(x)dx - S(h) = '' |N|) "// ; / ; (c) 



5.4 Probleme de Cauchy 



Soit [a, 6] un intcrvalle ferme de K. Soit / unc application deRx 
On appelle systeme differentiel du premier ordre la relation : 



dans 



dy 
dx 



(x) = f{x,y{x)) 



(5.1) 



avec y : [a,b] — ► M. n une application derivable sur ]a,b[. On appelle probleme 
de Cauchy le systeme defini par les equations (5.1) et la condition initiale : 



y(x ) = y 



(5.2) 



avec Xq € [a, b] 
et j/o G K n - 

Une application y de [a, 6] dans 



est dite solution du probleme de Cauchy 



(5.1 )+( 5.2 1 si elle est derivable sur ]a,b[ et verifie les equations (5.1) et (5.2|. 
Nous admettons le theoreme suivant : 
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THEOREME 5.4.1 Soit / : R x E" — ► W l une application continue. Si f 
verifie la condition de Lipschitz par rapport a la deuxieme variable : il existe 
L > tel que 



f(%,y) ~ f(x,z) ||< L \\ y- z ||, Vy, z£ 



Vx €}a, b[ 



avec || . || est une norme quelconque de W 1 . Alors, le probleme de Cauchy 
(5.1 )+(5.2) adraet une solution et une seule sur [a,b] pour tout xq € [a, b] et 
yo G K n - 



COROLLAIRE 5.4.1 Soit f : R x R n — ► E" une application continue, 
derivable par rapport a y et verifie : 

i=l ^ l 



^4Zors ; /e probleme de Cauchy (5.1)+(5.2) admet une solution et une seule sur 
[a, b] pour tout Xq € [a, 6] et yo € R n . 

Demonstration : Soicnt zefctyer. Pour tout x e [a, b] et tout 

j = 1, 7i, on applique lc lcmme dc Rollc a l'application Hj(t) — fj(x, y + t(z — y)) 
sur l'intervalle [0, 1] : 

3c G ]0, 1[ / f,(x, z) - fi{x,y) = H 3 (\) - ^(0) - H'Ac) 



or 



ce qui donne 



< Ik -ylUE^il^y + *(«-»))! 

< Afllz-i/Hoo 



d'ou 



||/(a;,z) - /(Xj^Hoo = max|/j-(x,z) - < M|jz- j/H^ 

j 

il suffit done d'appliquer le theoreme precedent. 

Dans toute la suite de ce chapitre nous nous interessons a la resolution numerique 
du probleme de Cauchy. Nous supposons done que le probleme de Cauchy 
( 5. 1 )+( 5.2 1 admet une solution unique definie sur [a, b]. Etant donne une sub- 
division dc l'intervalle [a, 6], {xq,xi,--- l'idee principale des methodes 
numeriques est d'approcher numeriquement les valeurs {y(%o), y(%i), ■ ■ • , y(%k)} 
par des valeurs approchees {yo,yx, ■ ■ ■ ,Vk}i (y( x o) — Vo et y(-) est la solution 
exacte du probleme de Cauchy (5.1 1+(5.2|). 
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5.5 Integration approchee 



5.5.1 Principe d'une methode numerique 

Notons par xq = a et par yo la condition initialc du probleme dc Cauchy. A 
l'etape m, < m < k — 1, on pose : 

2/m+i = Dm + h m $(x m ,y mi h rn ) (5.3) 

ou $ est l'application qui caracterise la methode d'integration et h m = x m +\ — 
x m . Le choix de la fonction $ depend de la precission avec laquelle on veut 
approcher les vraies valeurs {y{x m )) m= i^k- Pour simplifier les notations, on sup- 
posera dans toute la suite de ce chapitre que la subdivision de l'intervalle [a, b] 
est uniforme, c'est a dire : 

Vm = 1, k, h m = h = — - — 



5.5.2 Consistance d'une methode numerique 



La consistance est une notion theorique qui exprime la relation entre le systeme 
differentiel et la fonction <f>. 



DEFINITION 5.5.1 On dit que la methode d'integration (5.3) est consis- 
tante avec le systeme differentiel (5. W si : 



(5.4) 



l im max || y&rn+l) V(x m ) _ Q 
h— >0 m n 



pour toute solution y du probleme de Cauchy (5.1) + (5.2) 



THEOREME 5.5.1 Une condition necessaire et suffsante pour que la methode 



d'integration approchee (5.3) soit consistante est que 

®(x,z,0) = f(x,z), Va;e[o,6], VzeK" 



(5.5) 



DEMONSTRATION : Montrons que la condition (5.51 est necessaire; en 



effet, soit x € [a, b] et z e W 1 , d'apres le theoreme 5.4.1 il existe une solution 
unique y du systeme ( |5.1| avec la condition y(x) = z. On considere la subdivision 
de l'intervalle [a, b] definie par 

x — a 

h = , x a = a 

P 

xx = xq + h, X2 = xi + h, ■ ■ ■ , x p = x, x p+ i = x p + h, 
■ ■ ■ ,x k -i = a + (k — l)h 

oil k est l'entier verifiant l'inegalite suivante : 

a + (k - l)h < b < a + kh. 



et on pose x^ — b. Utilisons le fait que y est une solution du systeme (5.1 1, d'ou 

f(t, y{t))dt = y{x m+ i) - y(x m ). 
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Si la methode est consistante, on obtient que 
1 [ x ™+i 

limmax||- / f(t, y(t))dt - ®(x m , y(x m ), ft) \\= 0. 



h— >0 m ft _ 

En particulier pour m = p 



lim || - / [/(t,i/(*))-$(i,y(i),fc)]dt||=0. 
fr->o ft 

Utilisons lc fait que les fonctions /, $ sont continues et que z — y(x), alors 

= $(x,z,0). 



Montrons que la condition ( |5.5| est suffisante; d'apres (5.5 1, la condition (5.3 1 
est la meme que la condition suivante : 



1 



lim max || - / [$(t, y(t), 0) - >S>{x m , y{x m ), h)]dt \\= 
fe^O m ft - / a , m 

or, cette condition est une consequence immediate de la continuite uniforme de 
<I> sur le compact {(t, y(t))/t € [a, b]} x [0, b — a]. 

5.5.3 Stabilite d'une methode numerique 

Soient (y m )m=l,ki { z m)m=i.k les solutions respectivement de : 



Vm+i = Vm + h$(x m ,y m , ft) 
yo £ M. n quelconque donne 

et de 

z m +i = z rn + h[$(x m , z m ,h) + e rn ] 
zq £ W n quclconque donne 
avec 

zo = yo 

Xq = a 

Xm+l — X m -\- ft 



(5.6) 
(5.7) 



DEFINITION 5.5.2 On dit que la methode d' integration (5.S) est stable s'il 
existe M\ independant de h tel que : 

max || y m z m || < M\ max || e m \\ (5.8) 

m m 

La notion de stabilite veut dire en "gros" qu'une petite perturbation de la 
fonction e'est a dire de la methode cntrainc au plus une petite perturbation 
sur le resultat. Comme les valeurs (y m ) sont calculees d'une maniere approchee, 
il faut que la methode utilisee soit stable, sinon on risque d'obtenir des resultats 
completement faux. 

THEOREME 5.5.2 S'il existe M > independant de ft tel que pour ft assez 
petit <I> verifie 

|| $(x,y,h) - $(x,z,h) \\< M || y-z ||, Vj/,zeM", Vx£[a,b]. 



Alors, la methode d'integration (5.3) est stable. 
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DEMONSTRATION : Pour simplifier la demonstration, on suppose que le 

pas de la subdivision (x m ) m —i & est uniforme, c'est-a-dire que h = . On a 

k 

II Vm+i ~ z m+1 || <|| y m - z rn || + 

h || ®(x m ,y m ,h) - $(x m ,Zm,h) \\ +h \\ e m \\ 

soit encore 

|| y m+ i - z m+ i || < (1 + Mh) || y m - z m \\ +h \\ e m \\ . 
Montrons par recurrence que 

(1 + hM)^ 1 - 1 
II y m +i - 2m+i || < 77 max || e p \\ . 

IVl p=l,k 

En cffct 

II Vm+2 - z m+2 || < (1 + Mh) || y m+ i - z m+1 || +h || e m+ i ||, 

d'ou 

(1 + hM)' m+1 — 1 
II Vm+2 - z m+2 || < [(1 + /iM) — h /i] max || e p ||, 

M p=l,k 



soit encore 



(1 + hM) m+2 - 1 

M p=T"fe 



Vm+2 - z m +2 || < ^7 max || £ p 



D' autre part, pour x > on a : 

1 + x < 

d'ou 

(1 + Mf) m < e mhM < e [b - a ^ M , Vm = 1, • • • , k. 

Par consequent 

max || y m - z m ||< M x max || e m ||, 

ra — 1 , fc 7n=l, K 

e (b-a)M _ l 

avec Mi = . 

M 

5.5.4 Convergence d'une methode numerique 



DEFINITION 5.5.3 On dit que la methode d'integration (5.3) est conver- 
gente si : 

lim max || y rn - y{x m ) ||= 0. 

h— >0 m 



THEOREME 5.5.3 Supposons que la methode d'integration (5.3) est consis- 
tante et stable. Alors, elle est convergente. 
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DEMONSTRATION : Soit y la solution du systeme differentiel ([H7TJ) avec 
la condition initiale y(a) = yo, on pose 

e m = y(a:m+l) ~^ m) -Hx m ,y(x m ),h). 
h 

Si la methode est consistante, par definition on a 

lim max II e rn 11= 0. 
h->0 m 

D'autre part, si on applique la definition de la stabilite a 

(-^m)m— l.k (y(*^m) )m— 1,/c 

et 

(y?n)m— l.k j 

on obtient : 

max || y m - y{x m ) ||< Mi max || e m ||, 

m m 

d'oii la convergence en faisant tendre h vers zero. 
5.5.5 Ordre d'une methode numerique 

L'ordre d'une methode est un moyen theorique pour mesurer la precision avec 
laquelle les valeurs approchees sont calculees. 

DEFINITION 5.5.4 On dit que la methode d 'integration est d'ordre p si : 
max || y m - y(x m ) \\= 0(h p ). 

rn 

oil y(.) est la solution du probleme de Cauchy ( 5. 1) + JX^| ) pour y quelconque. 



THEOREME 5.5.4 Supposons que la methode d 'integration (5.3) est stable 
et que 

max || ~[y(x m+ i) - y(x m )} - $(x m ,y(x m ),ti) ||= 0(h p ). 



pour toute solution y(.) du probleme de Cauchy (5.1)+(5.2). Alors, la method 



d'integration (5.3) est d'ordre p 



DEMONSTRATION : Soit y(.) la solution exacte du probleme avec 
la condition initiale y(xo) = yo- Soit (y m )m=i,k la solution de la methode 



d'integration (|5.3|. On a 



y{x m+ i) - y{x m ) = h[$(x m , y(x m ),h) + ( y(3: " ,+1 , ) i y< " Xm) - $(x m , y{x m ), h))] 
Dm+i ~y m = h$(x m , y m ,h) 

a- 1j{Xm+l) y(Xm) » > , \ , / 

Si on pose e m = 9(x m ,y(x m ),h) et z m = y{x m ), par 

h 

definition de la stabilite on a : 

max || z rn - y m ||= max || y{x m ) - y m ||< Mi max || e m \\, 

mm m 

il suffit done de prouver que e m — 0(h p ), ce qui est vrai par hypothese. 

La proposition suivante est tres utile pour le calcul de l'ordre d'une methode 
numerique. 
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PROPOSITION 5.5.1 Soient x £ [a, b] et z £ K avec K un compact de R n . 
Supposons que f et ^ sont suffisamment regulieres, alors : 

1. la solution y du probleme de Cauchy \ 5.1)+{5.2) verifie : 

h 2 h p 
y(x + h) = y(x) + hy'{x) + —y" (x) + ■ ■ ■ + -y {p \x) + 0{h p+1 ) 

2. la fonction $ verifie : 

$(ar, z, h) = ${x, z, 0) + h — (x, z, 0) + ■ • • + h p -^(x, z, 0) + 0(/^ +1 ) 
et dans les deux cas le teste est une quantite uniformement par rapport a (x, z). 



La demonstration est une application simple de la formule de Taylor avec reste 
integrale a chaque fonction h — > yi(x + h) et a chaque fonction h — > <&i(x, z, h), 
?'=!,••• , n. 



5.6 Methode d'Euler 



Soit y(.) la solution exacte du probleme (5.1l. La methode d'Euler est une 



consequence de la remarque suivante : pour h assez petit, on a : 

y(x + h)- y(x) ~ h f(x, y{x)) ; 

par consequent, pour h = ^jf 2 , en tout point x m — a + mh, m = 0, 1, • • • , k, on 
obtient une approximation y m de la valeur y{x m ) de la solution exacte y{x) de 
la maniere suivante : 

2/0 = y{x ) 

Ce qui revient a prendre 

$(x,y,h) = f(x,y), Vx,y,h. 

II est clair que la methode d'Euler est consistante et si on suppose que / est 
lipschitziennc par rapport a y uniformement par rapport k x £ [a,b], on obtient 
une methode stable et convergente. 

Exercice d'application : Verifier que la methode d'Euler est d'ordre 1 lors- 
qu'on suppose que / est lipschitzienne par rapport a y. 

Reponse : Soit h > 0. La fonction t — > f(t,y(t)) est continue sur [a,b], done, 
elle est uniformement continue sur [a,b], d'ou : 

max \\f(t,y(t))-f(t',y(t')\\=0(h) 

\t-t'\<h 



D'autre part, d'apres (5.1l : 



y(x m+ i) - y(x m ) = / f(x,y(x))dx 
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d'ou 



\\[y{x m +\) - y{x m )} -${x m ,y{x m ),h)\\ 

n 1 f x m+1 



,11 , (f(x,y(x)) - f(x m ,y(x m )))dx\\ 
n J x 

1 f x ™+ 1 

<t \\f{x,y{x))- f{x m ,y(x m )))\\dx 

< - / 0(h)dx = 0(h) 



5.7 Methodes de Runge-Kutta 

Pour obtcnir des methodes d'ordre superieur a un, on introduit les methodes 
dites de Runge-Kutta : 

h = f(x,y), 

k 2 = f(x + 6 2 h,y + a 2 ihki), 

h = f(x + 9 3 h 1 y + a 3 ihki+a 32 hk 2 ), 



ct on pose <I>(x, y, h) = ^ Cjkj. 

= 1 

II est clair que 

r 

$(x,y,0) = f(x,y)J2 c f, 

= 1 

par consequent, les methodes de Runge-Kutta sont consistantes si et seulement 
si 

r 
3 = 1 

Si on suppose que / est lipschitzienne par rapport a y uniformement par rapport 
axe [a, b], on verifie facilement que pour h petit <fr(x,.,h) est aussi lipschit- 
zienne par rapport a y uniformement par rapport a x G [a, b], avec une constante 
de Lipschitz independante de h et de x, ce qui prouve que les methodes de Runge- 
Kutta sont stables. Les methodes de Runge-Kutta sont done convergentes. 

Etudions l'ordre de quelques unes de ces methodes. Un calcul simple nous 
donne : 

y'(x) = f(x,y(x)), 

y „ (x) = ^ = «M + «M/( IlB)l 

ct dans le cas <&(x, z, h) = J^ i=1 Ciki : 

®(x,z,0) = (ci + c 2 )f(x,z) 

, c\\ r „ df(x,z) df(x,z) ,. 

^(x,z,0) = c 2 [0 2 ^ + a 21 ^f(x,z)} 
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ce qui donne, en utilisant la proposition |5. 5. l| : 

y^ +h l~y^ = f(x , y{x)) + h[ d -i^l + '^"iu. y)] + o^) 

et 

$(x,y,h) = cif(x,y)+c 2 f(x + 9 2 h,y + a 2 ihf(x,y)), 
= (ci+c 2 )f(x,y) + c 2 h[0 2 ^^- 

+a 21 W^f(x >y )]+0(h 2 ). 

D'ou, pour que la methode de Runge-Kutta soit d'ordre deux, il faut que 

ci + c 2 = 1, c 2 9 2 = 1/2, c 2 a 21 = 1/2. 

Une solution de cette equation est 

ci = c 2 = 1/2, 9 2 = a 21 = 1, 

ce qui donne la methode de Heun (1900) : 

<S>(x,y,h) = -[f(x,y) + f{x + h,y+hf(x,y))], 

elle demande le calcul de / seulement deux fois par iteration. Une autre solution 
est 

ci = 0, c 2 = 1, 9 2 = 02i - 1/2, 
elle est connue sous le nom "Methode d'Euler modifiee" [Collatz 1960] 

$(x,y,h) = f(x+^,y+^f(x,y)), 
La methode de Runge-Kutta d'ordre quatre est la suivante 

k 2 = f{x + h/2,y+hk 1 /2), 
k 3 = f(x + h/2,y + hk 2 /2) 
fc 4 = f(x + h,y + hk 3 ), 

$(x,y,h) = l (h + 2k 2 + 2k 3 + hi) . 
6 

On montre que la methode de Runge-Kutta 4 est d'ordre 4. 



5.8 Exercices corriges 

Exercice 5.1 

Soit a < xq < X\ < ■ ■ ■ < x n < b une subdivision de l'intervalle [a, b]. Soit P{x) 
un polynome de degre < n. 

1) Montrer que le polynome d 'interpolation de Lagrange Q de degre < n tel 
que : Q{xi) — P{xi) pour i = 0, 1, • • • , n, coincide avec P(x). 

2) Montrer que les coefficients du polynome P(x) verifient un systeme de la 
forme Ax = b avec A une matrice a determiner. Deduire du theoreme d'inter- 
polation de Lagrange que A est inversible. 
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3) En utilisant l'cxprcssion du polynome d'interpolation de Lagrange, montrer 
qu'il existe 70,71, • • • ,7n des constantes independantes de P(x) tels que : 



r-b n 

/ P{x)dx = Y J liP{^i) 



4) En rcmplacant le polynome P(x), dans l'equation des 7^ par les monomes 
l,x,x 2 , ■ ■ ■ ,x n , montrer que le vecteur forme par les 7* est une solution d'un 
systeme de la forme A T x — d ( la matrice A de la question 2)). En deduire que 
les 7j sont uniques. 
Reponse 5.1 

Soit Q(x) le polynome d'interpolation de Lagrange de degre < n tel que : 

Q{xi) = P{xi), i = 0,l,- 
On sait que l'expression de Q(x) est la suivante : 



. n 



avcc 



Q(x) = ^P(x 2 )L 4 (x) 

i=0 

Li(x) - 



rij^iC^i x j) 

D'autre part, d'apres l'unicite du polynome d'interpolation de Lagrange et le 
fait que degre de P(x) < n , on obtient que P(x) = Q(x). D'ou : 



/b n r-b 

P(x)dx = ^P(xi) I Li(x)dx 
i=o •' a 



On pose 7i = J a Li(x)dx pour i = 0, 1, • • • , n, ce qui donne l'existence des 7$. 
Montrons l'unicite des 7,. Le polynome d'interpolation de Lagrange P(x) — 
ao + aix+- ■ - + a n x n est unique, done les coefficients do, oti, • • • ,a n sont uniques. 
Par consequent, le systeme suivant admet une solution unique : 



a + aix + 
a + a\X\ + 



+ a n XQ 
+ ffl ra x™ 



a + a!X n -\ ha„x™ = P(x„) 

quclquc soit les nombres P(xq),P(xi), • • • , P(x„). D'ou, la matrice : 



( I x xl 
1 Xi xf 



'-o 
x\ 



\ 1 2-" S'n ' ' ' x n ) 

est inversible. D'autre part, si on rcmplace, dans l'equation des 7,, le polynome 
P(x) par les monomes 1, x, x 2 , • • • , x™, on obtient le systeme : 



7o + 7i H ^7« 

70^0 + 71^1 H 1" 7n£r, 



J 6 ldx 
f xdx 

J a 



loxfi + 71 xi H h 7„x™ = J o x"dx 
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qui s'ecrit sous la forme matricielle : 

A T 1 = b 

avec 7 = (70 7i • • '7n) T - Comme la matrice A est inversible, alors, A T est 
aussi inversible. Ce qui prouve l'unicite des 7$. 
Exercice 5.2 

Si / G C 2 [a, 6] alors il cxiste x e [a, 6] tcl que l'erreur dc la mcthode du trapeze 
est egale a : 

b f(x)dx \{b a)(f(a) + f(b)) =^(b- «) 3 /" (*) 

Dcduire ce resultat en utilisant un resultat du cours du chapitre precedent. 
Reponse 5.2 

D'apres un theoreme du chapitre precedent et pour n = 1, xq = a et x\ = b : 
V.x G [a, b],3c(x) € [a, b]/f(x) = P(x) + (^^31^1 

avec P(x) = f(a) + (f(b) — f{a)){x — a)/(b — a) le polynome d'interpolation de 
Lagrange de degre < n = 1 et qui verifie P(a) = f(a) et P(b) — f(b). Montrons 
que f"(c(x)) est continue sur [a, b\. On a : 

2{f{x)-Pix)) 
1 {C[X)) ~ (x-a)(x-b) 

II est clair que la fonction f"(c(x)) est continue sur [a, b] sauf peut etre en a et 
b, mais il est facile de voir qu'on peut la prolonger en ces points car : 

Jim /■<«.)) - ,11m f^Sf - ^</» - P») 

et 

lim r(c(x))= lim 2(/(x) -P(x)) = 2 
i^fe- k^6- (x — a) (2; — 0) — a 

On peut done supposer que f"(c(x)) est continue sur [a,b]. D'autre part, (x — 
a)(x — b) nc change pas de signc sur [a,b], d'apres la formule de la moyenne : 

3a e [0,6]/ / {x-a){x-b)f'{c{x))dx = /"(c(a)) / (a; - a)(x - b)dx 

J a J a 

En calculant P{x)dx et f (x — a)(x — b)dx, on obtient ce qu'il faut. 
Exercice 5.3 

Si / G C 4 [a, b] alors il existe x £ [a, 6] tel que : 



J a 



' f{x)dx -\{b- o)(/(a) + /(&)) - (6 ^ )2 (/'(a) - /'(&)) - %^f«\x) 



Deduire ce resultat du cours du chapitre precedent. 
Reponse 5.3 
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D'apres un theoreme du chapitre precedent ct pour m = 1, £o = a, £1 = b et les 
conditionsP(a) = /(a), P'(a) = f(a), P(b) = f(b), P'{b) = f(b) : 

Vx e [a, 6], 3c(x) G [a, &]//(x) = P(x) + (* ~ ") 2 (* -j0 2 / (4) (c(*)) 

avec P(x) = cqo + Qn(x — a) + c 2(x — a) 2 + 033(2; — a) 2 (x — b) le polynome 
d'interpolation de Hermite de degre < 3 et les coefficients donnes par la methode 
des differences divisees generalised : 

coo = /(a) 

coi = /'(a) 

f(b)-f(a)-(b-a)f(a) 
C02 = 



C03 



{b -a) 2 

/'(6)+/'(o)-2(/(6)-/(a))/(6-a) 



(6 -a) 2 

Montrons que /( 4 )(c(x)) est continue sur [a, 6]. On a : 

7 1 [ >> ~ (x-a) 2 (x-b) 2 

II est clair que la fonction /( 4 )(c(x)) est continue sur [a, b] sauf peut etre en a 
et b, mais il est facile de voir qu'on peut la prolonger en ces points car : 

lim ./( 4 >(c(x)) = lim 4!(/(x)-P(x)) = _?_ {f , {a) p» (o)) 

a;^a+ rr^a+ (x — a) ^ (x — 6) ^ (a — Oj 2 

et 

lim /< 4 >( C (x))= lim 4!(/W -P(x)) = 2^ 

x^b- x-> b- (x — a) z (x — by [b—ay 

On peut done supposer que f^ 4 \c(x)) est continue sur [a, b\. D'autre part, 
(x— a) 2 (x— b) 2 ne change pas de signe sur [a, b], d'apres la formule de la moyenne : 

3a e [a,b]/ f \x-a) 2 (x-b) 2 f^(c(x))dx = f {A) {c{a)) Ax - a) 2 (x - b) 2 dx 

En calculant P{x)dx et f (x — a) 2 (x — b) 2 dx, on obtient ce qu'il faut. 
Exercice 5.4 

Soient / <G C 6 [— 1, 1] et P e R 5 [X] lc polynome d'interpolation de Hermite avec 
P( Xi ) = f( Xi ), P'(xi) = fix,), Xl = -1,0, +1 



1) Montrer que pour tout Q € R 5 [X] on a : 
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2) La formula do la question prcccdcntc est unc mcthodc pour approchcr l'integrale 
de la fonction / sur [—1,1] par l'integrale du polynome de Hermite P(x) : 



L 



1 mdx * + 15/(0) + + - 



Lorsque / est un polynome de degre < 5 la formule est exacte. Montrer que 
cette formule n'est pas, en general, exacte pour les polynomes de degre < 6. 
3) Utiliser un resultat du chapitre precedent pour determiner une estimation de 
l'erreur de la mcthodc dc la question 1). 
Reponse 5.4 

1) II suffit de verifier la formule pour les monomes : 1, x, x 2 , x 3 , x 4 , x 5 , car 
l'integrale d'une combinaison lineaire est une combinaison lineaire des integralcs 
et le polynome est une combinaison lineaire des monomes. 

2) La formule n'est pas exacte pour x 6 . 

3) D'apres un theoreme du chapitre precedent : 

w e [-1.1,. Mx) e - P(x) + (^DVt,-!)'/"'^)) 

De la meme maniere que l'exercice precedent, on vcrifie que la fonction (c(x)) 
est continue sur [—1, 1] ct on a que (x + l) 2 x 2 (x — l) 2 ne change pas de signe 
sur [—1, 1]. D'apres la formule de la moyenne, il cxistc a G [— 1, 1] tel que : 

J x 2 {x + l) 2 {x-l) 2 f 6 \c{x))dx= / (6) (c(a)) J x 2 {x + l) 2 {x-l) 2 dx 

Tout calcul fait, on obtient : 

J\.f(x)dx =£/(-!)+ "/((J) + 

Exercice 5.5 

Pour approcher la solution d'un systeme differentiel avec condition initiale, on 
utilise la methode numerique associee a la fonction <I> suivante : 

${x,y,h)= ^[fci+4fc 2 + fc 3 ] 

h = f(x,y) 

r, h ft, . 

h = f{x+ -,y+ -ki) 

k 3 = f(x + h,y + h{-k x + 2fc 2 )) 

Montrer que cette methode est au moins d'ordre 3. 
Reponse 5.5 

Pour (x, y) fixe, on pose : 

h(h) =f(x,y) 

k 2 (h) =f(x+%,y+$k 1 {h)) 

k 3 (h) = f(x + h,y + h(-k 1 (h) + 2k 2 (h))) 
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Pour simplificr les notations, on notera par / a la place de /(., .), par f x a la 
place de ^/(., ■), par /„ a la place de ^/(., .), par f xy a la place de gfg^/(., ■), 
par / xx a la place de §^f{., ■) et par f yy a la place de ^ /(., .). 



On a alors 



y"(z) =h + f.f v 

y"'(x) = f xx + 2f.f xy + f x .f y + f.{f y f + ; 2 ./ w 



et 



D'ou 



k[(h) =0 

= lifxx + 2h(h)f xy + kl(h)f yy 

k' 3 (h) =f x + (-fci(/i) + 2k 2 (h) + 2hk' 2 (h))f y 
h" (h) = f xx + 2{-h{h) + 2k 2 {h) + 2hk' 2 {h))f xy + (-fci(ft) 
+2fc 2 (ft) + 2hk' 2 {h))f yy + (4fc 2 (/i) + 2hk 2 "(h))f y 

k'M = \ (fx + /■/„) 

fc 2 "(0) =l(fxx + 2f.fxy + f 2 .fyy) 
k' 3 (0) =fx + f-fy 

fc 3 " (0) = / ra + 2/./ x „ + + 2./, (/ a + /./„) 

Un developpement limite a l'ordre 2 par rapport a ft, au voisinage de zero, (x, y) 
fixe, de $(x,y,h) nous donne : 

*(x,y,h) = + /•/„)+ 

f [L + 2/./ iy + / 2 / X!/ + (/* + /./ w )/ w ] + o(h 3 ) 

et si x € [a, 6] et y dans un compact, alors, le reste est une quantite 0(h 3 ) 
uniformement par rapport a (x,y). 

De meme, un developpement limite a l'ordre 3 par rapport a h au voisinage de 
zero, x fixe de y(x + h) nous donne : 

y(x + h)- y(x) = y'(x)h + ^y»(x) + ^y»'{x) + 0(ft 4 ) 

et si x £ [a, b] alors, le reste est une quantite 0(h A ) uniformement par rapport 
a x. 

Si on remplace les expressions de y'(x), y"(x), y"'(x), on obtient que : 

y{x + h l- y{x) -Hx,y{*),h) = 0(h>) 

avec le reste est une quantite 0(h 3 ) uniformement par rapport h x E [a,b] car 
y(x) e y([a, b]) — un compact. D'ou, pour toute subdivision (x m )o< m <N de pas 
uniforme h de l'intervalle [a, b], on a : 

max 1 1 ^ m+l) ,^ (Xm) -*(x m , y(x m ), h) \ \ = 0(h 3 ) 

m fl 

cc qui prouvc, d'apres un theoreme du cours, que la methode est au moins 
d'ordre 3. 
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Exercice 5.6 

Pour approcher la solution d'un systeme differential avee condition initiate, on 
utilise la methode numerique associee a la fonction $ suivantc : 

®{x,y,h) = f{x,y) + ^g(x + ^,y+ ^f{x,y)) 

avec 

d d 
9{x,y) = T^f(x,y) + f(x,y) — f(x,y) 

Montrer que cette methode est au moins d'ordre 3. 
Reponse 5.6 

Mcme raisonnement que l'exercice precedent. 
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Chapitre 6 

Programmation Lineaire 



6.1 Introduction 

Un problemc dc programmation lineaire (on dit aussi un programme lineaire) 
consiste a minimiscr ou a maximiscr une fonction lineaire sous contraintes 
lincaircs : 



Min c T x = Yh=i °i x i 
(P ) t ajx = dj, j = 1, • • • ,p (contraintes-egalites) 

[ ajx < dj , j = p + 1 , • • • , m (contraintes-incgalitcs) 

ou x, c, ai , a,2, ■ ■ ■ ,a m sont des vecteurs de M. n et d\ , cfe , • • • ,d m sont des nombrcs 
reels. 

On peut toujours exprimer un programme lineaire sous la forme suivante : 



(P) 



Min z(x) = c T x = YH=i °i x i 
Ax = d 

x = (xi x 2 ■ ■■ x n ) T > 



oii A est une matrice m lignes (to contraintes-egalites), n colonnes (n variables) 
ct d G l m . En effet, il est toujours possible de se ramener a ce cas : 

- une contrainte-inegalite se transformc en une contrainte-egalite en ajou- 
tant une variable dite d'ecart : 



ijx < dj 



l,j — > ajx + Vj = dj, Vj > 

une variable Xi de signe quelconque est remplacee par x\ — x~l , avec 
xf,xr>0 

On dit que (P) est la forme standard de (Po)- 
EXEMPLE 

Min z = hx\ — Zx-2. 



(Po) 



X\ — x 2 > 2 
2xi + 3x 2 < 4 
—x\ + 6x2 = 10 
xi > 0, x 2 e R 
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En introduisant les variables d'ecart, on obtient la forme standard de (P ) 



(P) 



Min z = 5xi — 3(x 2 — x 3 ) + 0x 4 + 0x 5 

x\ - (x 2 - x 3 ) - X4 = 2 

2xi + 3(x 2 — xs) + x 5 = 4 

-xi + 6(x 2 - x 3 ) = 10 

xi > 0, x 2 > 0, x 3 > 0, x 4 > 0, x 5 > 



(5 - 3 3 0),d T = (2 4 10) 

1-11 -10 
2 3 -3 1 
-16 -6 



REMARQUE 6.1.1 On peut toujours supposer que rang(A) = m. En ef- 
fet, si rang(A) < m, une ou plusieurs lignes de A peuvent s'exprimer comme 
combinaison lineaire des autres lignes, alors ou bien V ensemble des solutions 
de Ax = d est vide, ou bien les equations qui correspondent a ces lignes sont 
redondantes et peuvent etre eliminees. On obtient alors une nouvelle matrice 
telle que son rang est egal au nombre de contraintes. 



6.2 Formulation d'un programme lineaire 

6.2.1 Exemple : le probleme du portefeuille 

Supposons que c'est le debut de Fan 2000, et que vous etes un investisseur 
avec un portefeuille compose de trois types d'actions : 75 actions avec la STB, 
1000 actions avec la BNA et 25 actions avec la BIAT. Les prix actuels sont les 
suivants : 20D par action de la STB, 2D par action de la BNA et 100D par 
action de la BIAT. 

On suppose que l'etude du marche et de certains facteurs economiques vous 
pcrmcttent d'estimer les benefices de chaque action durant l'an 2000 et le prix 
de chaque action a la fin de l'an 2000. Par exemple : le benefice de la STB est 
5D par action, pas de benefice pour la BNA, le benefice de la BIAT est de 2D 
par action, et a la fin de Fan 2000, la valeur de la STB est 18D par action, de 
la BNA est de 3D par action et de la BIAT est 102D par action. 

Vous n'avez pas des frais a payer pour acheter ou vendre des actions et vous 
pouvez vendre ou acheter une fraction d'action. 

Vous devez ajuster votre portefeuille pour maximiscr vos benefices de Fan 
2000. Pour ce faire, vous devez satisfaire les restrictions suivantes : 

1. Apres ajustement, la valeur de votre portefeuille rcstc la mcmc ; e'est- 
a-dire, investir la valeur actuelle de votre portefeuille. 

2. II faut tenir compte de {'inflation, la valeur de votre portefeuille a la 
fin de Fan 2000 doit etre augmenter au moins de 5% de sa valeur actuelle. 

3. Pour eviter les risques du marche, il faut que votre portefeuille soit 
cquilibre : apres ajustement, la valeur de chaque type d'action doit etre supcricurc 
ou egale a 25% de la valeur totalc de votre portefeuille. 

4. L'ajustement des actions se fait une seule fois au debut de l'annee et 
le nombre d'actions par type d'action doit etre positive. 
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6.2.2 Programme lineaire du probleme du portefeuille 

On note par x\ la variation des actions de type STB, x-i la variation des 
actions de type BNA et par x 3 la variation des actions de type BIAT ; c'est-a- 
dire, si X\ = 50, votre investissement avec la STB devient 75 + 50 = 125 actions, 
si X2 = —30, votre investissement avec la BNA devient 1000 — 30 = 970 actions. 

Vous dcvcz maximiscr lc benefice de votre portefeuille, c'est-a-dire, maximi- 
ser la fonction affine : 

5(75 + Xl ) + 0(1000 + x 2 ) + 2(25 + x 3 ) 

ou encore 

5xi + 2x 3 + 425 
ce qui revient a minimiser la fonction lineaire 

min — 5xi — 2x 3 

Xl,X2,X3 

La premiere contrainte est donnee par la restriction 1 : la valeur actuelle est de 

20 x 75 + 2 x 1000 + 100 x 25 = 6000 
le portefeuille ajuste doit avoir la meme valeur 

20(75 + xi) + 2(1000 + x 2 ) + 100(25 + x 3 ) = 6000 

ou encore 

20xi + 2x 2 + 100x 3 = 

La seconde contrainte est donnee par la restriction 2 : la valeur du portefeuille 
a la fin de l'an 2000 doit etre superieure ou egale a 6000 + 5 x 6000/100 = 6300 

18(75 + xi) + 3(1000 + x 2 ) + 102(25 + x 3 ) > 6300 

ou encore 

18xi + 3x 2 + 102x 3 > -600 

La troisieme contrainte est donnee par la restriction 3 : la valeur de chaque type 
d'action est superieure ou egale a 1500 

20(75 + xi) > 1500 
2(1000 + x 2 ) > 1500 
100(25 + x 3 ) > 1500 

ou encore 

xi > 
x 2 > -250 
x 3 > -10 

La quatricmc restriction est contcnue dans la troisieme. 

La formulation du probleme du portefeuille est donnee par le programme 
lineaire suivant : 

Min z — — 5xi — 2x 3 
20xi + 2x 2 + 100x 3 = 
18xi + 3x 2 + 102x 3 > -600 
xi > 
x 2 > -250 
x 3 > -10 
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6.3 Definitions et proprietes 



Un ensemble convexe de la forme 

U = {x e R n /Ax = d, x>0} 

avec A une matrice to x n, d S K m et rang(A) = to, est appele un polyedre. Les 
elements de U sont appeles les solutions realisables du programme lineaire (P). 
On appelle sommet du polyedre U tout point de U qui ne s'ecrit pas comme une 
combinaison convexe d'autres points de U. On appelle base toute sous matrice 
B reguliere (to x to) extraite de A ( il y a au moins une puisque rang(A) = to). 
On note 1 < 61, b 2 , ■ ■ • ,b m < n les indices des colonnes, appelees colonnes de 
base, qui forment la sous matrice B, 1 < hi, h 2 , • • • , h n — m < n les indices des 
autres colonnes de A 7 appelees colonnes hors-base et H la sous matrice formee 
par les colonnes hors-base, c'est a dire, si on note par Aj, j = l,-- - , n, les 
colonnes de A, on a : 

la fcieme colonne de B — Ab k , 1 < k < m 
la fcieme colonne de H — Ah k , 1 < fc < n — to 
{1, 2, • • • ,n} = {bi,b 2 , • • • , b m , hi, h 2 , ■ ■ ■ , h n - m } 

Soit x — (xi x 2 ■ ■ ■ x n ) T E R n un vecteur quelconque et B une base de A, la 
notation x — [xb,x h] veut dire : 

xb = fabi x b 2 ■■■ x bm ) T , x H = (x hl x h2 ■■■ x hn _ m ) T 

d'ou : 

n m n—m 

Ax = ^2xjAj = ^2x bk A bk + ^2 x h k A hk 

■j = l k=l k=l 

Le systeme Ax = d est done equivalent a 

Bx B + Hx H = d (6.1) 
On appelle solution de base (associee a la base B), la solution particuliere de 



(6.1) obtenue en faisant xu — 0. Le vecteur xb est alors determine d'une fagon 
unique par : 

x B = S _1 d. 

Une base B est dite realisable si la solution de base est realisable e'est-a-dire 
xb > 0. Une base realisable est dite degeneree si l'une des composantes de 
xb = B~ x d est nulle. Une solution realisable (resp. une base realisable) est dite 
optimale si elle (resp. la solution de base associee) realise le minimum de (P). 
Pour x = (xi X2 ■ •■ x n ) T £ M. n , on note 

I*{x) = {j / Xj > 0}. 



THEOREME 6.3.1 L'ensemble des sommets du polyedre 

U = {xe W 1 /Ax = d, > 0} 

correspond d l'ensemble des solutions de base des bases realisables. En particu- 
lier, l'ensemble des sommets est fini. 
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LEMME 6.3.1 Soit u un sommet de U. Alors, la famille {Aj/j e I*(u)} est 
libre. 



Demonstration ( Lemme ) 

Supposons le contraire, c'est a dire : 

3 ("j)j€/*(«)/ a i A J =0 
3 ei'{u) 

avec aii moins un aj ^ 0. On pose 

aj = 0,Vj£/*(u) 

a T = (ai a 2 • • • a„) 
/i = {j/aj > 0}, I 2 = 07«i < °} 
II est facile de voir que 

A(u - 9a) = d, ye e R. 

D' autre part, si I\ ^ 0, alors : 

Vj e 7i, Uj - 0aj > Vj e h, 8a j < Uj 

Vj E Ii,6 < Uj/aj 
■^==> 9 < 0i = mirije/j (itj/ctj) 

sinon, on obticnt que w — 9a > pour tout 9 > 0, on pose dans ce cas 0i = +oo. 
De me me, si I2 7^ 0, alors : 

Vj e h, Uj - 6a j > Vj e h, 9a 3 < Uj 

Vj E I2, 9 > Uj/aj 
9 >9 2 = ma,Xj e i 2 (uj / aj) 

sinon, on obticnt que u — 9a > pour tout 9 < 0, on pose dans ce cas 02 = —00. 
II est clair que 9 2 < et 0i > et que 

V0 e]0 2 ,0i[,u-0a€ ^ 

D'autre part, soit G]02,0i[— {0} tel que —9 €]02,0i[ (un tel 9 existe car 2 < 
< 0i ), on a : 

u±9aeU 7 et u = ^(u - 9a) + ^(u + 9a) 

par consequent u s'ecrit comme combinaison convexe d'elements de U, ce qui 
contrcdit lc fait que u est un sommet. 

Demonstration ( Theoreme ) 

Soit u un sommet de U, d'apres le lemme precedent la famille {Aj / j e 
I*(u)} est libre, done, on peut completer cette famille par d'autres colonnes de 
A pour obtcnir une base B. II est clair que 

Bub = UiAi = d 

iei'(u) 
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d'ou u est la solution de base de la base B. 

Invcrsement, soit B une base realisable et soit u = [ub,uh] = [B~ 1 d, 0] la 
solution de base de la base B. Supposons que u n'est pas un sommet de U : 

3X e]0, 1[, 3v ^ w G £//u = Xv + (l- X)w 

d'ou 

n, rln1 M / c- TT I J U B = AU B + (1 - X)W B 

3X G 0, 1 [,3v w G (7/ < . . .< 

t w ff = Aw H + (1 - X)w H 

Or, u# = 0, ce qui prouve que vh — wh — 0, car vh > ct wh > 0. D'autre 
part : 

u G U =>• -Bub + Huh = d 
veU Bv B + Hv H = d 

w £ U Bwb + Hwh = d 

Par consequent 

u b — v b — w b — B~ d 

d'ou u = v = w. Done u est un sommet de U. Enfin, la matrice A est constituce 
de n colonnes, et une base est constitute de m colonnes, done, le nombre de 
bases realisables est au plus egal a C™. 

COROLLAIRE 6.3.1 Soit u G U . Alors, u est un sommet de U si et seule- 
ment si la famille {Aj; j € I*(u)} est une famille libre. 

DEMONSTRATION : 



Soit u £ U un sommet de U, d'apres le theoreme 6.3.1 il existe m- colonnes 
de A : A^ , Ai, 2 , • ■ • , Ab m , qui forment une base realisable B telle que : 



= [ub 0] avec ub — {ub 1 u b 2 ''' u b m ) = B d 



d'ou, I*(u) C {&i,&2j"' T°m}- Comme B est inversible, done, ses colonnes 
Af >1 , Ab 2 , • • • , A\, m sont lineairement independantes et par suite la famille {Aj ; j £ 
I*(u)} est libre ( car la sous- famille d'une famille libre est une famille libre ). 

Soit u G U tel que {Aj ;j £ I*(u)} est une famille libre, done, on peut 
completer cette famille par d'autres colonnes de A pour obtenir une base B = 
{A bl , A b2 , - ■ ■ , A bm } avec I*(u) C b 2 , • • • , b m }. D'autre part : 

TCI 

d = Au = u jAj = UbiAbi 

j£l*{u) 1=1 

ce qui prouve que u — [ub 0] est la solution de base associee a la base B formee 
par les colonnes {A bl , A b2 , • • • , A bm }. 

THEOREME 6.3.2 Un polyedre non vide 

U = {x£ R n /Ax = d,x>0} 

possede au moins un sommet. 
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DEMONSTRATION : Soit u € U. Supposons que la famille {Aj\ j £ I*(u)} 
est liee, done 

3(« 3 k;'(»)/ E a J A j = ° ( 6 - 2 ) 
3 ei'(u) 



avec au moins un aj ^ 0, j £ I*( u )- Quitte a multiplier (6.2 1 par -1, on peut 
supposer que I = {j S I*(u)/aj > 0} ^ 0. Soit 

6> = mm — = 

j&l aj a jo 

on considere le vecteur w £ E" : 

uj — 6a j, j £ I*(u) 
0, sinon 

Montrons que w S £/ ; en effet, on a : 



Aw = — 53_ygZ*(u) a 3 A j = ^ w = ^ 

w j = u j — — ^ ^ 

Wj = Uj — 6*<x, > Uj > 0, Vj € I*(u)\I 

Wj = 0,Vj <£l*(u) 



D' autre part, 



^io = u jo ~ 

r(w) c r(u) 



Par consequent 



I*(w)Cl*(u)\{j } 



Si la famille {Aj \ j £ J*(u>)} est liee, on recommence le memo raisonnemcnt avec 
w au lieu de u. A chaque etape le nombrc d'elemcnts dc la famille lice decroit au 
moins de 1 ; alors, on obtient necessairement apres un nombre fini d'etapes un 
vecteur w € U tel que la famille {Aj]j £ I*(w)} soit libre ou I*(w) = 0, dans 
ce dernier cas w — 0. II suffit done de prouver que si £ U, est un sommet de 
U. Supposons qu'il existe < A < 1 et vi ^ v 2 £ U tel que : 

Q = Xvi + (1 — \)v 2 

comme v\ > et v 2 > 0, necessairement vi = i>2 = 0. 

PROPOSITION 6.3.1 On suppose que U est non vide et que la fonction 
lineaire z(x) = c T x est minoree sur U . Soit v £ U . Si v n'est pas un sommet 
de U , alors, il existe v £ U tel que : 

z{v) < z(v) 

et 

card{I* (v)) < card(I* (v)) 
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DEMONSTRATION : Si v n'est pas un sommct, il existe m ^ u 2 e U et 
< A < 1 tcls que : 

v = Ami + (1 — A)«2 

L ; ensemble des indices des composantes non nulles de y = u\ — u 2 est non vide 
car u\ ^ u 2 et il est contenu dans I*(v) car si 7^ necessairement ^ 0. 
D'autre part, Ay = d'ou : 

z(v + 8y) = z(v) + 9z(y), A(v + 6y) = Av = d,M6 e R, 

On peut supposer que 

z(y) > et I*(y) ^ 

en effet, dans le cas z(y) < et y possede des composantes negatives, on rem- 
place y par —y, e'est a dire, on prend y = u 2 — U\ au lieu de U\ — u 2 . Le cas 
z(y) < et les composantes de y sont toutes positives est impossible car : 

v + e y e u, ye > 

et comme z(y) < 0, on obtient 

lim z(v + 9y) = lim z(v) + 9z{y) = —00 

^+00 >+oo 

ce qui contredit le fait que z est minoree. II reste seulement le cas I*(y) = et 
z(y) > qui est aussi impossible car : 

v - 9y e U, V6» > 

et comme z(y) > 0, on obtient 

lim z(v — 9y) = lim z(v) — 9z(y) = —00 

6 >+oo 6 >+oo 

ce qui contredit le fait que z est minoree. 
On pose alors : 

9 = mm{-/iel*(y)}= ^ > 

Vi Via 

et 

v — v — 9y 

Montrons que v e U ; en effet, on a Av = d, il suffit done de prouver que v > : 
pour i e I*{y) ■ 

Vt = Vi - 9yi 
= yi(v_i/yi_- 9) 
> yi{0-6) 
= 

et pour i £ I* (y) 

Vi = Vi- 9y, >v t >0,(9> 0) 
D'autre part, I*(y) C I*(v), d'ou : 

I*(v) C /*(«) 
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et on a 



iq € I*(v) et Uj = 0. 
z(v) — z(v — 9y) = z(v) — 9z{y) < z(v) 



THEOREME 6.3.3 On suppose que U ^ 0. Alors on a Vune des assertions 
suivantes : 

a ) inf{z(x)\ X G f/} = — oo 
ou bien 

b) il existe un sommet w de U tel que : 

z(uu) = min z(x). 

x£U 

DEMONSTRATION : L'ensemble des sommets de U est fini et non vide, il 
existe done un sommet w de U tel que 

z(w) — min{z(w) / v est un sommet de U} 

Si z est non minoree, e'est alors l'assertion a) qui est verifiee. Supposons que z 
est minoree. Soit e > 0. Par definition de l'inf., il existe v e U tel que 

z(v) < inf {z(x) ;xE U} + e 

Si v n'est pas un sommet, on applique plusieurs fois la proposition precedente 
jusqu'a ce que Ton obtienne un v G U tel que : 

J z(v) < z(v) < inf{z(a;); x e U} + e 
1 I*(v) = ou v est un sommet 

or, si I*(v) — 0, necessairement v — et done v est un sommet de U. Par 
consequent : 

z{w) < z{v) < inf{2:(a;); x G U} + e 
Comme le sommet w est independant de s, alors z{w) — inf{z(x); x G U}. 



6.4 Resolution des programmes lineaires : Algo- 
rithme du simplexe 

LEMME 6.4.1 Soit C un convexe ferme de R k , k > 1, et soit y G R k \C . 
Alors, il existe a G M. k tel que : 

a T y < inf{a T x; x G C} 

DEMONSTRATION : On pose 

S = inf || x - y || 2 (6.3) 



L'inf. de (6.3) est atteint, en effet, C est un ferme de R k et y ^ C, alors, S > ; 
d'ou, l'inf. sur C est aussi egal a l'inf. sur {x G C / \\ x — y ||2< 2<5} qui est un 
compact, ce qui prouve que l'inf. est atteint, soit 

S =|| x - y || 2 = inf \\x-y || 2 
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On pose a = x — y. Comme C est convexe, alors, pour tout < ct < 1 ct s e C, 
le vecteur x + a (x — x) est dans C ; d'ou 

S 2 <|| x + a(x — x) — y \\\— 5 2 + 2a(x ~ y) T (x — x) + a 2 || x — x \\\ 

soit 

2{x - y) T (x -x) + a\\x-x \\j> 
par consequent, lorsque a tend vers 0+, on obtient que : 

a T (x — x) — (x — y) T (x — x) > 

d'ou 

a T x > a T x = a 1 \x — y) + a T y = 5 2 + a T y 

ce qui prouve que 

inf a T x > a T y + S 2 > a T y. 

xec 

PROPOSITION 6.4.1 On suppose que le programme lineaire (P) possede 
une solution optimale x* G U, alors : 

DEMONSTRATION : Soit z* = c T x*. On pose 

C = {{r, w) £ R x R" / r = tz* - c T x, w = td- Ax ; t £ R+ et x £ M" } 

II est facile dc verifier que C est un convexe ferme, de plus, C est un cone : 

yt > 0, V(r, w)£C=^ t(r, w) = (tr, tw) £ C 

Montrons que (1,0) ^ C ; en effet, supposons le contraire, alors : 

3(t ,x ) £ R + x R'l/w = t d- Ax = 0, 1 = t z* - c T x 

Si to > 0, alors, xo/to est une solution realisable du programme lineaire (P) et 
on a 

c T (x /t ) = z* — l/t < z* = c T x* = minc T a; 

ce qui est absurde. Si t — 0, alors, 

Ax a = 0, x > et c T x = —1 

d'ou, pour tout x solution realisable du programme (P), x + axo est aussi une 
solution realisable du programme (P) pour tout a positif, ce qui est absurde, 
car : 

c T (x + ax ) = c T x — a — ► — oo 

quand a tend vers +oo. Par consequent, C est un ferme et y = (1,0) ^ C, 
d'apres le lemme precedent, il cxistc a = (s, A) 6Rx W n tel que 

a T y — s < S = inf{a T .x = s.r + \ T w, (r, w) £ C} 
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Le nombre S est necessairement positif ; sinon, 

3a; = (r, w) G C / a T x — s.r + X T w < 

et comme C est un cone, on obtient que l'inf. est egal a — oo ce qui contredit le 
fait que S > s. D'autre part, (0,0) G C, d'oii 

5 < 

ce qui prouve que 5 = 0. Done, s < et on a : 

s.r + X T w > 0,V(r,w) G C 
Si on note A* = A/(— s), on obtient : 

-r + (A*) T u; > 0,V(r,u;) G C 
soit encore, en utilisant que (r,w) G C : 

-(tz* - c T x) + (X*) T (td - Ax) > 

d'ou 

t(d T X* - z*) + (c - A T A*) T x > 0, Vt > 0, Vx > 
Pour t = 0, on obtient que 

A T X* < c (6.4) 

Pour x = et i = 1 , on obtient : 

d T X* > z* 

Utilisons le fait que x* > et que Ax* = d, et multiplions l'equation (6.4 1 par 
x* , on obtient : 

d'ou 



d A < c x =2 



d T A* = c T x* 

THEOREME 6.4.1 Soit B une base realisable. Une condition suffisante pour 
que B soit une base realisable optimale est que : 

c H = (c hk )fe=i, n -m = c H - {clB- 1 H) T > 0, (6.5) 

avec c — [cb,ch\- Si de plus, la solution de base est non degeneree, alors, la 



condition (6.5) est necessaire. 



DEMONSTRATION : 

Soit B une base realisable. Supposons que la condition ( |6.5| est verifiee. Par 
definition, la solution de base a; = = [B~ 1 d, 0] est realisable. Soit 

x = [xg, xh] G U, on a : 

Xb > 0, Xh > 0, et Bxb + Hxh = d 
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d'ou : 

z(x) = c B x B + cJjX H 

= c T B {B- l d- B^Hxh) +cJjX H 
= c T B B-H+{c T H -c T B B^H)x H 
> c T B B~ x d 
= z(x°) 

Supposons que la base realisable B est optimale, c'est a dire que la solution 
de base, qu'on note x* , est optimale et que x* est non degeneree. D'apres la 
proposition ci-dessus, il existe A* € R m tel que : 



c — A T X* > 

t * 

C X 



d T X* = 



de la premiere equation, on deduit que : 

c bk - AZ\* > Vfc = 1, • • • ,m 
c hk - A% k \* > Vfc = 1, • • • ,n - m 

et de la deuxieme equation, en remplagant d par Ax* , en utilisant les inegalites 
donnees par la premiere equation et le fait que la solution x* est non degeneree, 
on deduit que : 

Cb k — AT A* = Vfc = 1 , ■ • • , m 



d'ou 

par consequent 
d'ou 



c B = B T X* et c H - H T X* > 
(B T )-% = A* 



ch - H' 1 '(B T )~ 1 c B > 0. 



COROLLAIRE 6.4.1 Soit B une base realisable quelconque et x° la solution 
de base correspondante. On suppose qu'il existe une composante strictement 
negative de ch d'indice hors-base h s , ch est le vecteur defini par (6.5) et on 
pose 

d = B~ x d = (dk)k=i,m, Ah, = B~ l A hs = {a i<hs )i=i,m (6.6) 

alors : 

(a) ou bien les composantes de Ah s sont toutes negatives, et dans ce cas 
I'inf. sur U est egal a — oo. 

(b) ou bien 

Mi*(4)ci*(^ 

et dans ce cas on met en evidence une autre base B realisable et une solution 
de base x ( associee a B) telle que : z(x) < z{x ) 

(c) ou bien 

3i&I*{A K )/iiI*{d) 

et dans ce cas on met en evidence une autre base B realisable avec la meme 
solution de base Xq. 
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DEMONSTRATION : L'equation (6.51 s'ecrit sous la forme suivantc 



Ch 3 = c h] - ^ ak ' h i c b k A < j < n - m 
fe=l 



et par hypothese, on a 
Soit 



Ch s < 



d = B- 1 d=(d k ) k=hm >0 
On note par (ej)j = i )n la base canonique de B" et pour 9 £ R + on pose 



'(.0) = C^2( d k - 9a kJls )e bk ) + 9e hs 



(6.7) 



k=l 



on obtient que 

Ax(9) 



Aa; + 6 (Ah, — BA hs ) 
d 



(6.8) 



D' autre part 



z(x(6)) = z(x°) +0(z(e hs ) -Y,T=i a kM„z(e bk )) 
= z(x°) + 9(c hs - J2k=i d k ,h 3 c bk )) 
= z(x°) + 9c hs 



(6.9) 



(a) Si toutes les composantes de Ah 3 sont negatives, d'apres (6.7 1 et (6.8 1 x(9) 
est une solution realisable quclquc soit lc nombre 9 > et d'apres (6.9 1 z(x(9)) 
tend vers — oo quand 9 tend vers +oo. 

(b) Dans ce cas, la valeur de 9 ne peut pas etre augmentcc indefiniment, la 
plus grande valeur de 9 est donnee par : 



iei'(A hs ) ai,h s d r ,h s 

La nouvellc solution x = x(9) £ U et verifie : 

z(x) < z(x°) 
Xb r = 
x he = 9 > 

il suffit done de prouver que B = B + {A^ s } — {A^} est une base. Supposons 
que B n'est pas une base, soit I = {1, 2, • • • , m}\{r}, alors 



3(a k )kei/A hs = ^ a k A b 



(6.10) 



kei 



D'autre part, a r ,h a > et d'apres (6.6) on a : 



Ah B = ak,h B A bk 



(6.11) 



k=l 
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d'ou, en combinant ( |6.11[ ) et ( 6.10[ ), on obtient : 

/Xak.hs - a k)A bk + a r ,h s Ab r = 
fee/ 

ce qui est impossible car (Ab k )k=i, m est une base et a~ r ,h a ^ 0. 

(c) Meme raisonnement que (b), sauf que 6 = 0, on change de base sans que 
la solution de base change. 

L'interet du corollaire ci-dessus vient du fait qu'une des methodes de resolution 
des programmes lincaires en decoule directement : Valgorithme du simplexe 

Algorithme du simplexe 

On suppose qu'on dispose d'une base realisable de depart B° 

(a) B = B° base realisable de depart. Iteration k = 0. 

(b) k <- jfe + 1 

(c) a Viteration k soit B la base courante calculer : 

d = B- l d 
H = B- l H_ 

CR = C H — H T C B = (Ch k )k=l,n-mt 

(d) Si ch > 0, STOP : I 'optimum est atteint. 
Sinon, il existe h s tel que : Ch s < alors : 

(e) Poser Af lt = B~ 1 Ah, = {ck,h s )i=x,m = lQ sieme colonne de H. 
Si 3t,/i a < Q, Vi = 1, • • ■ , m, STOP : optimum non borne (—oo ). 
Sinon, calculer : 



a r,h s iei'(A hs ) ai,H a 



(f ) Soit x tel que 



Xb k = dk — 0ak,h a , Vfc = l,m 

Xh B = 

x hk = 0, Vfc ^ s 

(xb r quitte la base et Xh s entre en base), x est la solution de base de la nouvelle 
base B = B + — {Af, r }. Retourner en (b). 

La methode du simplexe, due a G.B.Dantzig, est une procedure qui permct 
de passer, le long de la frontiere de U, d'un sommet a un autre sommet adjacent, 
ce qui permet d'atteindre le sommet optimal. Algebriquement , elle permet de 
construire une suite de bases adjacentes 

f)0 pi r>k 

et des sommet s 

-rO r i . . . ... 

avec x k , k = 0, 1, 2, • • • , une solution de base de B k tel que : 
z(x°) > z{x 1 ) > z{x 2 ) >■■■> z{x k ) > ■ ■■ 
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REMARQUE 6.4.1 Dans Valgorithme du simplexe, le choix en (d) de la 
variable Xh 3 qui va entrer en base, dans le cas oil plusieurs composantes de ch 
sont strictement negatives, est quelconque. Dans la pratique, I'experience montre 
qu'il est plus interessant de choisir la variable qui correspond a la composante 
la plus negative de Ch- 



REMARQUE 6.4.2 Si la solution de I'etape k est non degeneree le nombre 
9 est strictement positif, dans ce cas z(x k+1 ) < z(x k ). Sinon, le nombre 9 peut- 



etre nul (c'est Vassertion (c) du corollaire 6.^.1), dans ce cas la base change 
mais la valeur de z ne change pas et il est possible done de retrouver, apres un 
certain nombre de changement de base a une base deja rencontree et de cycler 
indefiniment. On appelle ce phenomene "le phenomene de cyclage". L 'experience 
montre que le phenomene de cyclage est tres rare et la plupart des programmes 
utilises dans Vindustrie ne sont pas munis d'une strategic contre le cyclage. Une 
des strategies contre le cyclage est donnee par Bland (1977) : 

- parmi toutes les variables Xh s qui peuvent entrer en base ( qui cor- 
respondent a des composantes strictement negatives de ch ) choisir celle de plus 
petit indice. 

— parmi toutes les variables xi, r qui peuvent quitter la base, choisir celle 
de plus petit indice. 

II est clair que si le phenomene de cyclage ne se produit pas, Valgorithme 
du simplexe converge en un nombre fini d'iterations, car le nombre de som- 
mets de U est fini. On montre qu'avec la strategic de Bland, la methode du 
simplexe converge au bout d'un nombre fini d'iterations, meme dans le cas de 
degenerescence. 

REMARQUE 6.4.3 La base B de I'etape k et B de I'etape k + l ne different 
que d'une colonne : la colonne r de B (— Ab r ) est remplacee par la colonne h s 
de A (= Ah 3 ). Soit C la matrice m x m definie par : 

Ci,j = Si,j , VI < i < m, VI < j < m, j ^ r 
Ci.r = -a,i t hja r ,h a , VI < i < m, i ^ r 

II est facile de voir que la matrice C est inversible car ses colonnes sont lineairement 
independantes. Si on note par (/i); = i,m la base canonique de R m , il est clair 
que : 

Cfi = /i, VI < i <m,i ^ r 
CA hs - f r 



d'ou 

ce qui prouve que 
done 



Bf t = Bfi = BC^fc, i £ r 

Bfr = A hs = BA hs — BC~ 1 f r , 

B = BC'- 1 

B~ 1 d= CB~ 1 d = Cd 
B~ 1 A = CB^A 



et de la derniere equation, on deduit facilement H puisque les nouveaux indices 
de la matrice hors-base hi, h 2 , ■ ■ ■ , h n ^ m sont connus. 
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6.5 Base de depart 



Deux questions se posent avant dc lancer la mcthodc du simplcxc pour 
resoudre un programme lineaire : la premiere est de savoir si 1 'ensemble des 
solutions realisables 

U = {x g Rl/Ax = d} 

avec A une matrice mxmde rang met de R m , est vide ou non. La deuxieme 
question, dans le cas ou U ^ 0, est alors de determiner un sommet de U, ce 
qui est equivalent a une base realisable de A, pour initialiser la methodc du 
simplexe. 

Pour rcpondre aux questions precedentes, une des methodes consiste a construire 
un autre programme lineaire tel que l'un de ses sommets se calcule facilcmcnt, 
il possede une solution optimale et en fonction de sa solution on donnera une 
rcponse. Le cas m = 1 est evident. Pour m > 2, il y a au moins une colonne de 
A, Ai , qui n'est pas un multiple de d, sinon rang(A) = 1 < to. On pose 

r = d- A lQ 

A = [A r ], matrice to x (n + 1) 
c= (0 ••• 1) T € M. n+1 
x e R n+1 /x io = x n+1 = l,Xi = 0,Vi o ^ % + n+ 1 
On definit le programme lineaire (P) par : 

/ pN f min £(5) = c T x = S„+i 

1 j \ 5 e 17 = e = d, x > 0} 

On verifie facilement que x est un sommet de U et que z est minoree par 
zero sur U. La mcthodc du simplexe permet done dc determiner un sommet 
£ = (£i £2 • • • £n £n+i) T de [7 solution optimale du programme lineaire (P). II 
est clair que £ n +i = si et seulement si U 7^ 0. De plus, si £„+i = la famillc 
{Aj, j e /*(£)} est libre; d'ou, (£1 £2 ■ ■ ■ £«) T est un sommet de [/ (pour plus 
de details, consultcr les exercices de ce chapitre). 



6.6 Etude d'un exemple 

Pour resoudre, par la mcthodc du simplexe, un programme lineaire 



(P) 



min z(x) = c T x = X)"=i c i x i 

x e U = {x e W l /Ax = d et x = ( Xl x 2 ■■■ x n ) T > 0} 



on met toutes les informations de chaque iteration dans un tableau de la forme 
suivante : 

x\ x 2 x 3 xa x n 



B~ l A 



d 



Cl C 2 C 3 C4 c„ 

h 2 b 3 hi h 3 67 

Ch 2 * c hl Ch 3 * 
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avec B la base realisable de l'iteration en cours, 61,62, • • • ,b m lcs indices des 
colonnes de A qui formcnt la base realisable B et hi, h 2 , ■ ■ ■ , /i n -m l es indices 
des colonnes hors base de A qui forment la matrice H. On rcmarque que les 
indices 6^ et les indices hk ne sont pas necessaircment ordonncs, il suffit que 
les colonnes A bl , A h2 , • • • , A bm forment unc base realisable. II est facile de 
reconnaitre les colonnes de base dans ce tableau, car, a une permutation pres, 
les colonnes de B^ 1 A qui forment l'identite de R m sont les colonnes de base 
et les autres colonnes forment, a une permutation pres, la matrice H. Ce qui 
permet de remplir la derniere ligne du tableau : 

ch = (ch k ) = c H - H T c B - 



EXEMPLE 



min z = —Xi — 2x2 
sous contraintes : 
— Zxi + 2x 2 + X3 = 2 
—Xi + 2x2 + X4 = 4 
£1 + 2:2 + x 5 = 5 

Xi,X 2 ,X 3 ,X4,X 5 > 



On a 



A = 




I -1 \ 

-2 




V / 



On remarque que les trois dernieres colonnes de A forment une base. On peut 
done prendre : 

61 = 4, 6 2 = 3, 63 = 5, hi = 2, h 2 = 1 

soit 



B 




ce qui donne 









-1 


2 





1 




A = 


( 


-3 


2 


1 











1 


1 








Xi 


X2 X 3 


X4 


x 5 






-1 


2 





1 







4 


-3 


2 


1 










2 


+1 


+ 1 








1 




5 


-1 


-2 















h 2 


hi 


62 


61 


63 






-1 


-2 


* 


* 


* 






-2 






2 


2 






)-( 




-1 


-3 


I) 



ch 



Xhx — X2 entre en base car c/^ = —2 < Ch 2 — — 1 < 
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9 = min - — 

iei*(A 2 ) ai 2 



da 

Q>2,2 



-,(r = 2) 



done Xb 2 = x 3 quittc la base. La nouvelle base est done 

h = 4, 62 = 2, 6 3 = 5, fti = 3, h 2 = 1 

Pour calculer le nouveau vecteur Jet la nouvelle matrice _ff , il suffit de multiplier 
a gauche la matrice formee par les trois premieres lignes du tableau par la 
matrice C de changement de base : 



C 




ce qui revient a : 

Li 
L 3 

ce qui donne lc tableau suivant 



Li — L 2 

1/2L 2 

L 3 - 1/2L 2 





2:2 


£3 


x 4 


x 5 




2 





-1 


1 





2 


-3/2 


1 


1/2 








1 


5/2 





-1/2 





1 


4 


-1 


-2 













^2 


&2 


hi 


bi 


b 3 




-4 


* 


1 


* 


* 





. 

''II = I _J 
1 

-4 

£/i 2 = £1 entre en base car c/j 2 



mm. ^ — 
i€/*(Ai) aj,i 



-1 1/2 -1/2 
2 -3/2 5/2 



-4 < < c hl = 1 
di di 



ai,i 



-,(r = l) 



done x 6l = x 4 quitte la base. La nouvelle base est done 

h = 1, 6 2 = 2, 63 = 5, fti = 3, ft 2 = 4 
La matrice C de changement de base est egale a : 



ce qui revient a : 



C 



Li 
L 3 



1/2 
3/4 1 
-5/4 1 



l/2Lj 

L 2 + 3/4Li 

L 3 - 5/4Li 
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ce qui donne le tableau suivant 



Xi 


X 2 


X3 


X4 


«5 




1 





-1/2 


1/2 





l 





1 


-1/4 


3/4 





5/4 








3/4 


-5/4 


1 


3/2 


-1 


-2 













h 


b 2 


hi 


h 2 


h 




* 




-1 


2 







c# = - 



-1/2 -1/4 3/4 
1/2 3/4 -5/4 



Xfix = ^3 entre en base car c hl = — 1 < 



= min - — 

lei-(Ai) a i 3 



03,3' 



,(r = 3) 



done Xb 3 = X5 quitte la base. La nouvelle base est done 

61 = 1, 6 2 = 2, 63 = 3, fti = 5, ft 2 = 4 
La matrice C de changement dc base est egale a : 



ce qui revient a : 



C = 



Li 
L 3 



1 2/3 
1 1/3 
4/3 

Li + 2/3i 3 
L 2 + I/3L3 
- 4/3L 3 



ce qui donne le tableau suivant 



X\ 


x 2 


2^3 


X4 


x 5 




1 








-1/3 


2/3 


2 





1 





1/3 


1/3 


3 








1 


-5/3 


4/3 


2 


-1 


-2 













h 


b 2 


63 


h 2 


hi 




* 


* 


* 


1/3 


4/3 





Cff = - 



-1/3 1/3 -5/3 
2/3 1/3 4/3 



1/3 
4/3 

L 'optimum est atteint car ch > 0. La solution de (P) est done la solution 
de base de la dcrnicre base : 



Xb! = d\ = 2, a; b2 = d 2 = 3, a;b 3 = d 3 = 2, a;/^ = Xh 2 
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Or, &i = 1, b 2 = 2, &3 = 3, hi = 5, /12 = 4, ce qui donne : 



/ 2 \ 

3 

2 


v y 



solution optimale de (P). 



6.7 Dualite 



Considerons le programme lineaire sous la forme standard 



T 

C X 



(P) 



min z(x) 
sous contraintes 
Ax = d 
x > 



On associe a ce programme lineaire un autre programme lineaire : 
(D) 



max w(u) — d T u 
sous contraintes 
A T u < c 



Le programme lineaire (P) est appele le primal et le programme lineaire (D) 
est appele le dual de (P). 

Exercice d'application : Determiner la forme standard de (D) et montrer 
que le dual de (D) est (P). 

Reponse : On note par Aj € K m , j = l,n, lcs colonnes de la matrice A, par 
u T = (ui u 2 ■ ■ ■ u m ) G R m et par d T = (d\ d 2 ■ ■ ■ d m ) e R m . On peut toujours 
ecrire : 

/ vi - v m+1 \ 

V2 - V m +2 



J 



avec v T — (vi v 2 ■ ■ ■ v 2m ) € M+ m , ce qui donne : 

d T u = -f T v et A T u = Bv 

avec f T = (— di —d 2 ■■■ — d m d\ d 2 ■■■ d m ) et B est une matrice n lignes 2m 
colonnes telle que sa ieme lignc = (Aj — Af). On obtient alors un programme 
equivalent au programme (D) : 



Min f T v 
Bv <c 



v e 



p 2m 



Enfin, on obtient la forme standard de (D) en ajoutant n-variables positives a 
v, ce qui donne : 

Min f T ii 
(D) { Bv = c 
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avec f T = (f T H «), v T = (v T v 2m +i ■ ■ ■ v 2m +n) et B = [B I n \. Par consequent 
le dual de (D) est : 

{Max c T y 
B T y < f 
y e R n 

Si on remplace B et f par leurs valeurs en fonction de A et d et y par —x, on 
obtient : 

Max c T (—x) 
A{-x) < -d 
(P) { -A(-x) < d 
(-x) < 
x G W 1 

ce qui donne le programme (P). 

Les relations entre le probleme primal (P) et son dual (D) sont donnees par 
le theoreme suivant : 

THEOREME 6.7.1 (a) Six etu sont respectivement des solutions realisables 
du primal et du dual, alors : 

z(x) — c T x > w(u) = d T u 

(b) Si I'ensemble des solutions realisables du primal et I'ensem-ble des 
solutions realisables du dual sont non vides, alors, le programme primal admet 
une solution optimale. 

(c) Si x* et u* sont respectivement des solutions realisables du primal et 
du dual telles que : 

T * jT * 

c x = a u 

alors, x* est une solution optimale du primal et u* est une solution optimale du 
dual. 

(d) Si (P) admet une solution optimale x* , necessairement, (D) admet 
une solution optimale u* et on a : 

c T x * = d T u* 

DEMONSTRATION : 

(a) 

Ax = d =>■ u T Ax = u T d = d T u = (Ax) T u = x T (A T u) 

et comme x > 0, A T u < c, alors, d T u < x T c = c T x. 

(b) Si I'ensemble des solutions realisables du primal et I'ensemble des solu- 
tions realisables du dual sont non vides, d'apres (a), la fonction z est minoree 
et le theoreme 6.3.3 montre que (P) admet une solution. 

(c) On utilise (a) et (b). 

(d) Voir proposition 6.4. 1| 

COROLLAIRE 6.7.1 Etant donne un programme lineaire (P) et son dual 
(D): 

(a) Si (P) et (D) ont des solutions realisables, alors, chacun d'eux a une 
solution optimale et : 

z* = min(P) = max(P) = w* 
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(b) Si I'un des programmes (P) ou (D) admet une solution optimale, 
alors, I 'autre admet aussi une solution optimale et : 



z* = min(P) = max(P) = w* 



(c) Si I'un d'eux a un optimum non borne I'autre n'a pas de solution 



DEMONSTRATION : 

(a) D'apres (b) du theoreme precedent, lc programme (P) admet une solu- 
tion optimale x* et d'apres (d) du meme theoreme le programme lineaire (D) 
admet une solution optimale u* telle que : c T x* = d T u* . 

(b) Si (P) admet une solution optimale, on utilise (d) du theoreme precedent, 
on obticnt cc qu'il faut. Si (D) admet une solution optimale alors sa forme 
standard admet une solution optimale et d'apres l'excrcicc d' application, le dual 
de la forme standard de (D) n'est autre que (P), done encore une fois d'apres 
(d) du theoreme pecedent, on obtient que (P) admet une solution optimale. 

(c) Si l'optimum de (P) est non borne, e'est a dire : inf(P) = — oo, en 
utilisant (a) du theoreme precedent, on vcrific facilcmcnt que 1 'ensemble des 
solutions realisables de (D) est vide. Meme raisonnement si l'optimum de (D) 
est non borne, e'est a dire : sup(D) = +00. 

6.8 Exercices du chapire 6 : 

Exercice 1 

Soit 



1. Resoudre graphiquement le programme lineaire (Vo) ■ 

2. Donner un equivalent du programme lineaire (Vo) sous la forme standard. On 
notera par (V) le nouveau programme. 

3. Resoudre par l'algorithme du simplexe le programme lineaire (V). Conclure. 



1. Determiner le nombre 7 pour lequel l'ensemble des x £ R 4 qui verifient les 
contraintes soit non vide. 



realisable. 



(PoW 



min — xi — 2x2 
—Xi + x 2 < 2 
Xl + X2 < 4 
Xi > 0, X2 > 



Exercice 2 

Soit 




min 2xi — 3x2 + X3 + 3x± 
2x\ — X2 — X3 — X4 — 27 

— Xl +X2+ Xz + 2X4 = —7 

3a;i — X2 — xz — 27 + 1 

Xl > 0,X2 > 0,X3 > 0,X4 > 
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2. On notera par (V) le programme lineaire pour la valeur de 7 obtenue en 1.). 
Montrer que le programme (V) s'ecrit sous la forme : 



, . J min 2xi - 3x 2 + x 3 + 3 
1 j \ a; € U = {x € K 4 /Ar = 



3:r 4 
d;a; > 0} 



ou A et rf sont des donnees a determiner avec rang(A) — 2. 
3. Resoudre par l'algorithme du simplexe le programme lineaire (V). 

Exercice 3 

On considere le programme lineaire : 

fmin 2xi — X2 + x A 
xi + x 2 - x-s > -2 
-XI + X 2 + 2X 3 < 1 
xi + x :i = 1 
x\ > 0, x^ > 0, X3 > 

1. Ecrire le programme (Vo) sous la forme standard (V) : 

f min = c T a; 
1 J \ xeU = {x£ W 1 /Ax = d et Xi>0,i = !,-■■ ,n} 

ou n, c, d, ^4 sont des donnees a determiner. 

2. Verifier que x = (| | §0) T est un sommet de W. 

3. Resoudre par la methode du simplexe le programme lineaire (V) et deduire un 
vecteur de R 3 qui realise le minimum de (Vo)- 

4. Transformer le programme (Vo) en un programme lineaire de dimension 2. 

5. Resoudre graphiquement le programme lineaire de dimension 2. Comparer la 
solution obtenue par la methode graphique et la solution obtenue par la methode 
du simplexe. 

Exercice 4 

Une usine fabrique deux types de jouets en bois : des soldats et des trains. Les 
donnees de ce probleme sont representees dans le tableau suivant : 





P. vente 


Mat. prem. 


Frais gen. 


Menuiserie 


Finition 


1 soldat 


27DT 


10DT 


14DT 


fh de travail 


2h de travail 


1 train 


21DT 


9DT 


10DT 


fh de travail 


fh de travail 



Par semaine l'usine dispose de toutes les matieres premieres necessaires a la fa- 
brication et ne dispose que de fOOh de finition et 80h de menuiserie. La demande des 
trains et des soldats est illimitee. Determiner le plan de production qui maximise le 
profit de l'usine. 

Exercice 5 

Le self-service d'un hotel offre chaque jour a ses clients quatre plats : platf , plat2, 
plat3, plat4. Le prix d'une unite du platf vaut 0.5DT, du plat2 vaut 0.2DT, du plat3 
vaut 0.3DT et du plat4 vaut 0.8DT. Le tableau suivant nous donne la quantite de 
vitamines Vf , V2, V3 et V4 dans une unite de chaque plat : 



par unite 


Vf 


V2 


V3 


V4 


platf 


400 


3 


2 


2 


plat2 


200 


2 


2 


4 


plat3 


150 





4 


1 


plat4 


500 





4 


5 
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Un client suit un regime alimentaire doit manger au moins : 500 unites de VI, 6 
unites de V2, 10 unites de V3 et 8 unites de V4. Determiner le regime qui coute le 
moins cher. 

Exercice 6 

Soit (V) le programme lineaire : 

,~y. j min z{x) = c T x 

1 ' \ x<£U = {xe R n /Ax = d;x>0} 

avec A une matrice m x n de rang(A) = m, c T = (cic 2 • • • c„) £ 1" ct rf £ K m . On 
notera par A,- G E m , j = 1, • • • , n , les colonnes de la matrice A et on suppose que 
m > 2. 

1. Montrer qu'il existe jo € {1, • • • , n} tel que Aj n'est pas un multiple du vecteur 
d. 

2. On pose A n+1 = d — A, , A,- = Aj pour j — 1, - ■ ■ , n, ^4 la matrice m x (n + 1) 
dont les colonnes sont Aj, j = 1, • • • , n + 1 et c 2 " = (0---0 1) G K n+1 . On 
considere le nouveau programme lineaire : 

CP) / min = = Xn+l 

( '{ x G U = {x G E n+1 /^ = d; x > 0} 

(a) Montrer que x G R n+1 avec Xj = x n +i = 1 et x, = pour i / jo et 
i 7^ n + 1, est un sommet de U. 

(b) Montrer que le minimum de (V) est atteint en un sommet de U. On notera 
par £ T = (£i ^2 - ■ ■ in Cn+i) un sommet de hi qui realise le minimum de 

CP). 

(c) Montrer que : 

£n+i = W = {x G K n /^x = d; x > 0} / 

(d) Montrer que si £ n +i = alors : 

x T — (^i • • • £„) est un sommet de U = {x G E n / J 4x = d; a; > 0}. 



6.9 Corrige des exercices 



Reponse 1 

Soit 



CPo) 



min — xi — 2x2 

—X\ + X2 < 2 
Xl + X2 < 4 
Xl > 0, X2 > 



1. II est clair que la fonction lineaire z(x) — —xi — 2x2 est continue et que 
l'ensemblc U des points qui verifient les contraintes est un compact de R 2 , 
done le minimum est atteint. D'apres un theoreme du cours, le minimum est 
atteint en un sommet de W, or les seuls sommets de Li sont (0 2), (1 3), 
(0 0) et (4 0) et on a respectivement : -4, -7, et -4 comme images par z(x). 
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D'ou, le sommet (1 3) realise le minimum de (Vo) et sa valeur est egale a -7. 




-Xl + X2 = 2 



X'l + x 2 = 4 
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min — x\ — 2x2 
I -xi + x 2 + x 3 = 2 

' ' Xl + X2 + Xi = 4 

xi > 0,^2 > 0,2:3 > 0,2:4 > 
3. A = ^ ~* I q ° j , c T = (-1 - 2 0) et d = ( ^ j . On voit que 

(l 1 ) 6St Une base ^ = 2 > & 2 = 4 > hi = 1, h 2 = 3), = ^ 1 ^ 
et la solution de base associee est x = [xb 0] avec is = B~ 1 d = 
x T = (0 2 2) et «(a;) = -4. On calcule : i? = B^H, on trouvc 




1 



d'ou 



H 



-1 1 
2 -1 



= c H - 


-H t cb 




v 


)- 


(' 






x 2 


X.i 




d 


A 


-1 


1 


1 





2 




1 


1 





1 


4 




hi 


61 


fta 


62 




c 


-1 


-2 












-3 


* 


2 


* 





-1 2 
1 -1 



-2 




hi=l entre dans la base et Ai = B' 1 Ai = (-1 2) T , d'ou, J*(Ai) = {2} 

di d2 



mm 

ie{2} aii a 2 i 



done, 62 = 4 quitte la base. 
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4. La nouvelle base est done avec 61 = 2, 62 = 1, hi = 4, hi = 3 B = ( j ^ J 

et B' 1 = i ^ ^ I j d'ou x T = (1 3 0) et z(z) = -7. On calcule : 
_ff = B _1 H, on trouve 



1 1 
1 -1 



= c H -H T 


CB = 


(0 


H 


(I 




Xl 


£2 


£3 


2:4 


d 


A 


-1 


1 


1 





2 




1 


1 





1 


4 




hi 


61 


h 2 


&2 




c 


-1 


-2 










ch 


* 


* 


3/2 


1/2 





1 



-2 
-1 



3/2 
1/2 



ch > 0, alors le minimum est atteint. On garde seulement les variables de notre 
probleme initial et on laisse tomber les variables artificielles xz et X4. 
La conclusion : on trouve la meme solution avec les deux methodes, done, il vaut 
mieux, dans le cas de la dimension 2, utiliser la methode graphique. 

Reponse 2 

1. La troisieme equation est egale a 2 fois la premiere + la deuxieme. Pour que le 
systeme admet une solution, il faut que 2 fois le second membre de la premiere 
equation + le second membre de la deuxieme est egal au second membre de la 
troisieme equation, ce qui donne : 7 = 1. 

2. Pour 7 = 1, la troisieme equation est inutile. Le programme lineaire est done : 



CP) 



Ax = d 
x > 



avec A — 



-1 



-1 



1 



d 



2 

-1 



et c T = (2 - 3 1 3). II est 



clair que le rang de A est egal a 2 et on voit que B = 

(10 0) et z(x 



base realisable de solution de base x T 
b'2 = 2, hi = 3, h 2 — 4. Si on calcule : 



2 -1 
-1 1 
= 2, done, 61 



est une 
1, 



ch 



c h - H c B 



4 
10 



> 



done x T = (1 0) realise le minimum de (V) et la valeur du minimum est 2. 

Reponse 3 

1. Les variables xi,X2,X3 sont positives d'ou : 

Xl + X2 — Xz > — 2 <^=> 3x'4 > 0/iEl + X2 — Xz — X4 — — 2 
—351 + X2 + 2xz < 1 <=> 3X5 > 0/ — Xl + X2 + 2xz + X5 = 1 



D'ou : 



(To) 



min 2xi — x 2 + x 3 

Xl + X 2 — Xz — X4 = —2 
— Xl + X2 + 2X3 + %5 = 1 

Xl + x 3 = 1 
Xl > 0,X2 > 0, X3 > 0, X4 > 0, £5 > 
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ce qui est encore equivalent a : 



x e U = {x e R 5 /Ax = d,Xi>0,i = l,--- ,5} 



avec c T = (2 - 1 1 0), d T = (-2 1 1) et A 



1 

-1 
1 



II faut d'abord prouver que i£W 
x > et 



Ax = 



( 1/3 \ 


2/3 
5/3 

V o / 



= d 



Montrons maintenant que x est un sommet : d'apres un theoreme du cours, il 
faut et il suffit que {Ai/xi > 0} est une famille libre (Ai —la, ieme colonne de 
A). Dans notre cas la famille en question est {Ai, A3, A4}. Soit (ct,/3, 7) G E 3 , 
on a : 



olA x + P A 3 + -(A 4 = 




L'equation (3) + l'equation (2) nous donne 3/3 — 0, d'ou /3 = 0. De l'equation 
(3), on deduit que a = et de l'equation (1) on deduit que 7 = 0. D'ou la famille 
{Ai, A3, A4} est une famille libre. Done, x est un sommet de U. 

La matrice de base realisable associee au sommet x est : 



B = 



et la matrice hors-base est : 



H 





avec b\ = 1, 62 = 3, 63 = 4 



avec hi — 2, h2 — 5 



Une iteration de la methode du simplexe necessite le calcul de H = _B _1 _H" (ce 
qui est equivalent au calcul de A2 = B^ 1 A2 = (01,2 02,2 ffl3,2) T et A§ = 
B~ 1 Ar > = (01,5 0,2,5 a3j > ) T ) et le calcul de d — B~ 1 d . Pour ce calcul, on 
peut adapter la methode de Gauss- Jordan, a gauche la matrice B et a droite la 
matrice formee par les colonnes A 2 , A 5 et d : 
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-1/3 1/3 
1/3 2/3 



-2 5 



h + 2/3l 3 
h + 1/3J 3 



D'ou, les deux premieres colonnes de la matrice obtenue a droite represente H 
et la troisieme colonne est egale a d. On doit maintenant calculer le vecteur : 

ch — ch — U cb 



avec ch = ( °^ ) = ( n ] et cb = \ c 3 | = [ 1 | . D'ou : 



ch 




) 3 ( -1 1 -2 ) M ( 1/3 ) ( £i 



Ch 1 = — 2/3 < alors x ne realise pas le minimum de (V) et hi = 2 entre en 
base. D'autre part I*{A hl ) = r(A 2 ) = {2}, d'ou : 

n . dk d r d2 

t) = mm - — = - — = - — 

k£I(A 2 ) dk,2 a r ,2 0,2,2 

ce qui donne r = 2 et b r — 62 = 3 quitte la base. La nouvelle base realisable 
obtenue a la suite de la premiere iteration de la methode du simplexe est alors : 
bi = 1, 62 = 2, 63 =4,hi = 3, /i 2 = 5. 

La matrice de base realisable associee au nouveau sommet est : 



B= -1 1 avec 61 = 1, 62 = 2, 63 = 4 




et la matrice hors-base est : 




H — \ 2 1 avec hi — 3, h>2 — 5 



Le calcul de = B^ 1 !! (est equivalent au calcul de ^4 3 = B^ 1 A3 — (01,3 02,3 a3,3) T 
et ^5 = B^ 1 A$ = (ai,5 02,5 fl3,s) T ) et le calcul de d = B~ 1 d . Pour ce calcul, 
on peut adapter la methode de Gauss- Jordan, a gauche la matrice Bet a droite 
la matrice formee par les colonnes A2, A 5 et d : : 
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-3/2 -1/2 -3/2 

1 1 -1 

5/2 1/2 5/2 




D'ou, les deux premieres colonnes de la mat-rice obtenue a droite represente H 
et la troisieme colonne est egale a d. On doit maintenant calculer le vecteur : 



ch — ch — tl Cb 



avec ch = = n et eg = c 2 = —1 • D'ou : 



> 



Done, le sommet x associe a cette base realise le minimum de (V). Le vecteur 
x = [xb Xh] avec 




xb = 





= B 1 d — d— 2 , etx H = 



is / 
x 5 J i 



D'ou x = (1 2 5 0) T et = c T x = 2 - 2 = 0. Done, le vecteur 

x = (1 2 0) T realise le minimum de ("Po)- 

La troisieme contrainte du probleme (Vo) est : a;i + xz = 1 avec ii > et 
xz > 0. D'ou, a;3 = 1 — cci > ce qui donne xi < 1. Pour obtenir un probleme 
de dimension 2 equivalent a (Vo), il sufEt de remplacer la variable X3 par 1 — x\ 
et la contrainte x;i > par xi < 1, on obtient alors : 

min 2x\ — X2 + (1 — Ki) f min xi — X2 + 1 

xi +x 2 - (1 -xi) > -2 j 2a;i+x 2 >-l 

-xi +x 2 + 2(1 - xi) < 1 _ | -3xi+x 2 <-l 

a;i > 0, a;2 > 0, xi < 1 [ xi > 0,X2 > 0, x\ < 1 

5. On trouve la meme solution : 
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Reponse 4 \ 

On note par : 

xi = le nombrc de soldats produits chaque semaine , 
X2 = lc nombre de trains produits chaque semaine. 
D'apres les donnees, on a les contraintes suivantes : 

xi +x 2 < 80 
2xi + x 2 < 100 

X\ > 0, X2 > 

La fonction a maximiser est : 

(27 - 10 - U) Xl + (21 - 9 - 10)a;2 

Ce qui donne le programme lineaire suivant : 

max 3xi + 2x2 
xi + x 2 < 80 
2xi + x 2 < 100 

Xl > 0, X2 > 

Le programme est de dimension deux, done, on le resout graphiquement : 
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10 20 30 40 50 60 70 80 90 100 



le sommet (0 80) donne une valeur egale a:3x0 + 2x 80 = 160, 
le sommet (20 60) donne une valeur egale a:3x 20 + 2x 60 = 180, 
le sommet (50 0) donne une valeur egale a:3x 50 + 2x0 = 100, 
et le sommet (0 0) donne une valeur egale a:3x0 + 2x0 = 0. 
Le meilleur plan de production est done : 20 soldats et 60 trains. 

Reponse 5 

On note par : 

xi = la quantite du platl consommee par le client , 
X2 = la quantite du plat2 consommee par le client, 
x 3 = la quantite du plat3 consommee par le client, 
X4 = la quantite du plat4 consommee par le client, 
D'apres les donnees, on a les contraintes suivantes : 

' 400a;! + 200x 2 + 150x3 + 500x 4 > 500 

3xi + 2x2 > 6 
< 2xi + 2x 2 + 4a; 3 + 4:X 4 > 10 
2xi + 4x 2 + x 3 + 5x 4 > 8 

Xl > 0, X2 > 0,X3 > 0,X4 > 

La fonction a minimiser est : 

0.5xi + 0.2x2 + 0.3x3 + 0.4x4 
Ce qui donne le programme lineaire suivant : 
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min 0.5xi + 0.2x 2 + 0.3x 3 + 0.4x 4 

400xi + 200x 2 + 150x 3 + 500x 4 > 500 

3xi + 2x2 > 6 

2xi + 2x 2 + 4x 3 + 4x 4 > 10 

2xi + 4x 2 + x 3 + 5x 4 > 8 

Xl > 0, X2 > 0, X3 > 0, X4 > 

En utilisant l'applet Java de l'exempleinteractif, on trouve : 
Xi = 0, X2 = 3, X4 = 1, x 5 = 

Reponse 6 

1. On a rang(A) — m > 2, done il y a m colonnes de A qui sont lineairement 
independantes. Par consequent deux colonnes au moins sont lineairement 
independantes car m > 2. Done, il y a au moins une colonne qui n'est pas 
colineaire a d. 

2. a) II est clair que x > et on a : Ax = Aj + A n +\ = Aj + (d — Aj ) = d. D'ou, 
x € U. D'autre part, d'apres un theoreme du cours : 

x € U est un sommet de U -<=> {Aj /xj > 0} est une famille libre. 
Or {Aj/xj > 0} — {A jo , A n+ i} — {Aj ,d — Aj } avec Aj est non colineaire a 
d. Par consequent, d — Aj g est non colineaire a Aj (sinon d — Aj = aAj , cc 
qui donne d = (1 + a)Aj ce qui contredit le fait que d est non colineaire a Aj ). 
D'ou, {Aj/xj > 0} = {Aj ,d — Aj g } est une famille libre. 

b) La fonction lineaire a minimiser du probleme (V) verifie : 

Vx G U, z(x) = x n+ i > 

done, elle est minoree sur U et U / (d'apres 2)a)). D'apres un theoreme du 
cours, le minimum de (V) est atteint par un sommet de U. 

c) On a : £ = (£i £2 • • • £n £n+i) T est un point de Wqui realise le minimum de 
(V). D'ou, pour tout i = 1, • • • , n + 1 : 

> 

et 

d = A£ = 53 = X] + A n +itn+i =Ax + l n+ iC„+i 

avec x = (£1 £2 • • • Cn) T et la valeur du minimum de (V) est egale a z(£) = £ n +i- 
Alors, on a : 

Cn+i = O=^Ar = detx>O=>xe«=>W^0 
W ^ 3x = (xi x 2 • • • x n ) T G E n /Ax = d et x > 

On pose x = (xi X2 • • • x n 0) T , on obtient que x G et que 

z(x) = < «(x) Vx G W 

Ce qui prouve que x realise le minimum de (V) et que la valeur du minimum 
est egale a zero. D'ou, £ n +i = 0. 

d) £ = (£1 £2 • • • £n 0) T est un sommet de U , alors, la famille {Aj /£,• > 0} est 
libre. Or : 

{Aj/ti > 0} = {A,/^ > 0} 

car £n+i = et Aj = A, pour j ^ n + 1. Ce qui prouve que x = (£1 £ 2 • • • £n) T 
est un sommet de W car x G W et {Aj/£j > 0} est une famille libre. 
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